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The  authors  wish  to  acknowledge  several  persons  and  groups  who  contributed  to  the 
development  and  preparation  of  this  report.  Several  members  of  AFGWC  sections  SDDC  and 
SDDN,  especially  Major  Carol  Belt  and  Mr.  Ed  Carr,  carried  out  the  rather  arduous  special 
archiving  of  AFGWC  databases  needed  for  this  study.  Mr.  Don  Norquist  of  GL  derived  and 
provided  the  data  in  Figure  2  and  Table  4.  Mr.  Don  Aiken  of  GL  assisted  in  generating  some  of 
the  computer  graphics  as  well  as  managing  and  reading  the  many  data  tapes.  The  pursuit  of 
the  overall  study  benefited  from  the  many  discussions  with  and  advice  from  our  GL  colleagues 
Mr.  Don  Norquist,  Mr.  Don  Chisholm,  Dr.  C.-H.  Yang,  and  Dr.  H.  Stuart  Muench. 

The  new  GL  diagnostic  cloud  scheme  described  in  Section  3.1.3  was  first  presented  in 
December  1987  to  the  annual  meeting  of  the  GL  and  AWS  Numerical  Models  Technical 
Exchange  Group  (NUMTEG).  It  is  an  empirical  diagnostic  cloud  scheme  that  can  be  tuned  for 
any  given  forecast  model  and  can  optionally  include  a  real-time  update  feature,  made  possible 
by  the  real-time  operation  of  AFGWC's  automated  cloud  analysis  model,  known  as  the 
RTNEPH.  The  methodology  of  the  GL  cloud  scheme  evolved  from  steps  to  establish  a  critical 
value  of  relative  humidity,  RHc,  that  defined  a  range  of  model  forecast  humidity  (0  to  RHc) 
occurring  with  the  same  cumulative  frequency  as  zero  cloud  cover  in  the  verifying  RTNEPH 
cloud  analysis.  By  similarly  repeating  the  procedure  to  establish  critical  humidity  thresholds 
for  nonzero  RTNEPH  cloud  amounts,  quasi-continuous  empirical  curves  can  be  derived 
relating  the  full  range  of  forecast  humidity  to  the  full  range  of  fractional  cloud  cover.  The 
development  of  this  cloud  scheme  benefited  from  discussions  with  colleagues  Dr.  Ross 
Hoffman  of  AER,  Inc.  of  Cambridge,  MA  and  Lt  Dave  McDonald  of  AFGWC/SDDC. 

At  the  time  of  this  writing,  a  GL-funded  contract  effort  at  AFGWC  has  been  initiated  to 
perform  an  extensive  follow-on  study.  This  study  will  test  and  evaluate  the  new  GL  cloud 
scheme  in  the  AFGWC  operational  global  spectral  model  on  a  routine  dally  basis  over  a  period 
of  6-12  months.  In  support  of  this  contract  effort,  the  present  report  is  purposely  rather 
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detailed,  and  hence  probably  longer  than  preferred  by  most  readers.  To  this  latter  group  we 
recommend  Sections  1,  2.2,  3.1,  the  beginning  of  Section  3.2  (up  to  3.2.1)  and  Section  4. 

Finally,  the  latter  third  of  this  report  was  completed  after  the  lead  author  had  transferred 
from  GL  to  £  position  at  the  National  Meteorological  Center  beginning  In  November  1988.  As 
a  reserve  officer  with  AFGWC,  the  lead  author  thanks  AFGWC/SDD  for  permission  during  this 
period  to  complete  portions  of  this  report  while  on  reserve  duty. 
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Illustrations 


1.  Boundaries  of  the  NH  AFGWC  Octagon  Domain  and  the  U.S.  and  Europe  13 

Verification  Domains  Used  In  this  Study.  (Shading  key  at  lower  right 
applies  to  cloud  fields  in  Figures  30-47  and  shows  eight  shading 
Intensities,  including  blank,  for  eight  cloud  intervals). 

2.  (a)  Zonal  Time  Average  of  RH  for  the  FGGE  Level  II-B  Rawinsonde  16 

Observations  During  3-22  February  1979.  In  all  plots  of  zonal  average 
RH,  the  contour  interval  is  five,  (b)  For  the  Seven  Moist  Layers  of  the  GSM, 
the  Number  (X10)  of  Rawinsonde  Observations  of  RH  as  Summed  Over  10° 

Latitude  Bands  and  Ten  00Z  and  12Z  Time  Periods  During  15-20  February 
1979  (From  the  study  of  D.  Norquist12). 

3.  (a)  Zonal  Time  Average  of  RH  for  the  00Z  and  12Z  NMC  FGGE  III-A  18 

Objective  Analyses  During  3-22  February  1979.  (b)  Difference  Between 
Zonal  Time  Average  of  RH  In  Figures  3a  and  2a.  (In  all  difference  plots 
of  zonal  average  RH,  the  contour  interval  is  four  and  dashed  contours 
denote  negative  values). 

4.  (a)  Zonal  Time  Average  of  RH  for  the  NMC  Analysis  for  OOZ,  14  January  20 

1982.  (b)  Difference  Between  Zonal  Average  of  RH  in  Figures  4a  and  2a. 

5.  (a)  Zonal  Time  Average  of  RH  from  Eleven  OOZ  AFGWC  Analyses  During  22 

January  1985  (see  Table  1  for  dates),  (b)  Difference  Between  Zonal  Time 
Average  of  RH  in  Figures  5a  and  2a. 

6.  NH  AFGWC  30  kPa  RH  Analysis  for  OOZ,  17  January  1985.  (Shading  key  at  23 

lower  right  shows  five  shading  Intensities,  including  blank,  for  five  RH 
intervals). 

7.  Schematic  Diagram  of  the  AFGWC  Algorithm  Which  Infers  Initial  CPS  Moisture  24 

Values  for  the  AFGWC  5LAYER  Forecast  Model  from  the  AFGWC  RTNEPH 
Global  Cloud  Analysis.  (See  report  by  T.  Crum4  for  details). 

8.  Curves  of  the  AFGWC  CPS-cloud  Amount  Conversion  Scheme  Used  in  the  CPS-  25 

cloud  Conversion  Step  of  Figure  7. 
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9.  The  Critical  Relative  Humidity  (RH)  as  Derived  from  the  CPS-cloud  Curves  of 

Figure  8,  for  the  Lowest  Six  Mandatory  Pressures,  Assuming  U.S.  Standard 
Atmosphere  Temperatures  and  Applying  the  85  kPa  and  30  kPa  CPS-cloud 
Curves  Also  to  the  100  kPa  and  40  kPa  Levels,  Respectively. 

10.  (a)  Zonal  Average  of  RH  from  the  Experimental  Analysis  of  3DNEPH-inferred 

RH  Derived  by  Mitchell  and  Warburton6  for  00Z,  14  January  82.  (b) 

Difference  Between  Zonal  Average  of  RH  in  Figures  10a  and  2a. 

1 1 .  NH  Mean  RH  at  Six  Mandatory  Pressures  as  a  Function  of  Forecast  Length  in 

Four  GSM  Forecasts  Initialized  from  Four  Separate  Initial  Analyses  as 
Labeled. 

12.  Vertical  Profiles  of  the  NH  Mean  RH  at  the  (a)  0-hr  Initial  Time  and  (b)  96-hr 

Forecast  Time  for  the  Four  GSM  Forecast  Cases  in  Figure  1 1 .  Also  given 
as  a  reference  is  the  vertical  profile  of  NH  mean  RH  as  derived  from  the 
FGGE  II-B  zonal  average  RH  values  in  Figure  2. 

13.  (a)  Zonal  Average  of  RH  for  the  48-hr  GSM  Forecast  Initialized  from  the  NMC 

Analysis  of  00Z,  14  January  1982.  (b)  Difference  Between  the  Above  Average 
and  that  of  the  Corresponding  Initial  Analysis  (Figure  4a). 

14.  (a)  Zonal  Average  of  RH  for  the  48-hr  GSM  Forecast  Initialized  from  the  AFGWC 

Analysis  of  00Z,  17  January  1985.  (b)  Difference  Between  the  Above  Zonal 
Average  and  That  of  the  Corresponding  Initial  Analysis. 

15.  NH  NMC  RH  Analysis  for  00Z,  14  January  1982  at  Levels  (a)  100,  (b)  85,  (c)  70, 

(d)  50,  (e)  40  and  (f)  30  kPa.  (The  shading  key  in  Figure  6  applies  to  all 
shaded  displays  of  EH  fields). 

16.  NH  GSM  48-hr  RH  Forecast  at  Levels  (a)  100,  (b)  85,  (c)  70.  (d)  50,  (e)  40  and  (f)  30 

kPa  for  the  Forecast  Initialized  From  the  NMC  Analysis  of  Figure  15. 

17.  NH  AFGWC  RH  Analysis  for  00Z,  17  January  1985  at  Levels  (a)  100,  (b)  85,  (c)  70. 

(d)  50,  (e)  40  and  (f)  30  kPa.  (The  shading  key  in  Figure  6  applies  to  all  shaded 
displays  of  RH  fields). 

18.  NH  GSM  48-hr  RH  Forecast  at  Levels  (a)  100,  (b)  85,  (c)  70,  (d)  50,  (e)  40  and  (f)  30 

kPa  for  the  Forecast  Initialized  from  the  GSM  Forecast  Initialized  from  the 
AFGWC  Analysis  of  Figure  17. 

19.  Percent  Occurrence  Over  the  NH  Octagon  Domain  of  (a)  Moist  RH  Values  (RH>82) 

and  (b)  Dry  RH  Values  (RH<23)  at  the  50,  70,  and  85  kPa  Pressure  Levels  for 
the  0-,  48-,  and  96-hour  GSM  Forecasts  Initialized  from  the  NMC  Analysis 
of  00Z,  14  January  1982  (bottom)  and  the  AFGWC  Analysis  of  00Z,  17 
January  1985  (top). 

20.  Curves  of  the  ECMWF  Humidity-cloud  Conversion  Scheme  of  Eqs.  (l)-(2)  for 

Four  Mandatory  Pressures. 

21.  Critical  Relative  Humidity,  RHc,  as  a  Function  of  Pressure  for  the  Cloud  Amount 

Thresholds  of  (a)  0  Percent  and  (b)  30  Percent  for  the  ECMWF,  AFGWC,  NMC, 
and  GL  Humidity-cloud  Conversion  Schemes. 

22.  Curves  of  the  ECMWF,  AFGWC,  NMC,  and  GL  Humidity-cloud  Conversion  Schemes 

for  the  (a)  50  kPa  and  (b)  85  kPa  Mandatory  Pressures. 

23.  Curves  of  the  NMC  Humidity-cloud  Conversion  Scheme  of  Eq.  (4)  and  Table  6  for 

the  PBL,  Low  Middle,  and  High  Pressure  Regimes  of  Eq.  (3). 
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24.  Category  (top)  and  Cumulative  (bottom)  Frequency  of  Occurrence  (at  one  percent 

Intervals)  of  RTNEPH  Analyzed  Cloud  Amounts  (Part  a)  and  GSM  24-hour 
Forecast  RH  Amounts  (Part  b)  at  70  kPa  Over  the  NH  Octagon,  Valid  at  00Z, 

18  January  1985.  At  bottom,  the  one-to-one  mapping  of  the  RTNEPH  observed 
cloud  frequency  onto  the  GSM  forecast  RH  frequency  yields  one  of  the  GL 
humidity-cloud  conversion  curves  (solid)  In  Figure  25a. 

25.  Curves  of  the  GL  Humidity-cloud  Conversion  Scheme  Derived  for  the  GSM  24-hour  71 

Forecasts  at  (a)  70  kPa  and  (b)  30  kPa  by  Applying  the  Methodology  of  Figure 
24  to  Two  Separate  Samples  (and  their  aggregate)  of  the  Frequency  Distributions 
of  GSM  Forecast  RH  and  RTNEPH  Verifying  Cloud  Over  the  NH  Octagon. 

26.  Curves  of  the  GL  Humidity-cloud  Conversion  Scheme  as  a  Function  of  GSM  74 

Forecast  Length  for  (a)  70  kPa  and  (b)  30  kPa  Mandatory  Pressures.  (Curves 
based  on  the  aggregate  NH  octagon  GSM  RH/RTNEPH  cloud  frequency 
distributions  for  the  00Z,  17  and  24  January  1985  GSM  forecast  cases). 

27.  Curves  of  the  GL  Humidity-cloud  Conversion  Scheme  as  a  Function  of  Mandatory  76 

Pressure  for  the  GSM  0-hour  and  48-hour  Forecasts  (based  on  aggregate  sample 
as  in  Figure  26). 

28.  Curves  of  the  Four-layer  Version  of  the  GL  Humidity-cloud  Conversion  Scheme  (see  78 

text)  as  a  Function  of  Mandatory  Pressure  for  the  GSM  0-hour  and  48-hour 
Forecasts  (based  on  aggregate  sample  as  In  Figure  26). 

29.  The  ECMWF  Operational  NH  12Z  Analyses32  of  the  50  kPa  Height  Field  (top)  and  98 

Sea-level  Pressure  Field  (bottom)  for  (a)  17  January  and  (b)  19  January  1985. 

30.  The  RTNEPH  00Z  Analysis  of  NH  Octagon  Total  Cloud  Amount  for  (a)  17  January  99 

and  (b)  19  January  1985.  (The  shading  key  In  Figure  1  applies  to  all  shaded 
displays  of  cloud  amount). 

31.  The  RTNEPH  00Z  Analysis  of  NH  Octagon  Laver  Cloud  Amount  for  (a)  17  100-101 

January  1985  for  the  (a)  Gradient,  (b)  85,  (c)  70  and  (d)  30  kPa  Layers.  (The 
shading  key  In  Figure  1  applies  to  all  shaded  displays  of  cloud  amount). 

32.  The  RTNEPH  00Z  Analysis  of  NH  Octagon  Laver  Cloud  Amount  for  19  102-103 

January  1985  for  the  (a)  Gradient,  (b)  85,  (c)  70  and  (d)  30  kPa  Layers.  (The 
shading  key  In  Figure  1  applies  to  all  shaded  displays  of  cloud  amount). 

33.  The  5LAYER  45-hour  Forecast  of  NH  Octagon  Total  Cloud  Amount  Valid  at  OOZ,  105 

19  January  1985. 

34.  The  5LAYER  45-hour  Forecast  of  NH  Octagon  i^ysr  Cloud  Amount  Valid  at  OOZ,  106-107 

19  January  1985,  for  the  (a)  Gradient,  (b)  85,  (c)  70  and  (d)  30  kPa  Layers. 

35.  For  the  AFGWC  Cloud  Scheme,  the  GSM  (a)  0-hour  and  (b)  48-hour  Forecast  of  NH  112 

Total  Cloud  Amount  Valid  at  OOZ  on  17  and  19  January  1985,  Respectively. 

36.  For  the  AFGWC  Cloud  Scheme,  the  GSM  0-hour  Forecast  of  NH  Laver  Cloud  1 13-114 

Amount  Valid  at  OOZ  on  17  January  1985  for  Layers  Nominally  Centered 
at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa  Pressure  Levels. 

37.  For  the  AFGWC  Cloud  Scheme,  the  GSM  48-hour  Forecast  of  NH  Laver  Cloud  1 15-116 

Amount  Valid  at  OOZ  on  19  January  1985  for  Layers  Nominally  Centered 
at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa  Pressure  Levels. 

38.  Schematic  Illustration  of  the  Major  Static  Preprocessing  and  Postprocessing  Steps.  118 

which  when  Applied  (Along  with  the  Chosen  Cloud  Scheme  as  the  Last  Step) 
to  an  Input  Humidity  Analysis  on  Mandatory  Pressure  Levels  Yield  the 
So-called  "Synthesized"  GSM  0-hour  Cloud  Forecast. 
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39.  For  the  NMC  Cloud  Scheme,  the  GSM  (a)  O-hour  and  (b)  48-hour  Forecast  of  NH 

Total  Cloud  Amount  Valid  at  00Z  on  17  and  19  January  1985,  Respectively. 

40.  For  the  NMC  Cloud  Scheme,  the  GSM  (4d-hour  Forecast  of  NH  Laver  Cloud  122-123 

Amount  Valid  at  00Z  on  19  January  1985  for  Layers  Nominally  F entered 
r.t  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa  Pressure  Levels. 

41.  For  the  ECMWF  Cloud  Scheme,  the  GSM  (a)  0-hour  and  (b)  48-hour  Forecast  125 

of  NH  Total  Cloud  Amount  Valid  at  00Z  on  17  and  19  January  1985, 

Respectively. 

42.  For  the  ECMWF  Cloud  Scheme,  the  GSM  48-hour  Forecast  of  NH  Laver  Cloud  126-127 

Amount  Valid  at  00Z  on  19  January  1985  for  Layers  Nominally  Centered 
at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa  Pressure  Levels. 

43.  For  the  GL6  Cloud  Scheme,  the  GSM  (a)  0-houi  and  (b)  48-hour  Forecast  of  NH  129 

Total  Cloud  Amount  Valid  at  00Z  on  17  and  19  January  1985,  Respectively. 

44.  For  the  GL6  Cloud  Scheme,  the  GSM  0-hour  Forecast  of  NH  Laver  Cloud  130-131 

Amount  Valid  at  00Z  on  17  January  1985  for  Layers  Nominally  Centered 
at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa  Pressure  Levels. 

45.  For  the  GL6  Cloud  Scheme,  the  GSM  48-hour  Forecast  of  NH  Laver  Cloud  132-133 

Amount  Valid  at  00Z  on  19  January  1985  for  Layers  Nominally  Centered 
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Development  of  a  Cloud  Forecast 

GL  Baseline  Global 


Scheme  for  the 
Spectral  Model 


1.  INTRODUCTION 


Numerical  weather  prediction  groups  generally  incorporate  cloud  forecast  schemes  Into 
global  prediction  models  with  one  of  two  primary  objectives  In  mind.  The  first  and  most 
common  objective  Is  to  formulate  a  scheme  suitable  for  Interaction  with  the  global  model's 
parameterized  diabatlc  physics,  particularly  radiation.  These  schemes  are  ultimately  Judged 
by  their  Impact  on  the  parameterized  physics,  such  as  their  effect  on  the  earth's  radiation 
budget.1  In  this  context,  simple  cloud  schemes  are  often  utilized.  Some  examples  are 
prescribed  zonal  mean  cloud  amounts  at  every  grid  point  or  forecasts  of  ”on-ofF'  binary  cloud 
amounts  of  0  or  100  percent. 

The  second  objective,  which  is  the  foremost  objective  here,  is  to  forecast  cloud  cover  for  Its 
own  sake,  that  Is,  for  its  Importance  as  a  sensible  weather  element  to  end  users  of  numerical 


(Received  for  Publication  15  Dec.  1989) 

1  Geleyn,  J.-F.  (1981)  Some  diagnostics  of  the  cloud/ radiation  Interaction  in  the  ECMWF 
forecasting  model  ECMWF  Workshop  on  Radiation  and  Cloud-Radiation  Interaction  In 
Numerical  Modeling,  15-17  October  1980,  European  Centre  for  Medium  Range  Weather 
Forecasts,  Reading,  Berkshire,  U.K.,  135-162. 
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weather  prediction  (NWP)  products.*  Our  specific  goal  here  is  to  forecast  grid-area  values  of 
horizontal  fractional  cloud  cover  (both  total  and  layered)  in  the  0-4  day  range.  The  purpose  is 
to  develop  and  demonstrate  a  cloud  forecast  scheme  suitable  for  use  in  the  operational  global 
spectral  model  (GSM)  of  the  Advanced  Weather  Analysis  and  Prediction  System  (AWAPS)2  at 
Air  Force  Global  Weather  Central  (AFGWC).  The  cloud  cover  forecasts  must  be  sufficiently 
accurate  to  use  in  the  operation  of  numerous  Air  Force  missions  that  are  sensitive  to  clouds, 
particularly  their  obstruction  to  clear  lines  of  sight. 

Clearly,  the  two  cloud  forecast  objectives  cited  above  are  not  mutually  exclusive  and 
modeling  groups  consider  both  simultaneously.3  However,  the  state  of  the  art  in  modeling  the 
hydrological  cycle  and  other  diabatic  physics  in  global  models  is  not  sufficiently  advanced  to 
allow  pursuing  both  objectives  with  equal  priority.  Therefore,  in  developing  and  testing  cloud 
forecast  schemes,  modeling  groups  must  sometimes  adopt  methods  that  further  one  objective 
to  the  detriment  of  the  other.  Over  the  long  term,  as  modeling  advances  continue,  the 
approaches  to  achieve  these  two  cloud  forecast  objectives  should  converge.  In  the  meantime, 
the  near-term  emphasis  in  the  GL  effort  will  remain  accurate  forecasts  of  fractional  cloud 
cover. 

Before  the  implementation  of  the  AWAPS  at  AFGWC  in  October  1985,  the  AFGWC  model 
known  as  5 LAYER  was  the  only  moist  NWP  forecast  model  executed  by  AFGWC.  The  5 LAYER 
model  forecasts  temperature,  moisture,  and  precipitation  by  employing  quasi-Lagrangian 
numerical  methods,  in  which  3-hourly  3-D  air  parcel  displacements  or  trajectories  are 
computed  from  wind  velocity  forecasts  input  from  a  previously  executed  independent  NWP 
forecast  model.4  The  AWS  diy  PE  grid-point  model  served  as  the  independent  model  prior  to 
the  AWAPS  era.5  The  AWS  PE  model,  and  thus  also  5LAYER,  was  executed  only  out  to  48 
hours  and  only  on  extratropical  domains. 

With  the  advent  of  the  AWAPS  at  AFGWC,  the  concurrent  advent  of  a  moist  GSM  provided 
AFGWC  with  a  second  moist  model,  which  was  also  global  in  extent.  Since  the  beginning  of 
the  AWAPS,  AFGWC  has  continued  to  execute  the  5LAYER  model  as  its  primary  cloud  forecast 
model,  but  now  using  the  GSM  wind  velocity  forecasts  for  its  input  winds.  While  executing  the 
above  GSM/5LAYER  scenario,  AWS  and  AFGWC  have  naturally  considered  the  question  of 
whether  the  GSM,  being  a  moist  model,  could  be  utilized  directly  to  provide  cloud  forecasts,  thus 
eliminating  the  need  to  execute  5LAYER  as  a  follow-on  cloud  forecast  model.  The  present  study 


(Note:  all  acronyms  are  defined  in  List  of  Acronyms) 

2  Stobie,  J.G.  (1986)  AFGWC's  Advanced  Weather  Analysis  Prediction  System  (AWAPS), 
AWS/TN-86-001,  Air  Weather  Service,  Scott  AFB,  IL,  pp.  61. 

3  Slingo,  J.  and  Ritter,  B.  (1985)  Cloud  Prediction  in  the  ECMWF  model,  ECMWFTech.  Rep. 
No.,  46,  European  Centre  for  Medium  Range  Weather  Forecasts,  Reading,  Berkshire,  U.K., 
pp.  48. 

4  Crum,  T.D.  (1987)  AFGWC  Cloud  Forecast  Models,  AFGWC /TN-87-001,  Air  Force  Global 
Weather  Central,  Offutt  AFB,  NE,  pp.  66. 

5  Tarbell,  T.  and  Hoke,  J.  (1979)  The  AFGWC  Automated  Analysis / Forecast  Model  System, 
AFGWC -TN- 79/004,  Air  Force  Global  Weather  Central,  Offutt  AFB,  NE,  pp.  52. 
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addresses  this  question  In  the  context  of  the  resolution  and  physics  of  AFGWC’s  presently 
operational  GSM. 

The  GSM  cloud  forecast  experiments  presented  In  this  report  are  close  follow-ons  to  the 
GSM  experiments  described  in  the  paper  by  Mitchell  and  Warburton.6  We  urge  the  reader  to 
review  that  paper  (hereafter  denoted  as  MW)  as  a  forerunner  to  this  report.  The  GSM  cloud 
forecast  experiments  In  MW  were  performed  at  the  National  Meteorological  Center  (NMC)  using 
the  NMC  GSM. 

The  present  experiments  utilize  the  GL  global  spectral  model,  described  by  Brenner  et  al.7 
This  model  was  developed  at  GL  by  closely  following  the  NMC  spectral  model  developed  by 
J.  Sela,8  who  earlier  in  1980  with  other  NMC  staff  provided  the  NMC  GSM  codes  and 
substantial  consultation  to  GL.  In  the  GL  spectral  model,  the  formulation  for  the  adiabatic 
dynamics  of  the  NMC  model  was  utilized,  but  was  completely  recoded  by  GL  to  Increase  the 
code  modularity,  flexibility,  and  portability  desirable  in  R&D  applications.  The  GL  model  can 
be  configured  easily  to  run  In  numerous  horizontal  and  vertical  resolutions  and  on  various 
mainframe  computers  used  by  GL  and  GL  contractors  (see  Brenner  et  al.7  for  examples). 

Unlike  the  adiabatic  dynamics,  whose  formulation  is  now  fairly  fixed  and  has  remained 
essentially  unchanged  since  the  recoding  cited  above,  the  GL  GSM  parameterizations  for  the 
dlabatlc  physics  (for  example,  boundary  layer  fluxes,  precipitation,  radiation)  can  follow  any 
one  of  several  packages.  For  the  GSM  experiments  In  this  study,  we  used  GL's  "baseline" 
version  of  the  parameterized  physics.  This  version  utilizes  verbatim  the  parameterized 
physics  subroutines  presently  used  In  the  operational  GSM  in  AFGWC's  AWAPS.  The  latter 
operational  model  was  delivered  to  AFGWC  by  NMC  in  1983  and  represents  NMC's  operational 
GSM  as  it  stood  at  that  time. 

The  decision  to  utilize  the  baseline  physics  for  this  study  was  motivated  by  the  study 
objective  cited  earlier:  namely,  to  develop  and  demonstrate  a  cloud  forecast  scheme  suitable  for 
AFGWC's  currently  operational  GSM.  Thus  in  the  GSM  experiments  here,  which  were  executed 
on  the  Air  Force  Supercomputer  Center  -  Kirtland  (AFSCC-K)  CRAY-l/S  computer,  we  took  care 
to  parallel  the  configuration  of  operational  executions  of  AFGWC's  GSM.  Specifically,  in 
addition  to  using  the  baseline  (AFGWC)  parameterized  physics,  we  acquired  and  used  AFGWC's 
prescribed  GSM  surface  fields,  which  Include  terrain  height,  monthly  mean  sea-surface 
temperature,  and  surface-drag  coefficients.  Furthermore,  we  utilized  the  same  externally 
specified  control  parameters  (such  as  horizontal  diffusion  coefficients),  12-layer  vertical 
a-structure  (N12),  horizontal  rhomboidal  truncation  at  wavenumber  40  (R40),  and  nonlinear 
normal  mode  Initialization  package.  Lastly,  for  the  GSM  Initial  conditions,  we  used  AFGWC's 
operational  2.5°  latitude/longitude  global  objective  analysis  (uninitialized).  In  later 


6  Mitchell,  K.  and  Warburton,  J.  (1983)  A  comparison  of  cloud  forecasts  derived  from  the  NMC 
and  AFGWC  operational  hemisphere  forecasts  of  moisture,  Proceedings  of  AMS  Sixth 
Conference  or.  Numerical  Weather  Prediction,  6-9  June  1983,  Omaha,  NE,  pp.  66-73. 

7  Brenner,  S..  Yang,  C.-H.,  and  Mitchell,  K.  (1984)  The  AFGL  Global  Spectral  Model 
Expanded  Resolution  Baseline  Version,  AFGL-TR-84-0308,  ADA160370. 

8  Sela,  J.  (1980)  Spectral  modeling  at  the  National  Meteorological  Center,  Mon,  Wea,  Rev., 
108:1279-1292. 
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experiments,  we  also  used  several  NMC  global  2.5°  analyses  as  GSM  initial  states,  to  highlight 
several  serious  deficiencies  that  we  identified  in  the  AFGWC  global  moisture  analyses. 

The  possibility  of  future  AFGWC  application  of  the  GSM  cloud  forecast  algorithms  tested  In 
this  report  also  motivated  the  decision  to  evaluate  and  compare  the  experimental  GSM  cloud 
forecasts  against  the  operational  cloud  forecasts  from  the  AFGWC  5 LAYER  model.4  The  AFGWC 
global  objective  cloud  analysis  model,  known  as  RTNEPH  and  described  by  Kiess  and  Cox,9  Is 
used  by  AFGWC  to  routinely  verify  the  cloud  forecast  skill  of  the  5LAYER  model.  Therefore,  in 
this  study,  the  experimental  GSM  cloud  forecasts  as  well  as  the  operational  5 LAYER  cloud 
forecasts  for  the  same  valid  times  are  verified  in  parallel  against  the  RTNEPH  clord  analysis. 

It  is  appropriate  now  to  review  the  basic  components  of  the  "baseline"  physics  package  in 
the  GL  GSM  (or  present  AFGWC  GSM).  Following  that,  we  will  briefly  review  the  considerable 
efforts  at  GL  to  develop  substantially  improved  physical  parameterlzations  for  the  GSM.  In  so 
doing,  we  shall  clarify  the  context  of  the  present  study  within  the  larger  overall  GSM 
development  program  at  GL. 

The  baseline  physical  parameterlzations  Include  linear  horizontal  diffusion,  moist 
convective  adjustment  following  Kuo,10  dry  convective  adjustment,  and  a  simple  stable 
condensation  scheme  that  Immediately  condenses  out,  as  rain,  all  water  vapor  exceeding  a 
relative  humidify  saturation  criterion  of  1.0  (0.90  In  the  lowest  model  layer).  The  evaporation 
of  rainfall  is  Included  In  subsaturated  lower  levels.  The  planetary  boundary  layer  scheme  is 
very  simple  and  Includes  bulk  aerodynamic  surface  fluxes  of  momentum,  heat,  and  moisture 
in  the  first  model  layer  (the  surface  exchange  coefficients  are  not  stability  dependent).  Over 
land  the  surface  heat  and  moisture  fluxes  are  assumed  to  be  zero  and  over  oceans  they  are 
computed  using  a  specified  monthly  mean  surface  temperature.  Finally,  vertical  diffusion  Is 
Ignored  everywhere  above  the  first  model  layer  and  the  effects  of  radiation  are  entirely 
neglected. 

Since  1983  both  GL  and  NMC  have  pursued  comprehensive  improvements  to  the  baseline 
parameterized  physics.  The  GL  and  NMC  parameterization  efforts  are  proceeding  mostly 
Independently,  because  of  different  priorities  in  the  requirements  to  be  satisfied  by  the 
respective  GSMs  and  the  available  computer  resources.  However,  a  constant  exchange  of  model 
development  progress  takes  place  routinely  between  NMC,  GL,  and  AFGWC  through  periodic 
meetings  of  a  Numerical  Models  Technical  Exchange  Group. 

Within  the  GL  parameterization  program,  more  emphasis  is  being  placed  on  both  1)  the 
shorter  range  of  0-2  days  (Including  the  diurnal  cycle)  rather  than  the  medium  range  of  3-10 
days,  and  2)  the  hydrological  cycle,  to  provide  a  firm  foundation  for  forecasting  humidity, 
clouds,  and  precipitation.  The  latter  fields  are  notoriously  difficult  to  predict  In  global 
models,  owing  to  the  historically  coarse  resolution  of  such  models.  At  NMC  and  other  centers 
for  example,  emphasis  on  moisture-related  fields  is  foremost  In  regional  or  limited-area 
models,  but  somewhat  secondary  In  global  models.  GL  also  has  an  ongoing  development 


9  Kiess,  R.  and  Cox,  W.  (1988)  The  AFGWC  Automated  Real-Time  Cloud  Analysis  Model 
AFGWC/TN-88-001,  Air  Force  Global  Weather  Central,  Offutt  AFB,  NE.  pp.  82. 

10  Kuo,  H.L.  (1965)  On  formation  and  Intensification  of  tropical  cyclones  through  latent 
heat  release  by  cumulus  convection,  J.  Atmos.  Sd  22:40-63. 
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program  In  regional  and  mesoscale  modeling  (for  example,  Halberstam  et  al. 11  ).  However, 
because  the  USAF  operates  world-wide  missions  routinely,  it  may  not  always  be  feasible  to 
provide  explicit  regional  modeling  support  to  all  mission  areas  at  all  times.  This  emphasizes 
the  need  for  the  AFGWC  GSM  to  be  able  to  provide,  as  much  as  possible,  regional  support  on 
the  relatively  short  range  and  including  detailed  representation  of  the  moist  physical  pro¬ 
cesses  and  parameters  that  strongly  impact  Air  Force  operations. 

The  development  of  GL's  global  objective  analysis  and  4-D  data  assimilation  system 
(Norquist12)  provides  an  example  of  the  emphasis  on  humidity,  clouds,  and  precipitation  in 
GL's  global  modeling  program.  Through  contracts,  GL  is  pursuing  global  analysis/ assimilation 
techniques  tailored  for  moisture  and  alternative  global  data  sources  for  moisture  (Kaplan  et 
al.13).  Furthermore,  in  response  to  large  moist  biases  in  early  GL  assimilation  experiments, 

GL  developed  improved  methods  for  vertically  interpolating  moisture  in  4-D  assimilations 
(Mitchell  and  Yang14). 

With  a  similar  emphasis  on  moisture  and  clouds,  GL  has  pursued  its  R&D  in  GSM  physical 
parameterizatlons  through  contracts  with  three  universities,  one  each  in  the  areas  of 
radiation  (Liou  et  al.15),  deep  moist  convection  (Soong  et  al.16),  and  the  planetary  boundary 
layer  or  PBL  (Mahrt  et  al.17).  The  PBL  and  convection  initiatives  in  particular  have 
emphasized  moisture  and  hydrological  aspects.  In-house,  GL  is  testing  the  three  physical 
parameterizatlons  together  in  the  GL  GSM  and  a  first  evaluation  of  their  performance  and 
impacts  is  given  in  Yang  et  al.18.  Although  the  GSM  cloud  forecast  scheme  interacted  with  the 
physical  parameterizatlons  in  these  latter  tests,  the  first  evaluation  by  Yang  et  al.18  purposely 
did  not  examine  point-wise  cloud  forecast  accuracy.  Nevertheless,  as  the  new  physical 
parameterizatlons  are  further  extended  and  refined,  ultimately  a  chief  GL  criterion  for  judging 
the  new  parameterizatlons  will  be  the  impact  on  accurate  forecasts  of  fractional  cloud  cover 
(along  with  the  closely  associated  relative  humidity  distribution).  Thus,  the  present  effort 
here,  to  determine  the  cloud  forecast  skill  achievable  in  the  GL  GSM  using  the  original 


11  Halberstam,  I.M.  (1988)  Objective  Nephology,  Section  1:  Improved  Regional  Cloud  Forecast 
Model,  AFGL-TR-88-0109,  pp.  1-46,  ADA200500. 

12  Norquist,  D.  (1986)  Alternative  forms  of  moisture  information  in  4-D  data  assimilation, 
AFGL-TR-86-0194,  pp.  139,  ADA1 79792. 

13  Kaplan,  L.D.,  Hoffman,  R.N.,  Isaacs,  R.G.,  Rosen,  R.D.,  Salstein,  D.A.,  and  Wang,  W.-C.  (1983) 
Outlook  for  Improved  Numerical  Weather  Prediction  Using  Satellite  Data  with  a  Special  Emphasis 
on  the  Hydrological  Variables,  AFGDTR-83-0305,  pp.183,  ADA141233. 

14  Mitchell,  K.  and  Yang,  C.-H.  (1985)  A  Comparison  of  Moisture  Variables  in  the  Vertical 
Interpolations  of  a  4-D  Assimilation  System,  AFGL-TR-85-0090,  pp.  77,  ADA1 60464. 

15  Llou.  K.-N.,  Ou,  S.C.,  Kinne,  S.,  and  Koenig,  G.  (1984)  Radiation  Parameterization  Programs 
for  Use  in  General  Circulation  Models,  AFGL-TR-84-0217,  ADA148015. 

16  Soong,  S.-T.,  Ogura,  Y.,  and  Kau,  W.S.  (1985)  Study  of  Cumulus  Parameterization  in  a 
Global  Circulation  Model,  AFGL-TR-85-0160,  pp.  122,  ADA170137. 

17  Mahrt,  L.,  Pan,  H.-L.,  Paumier,  J.  and  Troen,  I.  (1984)  A  Boundary  Layer  Parameterization 
Scheme  for  a  General  Circulation  Model  AFGL-TR-84-0063,  pp.  169,  ADA144224. 

18  Yang,  C.-H.,  Norquist,  D.,  Yee,  S.Y.K.,  and  Mitchell,  K.  (1989)  Diagnostics  for  and 
Evaluations  of  New  Physical  Parameterization  Schemes  for  Global  NWP  Models,  GL-TR-89-0158. 
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baseline  physics,  will  provide  a  basis  for  comparing  future  cloud  forecast  experiments  using 
the  newer  physical  parameterizations. 

We  conclude  this  introduction  with  an  overview  of  the  sections  that  follow.  In  Section  2  we 
examine  the  GSM  humidity  forecasts.  To  do  this  in  a  complete  sense,  we  first  examine  the 
AFGWC  global  humidity  analyses  used  to  initialize  the  GSM  humidity  forecasts.  To  document 
several  substantial  weaknesses  in  the  AFGWC  humidity  analyses,  key  features  of  global 
relative  humidity  are  derived  from  global  radiosonde  observations  (RAOBs)  and  from  NMC 
operational  global  moisture  analyses.  The  above  emphasis  on  humidity  analyses  and 
forecasts  is  a  necessary  forerunner  to  understand  several  central  features  of  the  experimental 
GSM  cloud  forecasts  presented  later  in  Section  3. 

Section  3  begins  with  a  description  of  several  candidate  empirical  cloud  schemes  that 
relate  large-scale  humidity  fields  to  fractional  cloud  cover.  Three  of  these  humidity-to-cloud 
schemes  are  previously  known  schemes  obtained  from  operational  NWP  groups.  A  fourth 
empirical  cloud  scheme  is  highlighted  here  as  it  is  a  new  scheme  recently  developed  at  GL  as  a 
part  of  this  study.  This  new  GL  cloud  forecast  scheme  is  unique  in  that  it  is  the  first  known 
scheme  to  be  derived  objectively  from  coincident  samples  of  global  gridded  NWP  humidity 
forecasts  and  global  gridded  objective  cloud  analyses  (in  this  case  AFGWC  RTNEPH  cloud 
analyses).  Because  the  RTNEPH  is  a  daily  operational  cloud  analysis  model  at  AFGWC,  the 
method  of  deriving  the  empirical  humidity-to-cloud  relations  in  the  new  GL  cloud  scheme  can 
optionally  include  real-time  updates  of  the  relations  as  part  of  the  daily  AFGWC  analysis  and 
forecast  cycle.  As  a  by-product  of  such  a  real-time  update,  the  GL  cloud  scheme  can  adjust  to 
changes  in  resolution,  physics,  or  systematic  trends  in  either  the  NWP  model  humidity 
forecasts  or  the  RTNEPH  cloud  analyses.  Also  in  Section  3,  displays  of  the  experimental  GSM 
cloud  forecasts  are  presented,  as  obtained  by  employing  each  of  the  four  empirical  cloud 
schemes  in  conjunction  with  4-day  GSM  humidity  forecasts.  Finally,  Section  3  presents  the 
comparative  5LAYER  and  GSM  cloud  forecast  accuracy  statistics  obtained  by  grid-to-grid 
objective  verifications  of  the  various  cloud  forecasts  against  the  AFGWC  RTNEPH  cloud 
analysis.  Section  4  closes  out  this  report  with  conclusions  and  recommendations. 


2.  GLOBAL  HUMIDITY  ANALYSES  AND  FORECASTS 


The  various  cloud  forecast  schemes  documented  in  the  literature  fall  into  the  two  broad 
categories  of  1)  diagnostic  schemes  and  2)  prognostic  schemes.  In  prognostic  schemes  (for 
example,  Sundqvist19)  predictive  equations  are  formulated  and  integrated  to  predict  cloud 
water/ice  content  explicitly.  In  diagnostic  schemes  (for  example,  Slingo  and  Ritter3),  desired 
cloud  properties  are  diagnosed  empirically  or  statistically  from  other  model  output  fields,  such 
as  humidity,  precipitation  rate,  and  vertical  velocity.  In  the  present  study,  only  diagnostic 


19  Sundqvist,  H.  (1981)  Prediction  of  stratiform  clouds:  results  from  a  5-day  iorecast  with  a 
global  model,  Tellus,  33:242-253. 


cloud  schemes  will  be  considered  and  tested.  Although  prognostic  schemes  are  physically 
more  pleasing,  they  aic  in  the  early  stages  of  development,  more  difficult  to  verify,  and  more 
expensive  computationally.  More  Importantly,  they  typically  Introduce  numerous  arbltraiy 
parameters  to  relate  prognostic  cloud  quantities  (cloud  water)  to  desired  cloud  quantities 
(geometric  cloud  cc-  cij.  Thus,  there  may  be  as  much  or  more  empirical  or  diagnostic  character 
in  predictive  as  in  purely  diagnostic  cloud  schemes 

Virtually  all  diagnostic  cloud  schemes  depend  foremost  on  the  prediction  model’s  relative 
humidity  simulation  To  obtain  a  good  diagnostic  clot  id  forecast  then,  it  is  first  imperative  to 
achieve  realistic  humidity  forecasts,  regardless  of  the  chosen  diagnostic  cloud  scheme.  Thus, 
before  applying  several  different  diagnostic  cloud  schemes  in  Section  3,  we  first  examine  in 
this  section  the  crucial  characteristics  of  the  initial  global  humidity  analyses  and  follow-on 
GSM  humidity  foiecasts. 


2.1  Global  Humidity  Analyses 

At  the  beginning  of  this  study,  the  intent  was  to  utilize  only  AFGWC  global  objective 
analyses  to  initialize  and  execute  4-day  forecasts  from  GL's  GSM  (configured  to  be  parallel  to 
that  of  AFGWC).  As  the  study  progressed,  unrealistic  features  were  identified  in  the  AFGWC 
humidity  analyses.  These  poor  initial  humidity  features  and  tire  GSM's  large  short-range 
humidity  trends  (shown  later)  that  developed  in  response  to  these  initial  features  necessitated 
the  use  of  and  comparison  with  alternative  global  humidity  analyses. 

2.1.1  SOURCES  OF  THE  HUMIDITY  ANALYSES 

The  various  global  analyses  that  were  eventually  considered  ate-  summarized  In  Table  1 
and  discussed  further  in  this  subsection  The  collection  of  non  AFGWC  analyses  In  Table  1  is 
not  exhaustive,  but  intentionally  rather  limited  in  number  and  source  (NMC).  The  non- 
AFGWC  analyses  that  were  chosen  here  fot  consideration  were  largely  determined  by  (1)  the 
current  availability  of  these  analyses  at  GL  and  software  to  process  them,  and  (2)  the  adequacy 
of  these  analyses  in  providing  an  Inter comparison  set  sufficient  to  establish  the  major 
weaknesses  in  the  AFGWC  analyses. 


Table  1.  Attributes  of  the  Global  Analyses  Used  In  Study 
(All  analyses  performed  on  mandatory  pressure  levels). 


DESIGNATION 

SOURCE 

DATES 

ANALYZED 

VARIABLES 

HORIZONTAL 

RESOLUTION 

FGGE  III-A 

NMC 

3  22  FEB  79 

u,  v,  z,  Q+ 

2.5° 

NMC/OP 

NMC 

00Z,  14  Jan  82 

u,  v,  z,  RH 

24-wave 

NMC  /  3DNEPH 

NMC/GL 

00Z,  14  Jan  82 

u,  v,  z,  RH 

24-wave 

AFGWC 

AFGWC 

00Z,  Jan  85 
(10-12,  16-19, 
23-26) 

u,  v,  z *,  Q+ 

2.5° 

DESIGNATION 

FIRST-GUESS 

FORECAST  MODEL 

MASS/MOTION 

ANALYSIS  METHOD 

HUMIDITY 

ANALYSIS  METHOD 

FGGE  III-A 

9-layer  grid  point 

OI 

OI 

NMC/OP 

12-layer  spectral 

(R30) 

Hough 

Hough 

NMC  /  3DNEPH 

12-layer  spectral 

(R30) 

OI 

3DNEPH/CPS 

AFGWC 

9-layer  spectral 

(R20) 

OI 

Cressman 

The  primary  reference  set  of  non-AFGWC  analyses  In  Table  1  is  the  FGGE  III-A  analyses. 
These  are  the  optimal  interpolation  (OI)  analyses  (see  Bergman20)  that  were  products  of  the 
"final  update"  (that  is,  run  off  of  the  latest  data-cutoflf  time  and  using  the  global  forecast  model 
of  Stackpole21  as  the  first-guess)  of  the  NMC  operational  Global  Data  Assimilation 

20  Bergman,  K.H.  (1979)  Multivariate  analysis  of  temperature  and  winds  using  optimum 
interpolation,  Mon.  Wea.  Rev.,  107:1423-1444. 

21  Stackpole,  J.D.  (1978)  The  NMC  9-Layer  Global  Primitive  Equation  Model  on  a  Latitude- 
Longitude  Grid,  NMC  Office  Note  178,  National  Meteorological  Center,  Washington,  D.C. 
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System22  (GDAS),  as  configured  on  the  given  analysis  dates.  The  FGGE  designation  denotes 
First  GARP  Global  Experiment,  which  refers  to  the  internationally  sponsored  period  (in  1979) 
of  globally  augmented  conventional  and  satellite  weather  observing  systems.  The 
augmentation  of  observations  during  FGGE  was  especially  significant  during  the  Special 
Observing  Period  1  (5  January  1979  -  5  March  1979)  and  Period  2  (5  May  1979  -  5  July  1979) 
which  encompass  the  valid  times  of  the  FGGE  III-A  analyses  in  Table  1 .  Hence  the  NMC  GDAS 
during  these  periods  benefited  from  an  expanded  set  of  input  observations. 

The  secondary  set  of  non-AFGWC  analyses  in  Table  1  is  given  by  two  NMC-based  analyses 
for  00Z,  14  January  1982.  These  analyses  were  used  and  discussed  in  the  previous  study  of 
MW.6  The  analysis  designated  NMC /OP  corresponds  to  the  Hough-function  based  NMC 
operational  global  analysis  used  as  the  initial  state  for  the  NMC  global  forecast  cycle  (then 
called  the  "Large-Scale  Cycle").  This  analysis/forecast  cycle  was  executed  with  an  earlier  data 
cutoff  time  than  the  Ol-based  final  update  cycle  in  the  GDAS.  The  6-hour  GSM  forecast  from 
the  GDAS  OI  analysis  of  the  previous  6-hour  cycle  provided  the  first-guess  for  the  Hough 
analysis.  The  Hough  analysis  is  a  spectral  analysis  (rhomboidally  truncated  at  wavenumber 
24)  that  uses  spherical  harmonic  basis  functions  to  fit  available  data  by  means  of  a  global 
least  squares  technique.23  Further  details  and  limitations  of  the  Hough  analysis,  especially  as 
regards  the  humidity  fields,  are  given  in  MW.  With  respect  to  the  FGGE  III-A  analyses,  it  is 
Important  to  emphasize  that  the  NMC/OP  analysis  is  from  a  different  year  and  utilizes  1)  a 
different  analysis  method,  2)  an  earlier  data  cutoff  time,  and  3)  a  later  generation  GDAS24 
(including  the  newer  NMC  GSM8  as  the  first-guess  forecast  model).  Thus  the  NMC/OP  analysis 
will  serve  as  a  useful  cross  reference  to  validate  the  robustness  of  the  major  zonal-mean 
humidity  features  that  we  shall  first  infer  from  the  earlier  FGGE  Ili-A  analyses. 

The  analysis  designated  NMC/3DNEPH  is  identical  to  the  NMC/OP  analysis  for  the 
variables  u,v,  and  z.  In  NMC/3DNEPH  however,  the  operational  NMC  relative  humidity  (RH) 
analysis  of  NMC/OP  was  replaced  in  the  study  of  MW6  by  an  experimental  RH  analysis.  The 
latter  was  derived  by  applying  AFGWC's  empirical  humidity/cloud  conversion  scheme  (known 
as  the  CPS  scheme  detailed  later)  to  derive  humidity  profile  estimates  globally  from  the 
AFGWC  3DNEPH  global  cloud  analysis.  (The  predecessor  to  the  RTNEPH  was  the  3DNEPH,25 
which  was  executed  operationally  at  AFGWC  until  August  1983.)  The  ultimate  result  of  this 
application  was  a  conversion  of  global  3DNEPH  cloud  amounts  to  a  global  2.5°  RH  analysis  at 
the  lowest  six  mandatory  pressure  levels.  Further  details  on  the  derivation  of  this  humidity 
analysis  are  given  in  MW  as  well  as  below,  when  discussing  the  source  of  the  AFGWC  first- 
guess  field  for  the  AFGWC  humidity  analysis  in  Table  1 .  Both  derivations  are  closely 

22  McPherson,  R.D.,  Bergman,  K.H.,  Kistler,  RE.,  Rascli,  G.E.,  and  Gordon,  D.S.  (1979)  The 
NMC  operational  global  data  assimilation  system,  Mon.  Wea  Rev.,  107:1445-1461. 

23  Flattery,  T.  (1971)  Spectral  methods  for  global  analysts  and  forecasting,  Proc.  Sixth 
AWS  Technical  Exchange  Conference,  U.S.  Naval  Academy,  21-24  September  1970,  AWS-TR- 
242,  Air  Weather  Service,  Scott  AFB,  IL,  42-54. 

24  Dey,  C.H.  and  Morone,  L.L.  (1985)  Evolution  of  the  National  Meteorological  Center  global 
data  assimilation  system:  January  1982-December  1983.,  Mon.  Wea  Rev.,  113:304-318. 

28  Fye,  F.K.  (1978)  The  AFGWC  Automated  Cloud  Analysis  Model,  AFGWC-TM-78-002,  Air 
Force  Global  Weather  Central,  Offutt  AFB,  NE.  pp.  97. 
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analogous.  The  experimental  NMC/3DNEPH  humidity  analysis  serves  a  crucial  role  In 
Section  2.1.2,  as  it  illustrates  the  cause  of  a  pervasive  weakness  in  the  AFGWC  humidity 
analysis. 

Lastly  in  Table  1  are  the  analyses  designated  AFGWC.’.  These  analyses  were  operational 
products  of  the  AFGWC  global  4-D  assimilation  system  as  configured  during  the  January  1985 
period  denoted  for  them  in  Table  1.  Except  for  the  humidity  analyses,  the  production  AFGWC 
global  analyses  during  this  period  were  products  of  AFGWC's  global  Ol  analysis  model  known 
as  HIRAS2  and  the  AFGWC  GSM  first -guess  model  (obtained  from  NMC  in  1983  and  virtually 
identical  to  the  GSM  described  by  Sela8). 

The  AFGWC  global  4-D  data  assimilation  system  is  henceforth  referred  to  simply  as  the 
HIRAS/GSM  system.  The  HIRAS  OI  analysis  model  was  developed,  uoued,  and  tested  by 
AFGWC,  but  is  a  fairly  close  adaptation  of  the  NMC  Ol  analysis  model  described  by  Dey  and 
Morone.24  The  AFGWC  HIRAS/GSM  system  is  well  documented  in  the  report  by  Stobie.2 
Furthermore,  in  a  separate  report  Stobie  et  al.26  describe  (1)  how  the  AFGWC  OI  model  differs 
from  the  NMC  version  and  (2)  how  the  HIRAS/GSM  system  was  Implemented  in  two  major 
phases.  The  dates  of  the  AFGWC  analyses  In  Fable  1  are  during  the  period  of  Phase  I.  It  is 
important  for  this  study  to  next  describe  the  major  differences  between  Phase  I  and  Phase  II. 
especially  as  regards  the  humidity  fields.  The  period  of  Phase  t  extends  from  early  1984,  when 
the  HIRAS/GSM  system  was  first  implemented  until  October  1985,  when  AFGWC  acquired  Its 
CRAY  X-MP/12  supercomputer.  Thus  Phase  I  is  a  pre-CRAY  period  wherein  the  HIRAS/GSM 
system  operated  in  a  substantially  "scaled-down"  mode  as  dictated  bv  constrained  computer 
resources.  The  important  differences  between  Phase  1  and  Phase  11  are  listed  In  Table  2.  These 
differences  are  self  explanatory,  except  for  the  differences  regarding  the  humidity  analysis. 
This  latter  difference  is  explained  below  --  and  necessarily  in  some  detail. 


26  Stobie,  J.G.,  Lewis,  M.D.,  I-angford.  M.A.,  Tuell,  J.P.,  and  Carr.  E.L..  (1985)  The  use  of 
optimum  interpolation  at  AFGWC.  Preprints  Seventh  Conference  on  Numerical  Weather 
Prediction,  Montreal,  P.Q..  Panada,  Amec  Meteor.  Soc..  45 
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Table  2.  Main  Differences  Between  Phase  I  and  Phase  II  of 
AFGWC  AWAPS  Global  4-D  Assimilation  System. 


CHARACTERISTIC 

PHASE  I 

PHASE  II 

Max.  no.  of  obs 
analyzed  per 
grid  point 

8  (u,  v,  z,  Q) 

35  (u,  v,  z) 

20  (Q) 

Spectral  Resolution 
of 

Analysis 

Rhomboidal-24 

Rhomboidal-30 

Resolution  of 

GSM  6-hour 
first  guess 

Rhomboidal-20 

9  vertical  layers 

Rhomboldal-30 

12  vertical  layers 

Quality  Control  Steps 

Gross  check 

(No  buddy  check) 

Gross  check  and 

buddy  check 

Humidity  Analysis 

Not  stored 

into 

production  database 

Stored 

into 

production  database 

Resolution  of 

4-day  production 

GSM  forecast 

Rhomboidal-20 

9  vertical  layers 

Rhomboidal-40 

12  vertical  layers 

During  Phase  I,  the  global  moisture  analysis  produced  by  the  GSM/HIRAS  system  was  not 
stored  into  the  production  AFGWC  database.  That  is,  although  the  GSM/HIRAS  cycle  produced 
a  global  moisture  analysis  and  although  the  96-hour  GSM  moisture  forecast  was  initialized 
with  this  moisture  ana'vsis,  neither  the  HIRAS  moisture  analysis  nor  the  GSM  moisture 
forecast  were  stored  into  the  production  databases,  that  is,  those  databases  accessed  by  follow- 
on  applications,  display,  dissemination  and  archiving  programs.  This  decision  at  AFGWC  was 
based  on  the  following  two  main  reasons: 


1  -  Lack  of  evaluation  during  Phase  I  of  the  HIRAS/GSM  moisture  fields,  owing  to  the 

higher  priority  evaluation  requirements  for  the  wind,  height,  and  temperature 
fields. 

2  -  The  poor  Phase  I  spectral  resolution  of  both  the  6-hr  GSM  first-guess  and  the 

96-hour  GSM  forecast.  (Because  of  their  greater  spatial  variance,  upper-air 
humidity  fields  suffer  greater  deterioration  from  low  spectral  truncation  than 
upper-air  wind,  height,  or  temperature  fields.) 
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For  the  two  reasons  above,  the  AFGWC  production  2.5°  global  moisture  analyses  during 
Phase  I  continued  to  be  provided  by  the  AFGWC  moisture  assimilation  system  that  existed 
prior  to  the  HIRAS/GSM  era.  There  Is  no  accepted  name  for  this  system,  because  It  Involved 
an  array  of  models.  Nevertheless,  the  system  is  documented  in  the  report  by  Tarbell  and 
Hoke5  and  we  shall  refer  to  the  system  here  simply  as  the  pre-HIRAS  moisture  assimilation. 
The  latter  is  a  system  that  evolved  over  15-20  years  and  accordingly  Is  a  system  whose  domain 
expanded  In  stages  as  computer  resources  expanded.  Thus  the  pre-HIRAS  moisture 
assimilation  is  not  a  sti  alghtforward  unified  global  system;  but  rather  it  consists  of  separate 
surface  and  upper- air  analysis  models  over  separate  tropical  and  extra  tropical  domains. 

In  the  tropical  domain  (on  a  Mercator  projection  extending  from  30°  N  to  30°  S  latitude), 
the  first-guess  field  for  the  moisture  analysis  is  merely  persistence  of  the  previous  analysis 
(hence  the  tropical  first-guess  moisture  field  is  substantially  dependent  on  a  tropical  moisture 
climatology,  which  was  the  first-guess  for  the  first  tropical  moisture  analysis).  In  the  NH  and 
SH  extratropical  domains  (on  a  polar  stereographic  projection),  the  first-guess  field  for  the 
moisture  analysis  is  the  3-hour  moisture  forecast  of  AFGWC's  NH  and  SH  5LAYER  model.4 
Figure  1  shows  the  boundaries  of  the  NH  octagon  domain,  mirrored  by  a  companion  SH 
octagon  domain.  Given  their  respective  first-guess  moisture  fields,  the  separate  tropical  and 
extratropical  moisture  analysis  models  (known  as  TROPUA  and  MULTAN  respectively5) 
applied  a  standard  Cressman  analysis  scheme27  to  all  available  RAOBs  of  humidity.  During 
Phase  I,  a  utility  program  called  DATFCH  merged  the  separate  tropical  and  extratropical 
humidity  analyses.  The  latter  program  interpolated  the  several  humidity  analyses  to  a  single 
global  2.5°  latitude/longitude  grid  at  the  lowest  six  mandatory  pressure  levels,  defaulting  to  the 
extratropical  humidity  analysis  in  areas  of  domain  overlap  (except  along  a  fairly  narrow 
latitude  belt  inside  the  octagon  boundary  where  some  blending  with  the  tropical  analysis 
occurred).  The  essential  details  of  the  above  discussion  are  summarized  In  Table  3. 


27  Cressman,  G.P.  (1959)  An  operational  objective  analysis  system,  Mon.  Wea.  Rev., 
87:367-374. 
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Figure  1 .  Boundaries  of  the  NH  AFGWC  Octagon  Domain  and  the 
U.S.  and  Europe  Verification  Domains  Used  In  this  Study.  (Shading 
key  at  lower  right  applies  to  cloud  fields  In  Figures  30-47  and  shows 
eight  shading  intensities,  including  blank,  for  eight  cloud  intervals). 
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Table  3.  Outline  of  AFGWC  Pre-HIRAS  Operational 
Global  Moisture  Analysis. 


DOMAIN 

RESOLUTION 

AFGWC 

NAME 

ANALYSIS 

METHOD 

SOURCE  OF 

FIRST-GUESS 

SOURCE 

OFOBS 

N.H.  Octagon 

whole-mesh 

MULTAN 

Cressman 

3-hour 

5LAYER  fcst 

RAOBs 

S.H.  Octagon 

whole-mesh 

MULTAN 

Cressman 

3-hour 

5LAYER  fcst 

RAOBs 

Tropics 

30S-30N 

whole-mesh 

TROPUA 

Cressman 

Persistence 

and  Climo. 

RAOBs 

The  three  separate  analysis  domains  above  are  merged  and  interpolated  to  a  global  2.5° 
lat/lon  grid. 

In  addition  to  the  January  1985  cases  for  the  AFGWC  analyses  in  Table  1,  AFGWC  archived 
April.  July,  and  October  analyses  during  1984  for  GL  use  in  this  study.  However,  until 
corrected  in  late  October  1984,  a  logic  error  in  the  DATFCH  program  cited  above  inadvertently 
resulted  in  zero  humidity  values  being  stored  at  all  tropical  points  (that  is,  outside  the  NH  and 
SH  octagon  domains)  in  the  Phase  I  global  humidity  analyses.  Hence  the  Phase  I  analyses 
archived  before  November  1984  were  unusable  in  this  study  and  are  not  listed  in  Table  1. 

Despite  yet  other  weaknesses  in  the  January  1985  Phase  I  AFGWC  global  moisture  analyses 
to  be  documented  in  Section  2,1.2.  Section  3.2  demonstrates  that  beyond  the  12-  to  24-hour 
GSM  forecast  range,  competitive  and  skillful  global  cloud  forecasts  are  obtained  from  the  GSM 
moisture  forecasts  initialized  from  these  analyses.  To  fully  measure  the  potential  skill  in 
GSM-based  cloud  forecasts  from  AFGWC's  present  AWAPS  system,  one  obvious  consideration 
was  to  extend  the  present  study  to  the  use  of  Phase  II  AFGWC  global  analyses,  following 
October  1985  when  the  higher  resolution  AWAPS  assimilation  cycle  was  Implemented  and  the 
GSM/HIRAS  global  moisture  analysis  was  stored  and  archived  as  the  production  AFGWC 
global  moisture  analysis.  However,  the  GL  study  of  Mitchell  and  Yang14  indicated  the  strong 
likelihood  of  substantial  systematic  moist  biases  in  the  Phase  II  GSM/HIRAS  moisture 
assimilation  due  to  poor  methods  of  vertically  interpolating  analyzed  specific  humidity  (full 
values,  not  residuals)  from  the  pressure  surfaces  of  the  HIRAS  analysis  to  the  o-coordinate 
surfaces  of  the  GSM  forecast  model  and  vice  versa.  This  foreseen  bias  was  fully  substantiated 
by  subsequent  AFGWC  evaluation  of  the  Phase  II  moisture  analyses,  but  the  problem  was  not 
corrected  until  March  1988.  All  experiments  for  the  present  study  had  to  be  completed  before 
the  latter  date.  A  second  follow-on  GSM  cloud  forecast  study  using  Phase  II  AFGWC  moisture 
analyses  subsequent  to  March  1988  was  initiated  in  July  1989. 

There  is  another  crucial  reason  why  the  present  GSM  cloud  forecast  study  was  carried 
through  with  AFGWC  Phase  I  global  moisture  analyses.  Unlike  Phase  II,  the  Phase  I  moisture 


analyses  implicitly  utilized  RTNEPH  global  cloud  analyses  as  a  source  of  estimated  humidity 
observations.  This  was  by  virtue  of  the  previously  cited  use  of  the  5IAYER  model's  moisture 
forecast  as  the  Phase  1  first-guess  humidity  field  tn  the  extratropics.  As  detailed  further  in  the 
next  subsection,  5LAYER  initial  humidity  fields  are  almost  entirely  inferred  from  the  RTNEPH 
cloud  analysis.  In  contrast.  Phase  II  global  moisture  analyses  utilize  only  RAOBs  as  a  source 
of  humidity  observations.  Therefore,  to  maximize  the  potential  GSM  cloud  forecast  skill  in 
comparison  to  that  of  5LAYER,  especially  in  the  short  range,  we  felt  it  imperative  that  the 
GSM  moisture  forecast  be  initialized,  like  SLAYER,  from  a  moisture  analysis/assimilation 
system  that  utilized  the  RTNEPH  cloud  analysis  (and  hence,  indirectly,  satellite  cloud  imagery 
and  surface  cloud  reports)  as  a  data  source.  This  remains  a  real  possibility,  but  not  (as  we 
shall  show)  using  the  faulty  eloud-to-humidity  inference  methods  employed  in  Phase  I.  In  an 
important  related  study,  Norquist12'28  demonstrated  that  the  quality  of  global  humidity 
analyses  can  be  improved  by  assimilating  humidity  estimates  inferred  from  RTNEPH  cloud 
analyses,  provided  an  appropriate  cloud-to-humidity  inference  scheme  is  employed. 

2.1.2  CHARACTERISTICS  OF  THE  HUMIDITY  ANALYSES 

The  previous  section  established  that  the  humidity  fields  in  the  various  global  analyses  in 
Table  1  were  derived  with  several  diflerent  first-guess  forecast  models,  analysis  methods,  and 
data  sources.  Not  suxprislngly  then,  tills  section  will  demonstrate  that  the  gross  features  of 
these  humidity  fields  have  some  substantial  differences  from  each  other. 

The  differences  are  best  understood  by  grouping  the  analyses  in  Table  1  into  two  sets  --  Set 
A  consisting  of  the  NMC  operational  analyses  (FGGE  III  A  and  NMC/OP)  and  Set  B  consisting 
of  the  analyses  (AFGWC  and  NMC/3DNEPH)  that  directly  or  indirectly  utilize  the  3DNEPH  or 
RTNEPH  cloud  analysis  database.  Between  the  two  sets,  the  general  characteristics  of  the 
humidity  fields  show  substantial  differences,  but  within  each  set  the  characteristics  are  quite 
consistent.  The  emphasis  in  these  comparisons  will  be  on  zonal  mean  humidity  fields.  In 
particular,  the  zonal  mean  features  of  the  humidify  analyses  In  Set  A  will  be  found  to  agree 
well  with  observations,  despite  the  fact  that  the  FGGE  III  A  and  NMC/OP  analyses  utilized  a 
different  analysis  method  and  first-guess  model.  The  zonal  mean  features  of  the  humidity 
analyses  in  Set  B  show  substantial  departures  from  the  zonal  means  derived  from 
observations  and  from  Set  A. 

The  background  truth  fields  for  this  section  arc  provided  by  the  FGGE  II-B  RAOBs  gathered 
during  the  period  of  3-22  February  1979.  which  (alls  within  FGGE  Special  Observing  Period  1. 
The  zonally  and  temporally  averaged  relative  humidity  (RH)  field  obtained  from  this  nearly 
three-week  dataset  is  given  in  Figure  2.a.  The  corresponding  depiction  of  RH  observation 
frequency  during  five  days  within  that  period  is  given  in  Figure  2,b.  Figures  2a  and  2b  were 
taken  from  Figures  16  and  17  In  the  study  by  Norquist.12  in  that  study,  to  obtain  the  zonal 
average  in  Figure  2a.  all  RAOBs  were  n»sl  grouped  according  to  nearly  equal -area  grid  cells  of 
10  degrees  width  In  latitude,  then  the  mean  RH  value  fm  each  grid  cell  was  calculated,  and 

28  Norquist,  D.  (1988)  Alternative  forms  of  humidity  information  in  global  data  assimilation, 
Mon.  Wea.  Rev..  116  452  471 
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lastly  the  zonal  mean  was  computed  from  the  mean  grid-cell  values.  In  this  fashion,  areas  of 
high  observation  density  were  not  given  unduly  large  weight. 


Z0NAL  MEAN  (RH)  FGGE  LEVEL  II-B  3-22  FEB  1979 


N0.  RH  RA0BS  (X10)  2/15. 12Z-20. 00Z/79 


Figure  2.  (a)  Zonal  Time  Average  of  RH  for  the  FGGE  Level  n-B 
Ra  win  sonde  Observations  During  3-22  February  1979.  In  all  plots 
of  zonal  average  RH,  the  contour  interval  is  five,  (b)  For  the 
Seven  Moist  Layers  of  the  GSM,  the  Number  (XI 0)  of  Rawinsonde 
Observations  of  RH  as  Summed  Over  10°  Latitude  Bands  and  Ten 
00Z  and  12Z  Time  Periods  During  15-20  February  1979  (From  the 
study  of  D.  Norquist12). 
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The  distinctive  features  of  Figure  2a  are  1 )  the  rather  low  mean  humidity  in  the  subtropical 
descending  branches  of  the  Hadley  circulation,  2)  the  rather  high  mean  humidity  in  the  region 
of  mean  ascending  motion  around  the  Inter-Tropical  Convergence  Zone  (TTCZ),  and  3)  the  high 
mean  humidities  poleward  of  40°  latitude  In  both  hemispheres,  particularly  in  the  NH  (the 
hlgh-latitude  low-level  anomalies  of  the  SH  are  suspect  because  of  the  scarcity  of  observations 
there).  To  compare  with  Figure  2a,  we  first  consider  the  corresponding  FGGE  III-A  objective 
analyses,  which  are  listed  first  in  Table  1 .  Figure  3a  is  the  zonal  and  temporal  average  of  the 
dally  00Z  and  12Z  NMC  operational  global  objective  analyses  (FGGE  m-A)  for  the  same  period 
as  Figure  2a.  Figure  3a  agrees  well  with  Figure  2a,  showing  the  same  three  general  features 
cited  above.  Figure  3b  provides  a  plot  of  the  differences  between  Figures  3a  and  2a.  Except  for 
the  differences  poleward  of  latitude  50°  S  (where  Figure  2a  is  suspect),  the  zonal  mean  RH 
values  in  Figure  3a  generally  agree  with  those  In  Figure  2a  to  within  8  percent  RH. 
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PRESSURE  LEVEL  (CB)  PRESSURE  LEVEL  tCB ) 


Z0NAL  MEAN  (RH)  FGGE  LEVEL  III-fl  3-22  FEB  1979 


Figure  3.  (a)  Zonal  Time  Average  of  KH  for  llie  d()Z  and  12Z 
NMC  FGGE  III -A  Objective  Analyses  During  3-22  February  1979. 

(b)  Difference  Between  Zonal  Time  Average  of  RH  in  Figures  3a  and 
2a.  (In  all  difference  plots  of  zonal  average  RH,  the  contour  interval 
is  four  and  dashed  contours  denote  negative  values). 
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Considering  next  the  second  source  of  objective  analyses  in  Table  1 ,  we  turn  to  Figure  4a, 
which  provides  the  zonal  average  of  the  NMC  operational  global  analysis  for  00Z  on  14 
January  1982  (this  is  the  case  highlighted  in  MW6).  The  similarity  between  Figure  4a  and 
Figure  3a  is  striking,  considering  that  Figure  4a  is  from  a  different  year  and  from  an  earlier 
winter  month.  Moreover,  the  NMC  analysis  represented  in  Figure  4a  utilized  a  different 
analysis  method  and  a  later  generation  NMC  global  data  assimilation  system  (see 
Section  2.1.1)  that  included  a  different  global  forecast  model  for  the  first  guess.  The  similarity 
between  Figures  3a  and  4a  is  further  illustrated  by  comparing  Figures  3b  and  4b,  where 
Figure  4b  is  the  plot  of  the  difference  between  Figure  4a  and  Figure  2a.  The  key'  purpose  of 
Figure  4a  is  to  demonstrate  the  interannual  robustness  of  the  pattern  of  NH  winter  season 
zonal  mean  RH  depicted  in  Figures  2a  or  3a. 
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Figures  2-4  provide  a  sufficient  reference  to  compare  against  the  zonal  mean  RH  fields 
obtained  from  the  AFGWC  analyses  In  Table  1.  For  purposes  of  this  comparison.  Figure  5a 
provides  the  zonal  and  temporal  average  of  the  eleven  00Z  January  1985  AFGWC  RH  analyses 
cited  In  Table  1.  The  plot  of  the  difference  between  Figure  5a  and  Figure  2a  Is  given  in 
Figure  5b.  Figure  5  shows  large  departures  from  Figures  2-4.  The  most  outstanding  feature  of 
Figure  5  is  the  very  dry  and  erroneous  character  of  the  AFGWC  analysis  In  the  upper-half  of 
the  troposphere  In  the  extratropics  of  both  hemispheres.  The  discussion  below  establishes  and 
addresses  the  source  of  this  error. 
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Figure  5.  (a)  Zonal  Time  Average  of  RH  from  Eleven  OOZ  AFGWC 
Analyses  During  Januaiy  1985  (see  Table  1  for  dates),  (b)  Difference 
Between  Zonal  Time  Average  of  RH  in  Figures  5a  and  2a. 
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Figure  6  provides  a  shaded  display  of  the  NH  30  kPa  AFGWC  RH  analysis  for  00Z  on  1 7 
January  1985.  Superimposed  on  Figure  6  Is  the  AFGWC  octagon  domain  boundary  of  Figure  1. 
It  is  evident  that  the  widespread  region  of  uniformly  low  RH  values  in  Figure  6  occurs  over 
that  region  where  the  exclusive  source  of  the  AFGWC  analysis  is  the  extratropical  octagon 
humidity  analysis  that  utilizes  the  humidity  forecast  of  the  AFGWC  5LA.YER  model  as  the 
first-guess  field.  (See  Table  3  and  the  earlier  discussion  of  the  AFGWC  DATFCH  program 
toward  the  end  of  Section  2.1.1.)  It  is  also  noteworthy  in  Figure  5a  that  the  extratropical 
upper-level  dry  bias  is  more  extreme  in  the  SH  (versus  the  NH)  where  Figure  2b  shows 
considerably  fewer  RAOBs  are  available  to  correct  the  5LAYER  forecast  first-guess  field.  These 
facts  point  to  the  5LAYER  first-guess  field  as  the  source  of  the  upper-level  extratropical  dry 
bias  in  Figure  5.  Additional  facts  supporting  this  supposition  follow. 


Figure  6.  NH  AFGWC  30  kPa  RH  Analysis  for  00Z,  17  January  1985. 
(Shading  key  at  lower  right  shows  five  shading  intensive,  including 
blank,  for  five  EH  intervals). 
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As  shown  in  Table  3,  the  AFGWC  5LAYER  first-guess  humidity  forecast  for  the 
extratropical  AFGWC  humidity  analysis  (MULTAN)  is  a  three-hour  short-range  forecast. 

Hence  the  character  of  this  first-guess  humidity  forecast  is  very  similar  to  the  initial  O-hour 
5LAYER  humidity  analysis.  The  5LAYER  initial  or  O-hour  humidity  analysis  is  unique  in  that 
it  utilizes  the  AFGWC  RTNEPH  3-D  global  cloud  analysis  as  its  primary  humidity  "data" 
source.  Figure  7  provides  a  summary  schematic  of  the  procedure  as  applied  in  January  1985 
(and  as  still  applied  at  the  time  of  this  writing  except  for  some  changes  in  the  details  of  the 
vertical  compaction4).  First,  the  high  resolution  3-D  RTNEPH  cloud  analysis  is  vertically 
compressed  and  horizontally  averaged  to  yield  a  "compacted"  3-D  cloud  analysis  on  the  coarser 
resolution  grid  of  the  5LAYER  forecast  model.  Secondly,  the  empirical  cloud-to-humidity 
conversion  curves  in  Figure  8,  known  as  the  "CPS"  curves,  are  used  to  convert  the  compacted 
3-D  RTNEPH  cloud  amounts  to  inferred  pseudo  or  surrogate  initial  humidity  values.  At  grid 
points  where  the  compacted  3-D  RTNEPH  cloud  amounts  are  strictly  zero  (clear),  the  humidity 
analysis  defaults  to  the  drier  of  1)  the  zero-cloud  threshold  value  or  critical  value  from 
Figure  8  or  2)  the  humidity  value  from  the  previous  cycle’s  conventional  MULTAN  humidity 
analysis  (Table  3). 


CPS  INITIALIZATION  IN  CLOUDY  AREAS 

3DNEPH  CLOUD 

25  NM 

1 5  LAYERS 


5LAYER  CLOUD 


5 LAYER  CPS 


100  NM 
5  LAYERS 


100  NM 
5  LAYERS 


COMPACT  CONVERT 


Figure  7.  Schematic  Diagram  of  the  AFGWC  Algorithm  Which 
Infers  Initial  CPS  Moisture  Values  for  the  AFGWC  5LAYER 
Forecast  Model  from  the  AFGWC  RTNEPH  Global  Cloud  Analysis. 
(See  report  by  T.  Crum4  for  details) . 
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Figure  8.  Curves  of  the  AFGWC  CPS-cloud  Amount  Conversion 
Scheme  Used  in  the  CPS-cloud  Conversion  Step  of  Figure  7. 


The  humidity-to-cloud  curves  in  Figure  8  are  defined  in  terms  of  the  moisture  variable 
called  condensation  pressure  spread  {CPS}.  For  an  air  parcel  having  a  given  temperature  and 
dewpoint.  CPS  is  defined  as  the  amount  of  upward  dry  adiabatic  vertical  displacement  in 
millibars  required  for  the  parcel  to  reach  saturation.  Thus  a  zero  CPS  value  denotes 
saturation  (RH  =  100  percent)  and  increasingly  larger  values  of  CPS  denote  Increasingly  drier 
air  (analogous  to  dewpoint  depression).  The  AFGWC  CPS  curves  were  originally  derived  by  Air 
Weather  Service  in  the  early  1960's  from  a  coincident  sample  of  surface  cloud  reports  (note: 
not  3DNEPH  or  KTNEPH  cloud  amounts)  and  RAOBs.29  Unfortunately,  no  information 
remains  on  the  nature  of  this  coincident  cloud/  humidity  observation  sample  such  as  the  time 
period,  the  region,  or  scatter  plots. 


29  Edson.  H.  (1965)  Numerical  Cloud  and  Icing  Forecasts,  Scientific  Services  Technical  Note  13. 
3rd  Weather  Wing,  Offutt  AFB,  NE,  pp.  44. 
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Given  a  CPS  value  and  a  coincident  temperature  value,  any  of  several  other  common 
moisture  values,  such  as  RH.  may  be  derived.  Using  U.S.  Standard  Atmosphere  temperatures, 
the  critical  CPS  values  for  zero  cloud  amount  in  Figure  8  have  been  converted  to  their 
corresponding  critical  RH  values  (RHc)  and  plotted  for  the  lowest  six  mandatory  pressure  levels 
in  Figure  9.  (To  obtain  the  critical  RH  values  for  the  40  and  100  kPa  levels  in  Figure  9,  the  30 
and  85  kPa  critical  CPS  values  from  Figure  8  were  assumed  respectively,  as  the  AFGWC  CPS 
curves  do  not  include  separate  curves  for  the  40  and  100  kPa  levels.)  The  implied  critical  RH 
values  in  Figure  9  at  and  above  50  kPa  are  unexpectedly  low. 


0  50  100 

RHc  (Pet) 


Figure  9.  The  Critical  Relative  Humidity  (RH)  as  Derived  from 
the  CPS-cloud  Curves  of  Figure  8,  for  the  Lowest  Six  Mandatory 
Pressures,  Assuming  U.S.  Standard  Atmosphere  Temperatures 
and  Applying  the  85  kPa  and  30  kPa  CPS-cloud  Curves  Also  to 
the  100  kPa  and  40  kPa  Levels,  Respectively. 


On  a  separate  but  relevant  matter,  the  studies  by  Hughes  and  Henderson-Sellers30,31  of  the 
FGGE  year  climatology  of  the  January  and  July  1979  3DNEPH  cloud  databases  show  that  the 
mean  amounts  of  high  cloud  in  the  3DNEPH  database  are  lower  than  several  other 
independent  cloud  climatologies.  (However,  the  3DNEPH  low,  middle,  and  total  cloud  amounts 


30  Hughes,  W. A.  and  Henderson-Sellers,  A.  (1985'  Global  3-D  nephanalysis  of  total  cloud 
amount:  climatology  for  1979,  J.  Climate  Appl  Meteor.,  24:669-686. 

31  Henderson-Sellers,  A.  (1986)  Layer  Cloud  Amounts  for  January  and  July  1979  from 
3D-Nephanalysis,  J.  Climate  Appl.  Meteor.,  25:1  18- 132. 
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did  compare  favorably  with  other  climatologies,  except  around  the  two  polar  regions  and  the 
summer  Asian  monsoon  region.) 

When  the  CPS  curves  are  applied  to  3DNEPH  or  RTNEPH  cloud  amounts  to  infer  humidity 
estimates,  the  above  two  independent  factors  of  1)  surprisingly  low  equivalent  critical  RH 
values  from  the  upper -level  CPS  curves  and  2)  low  mean  amounts  of  high  cloud  In  the  3DNEPH 
or  RTNEPH  database  both  contribute  to  yield  erroneously  low  inferred  humidity  estimates. 
Most  importantly,  the  conclusion  that  the  3DNEPH/RTNEPH-inferred  humidity  estimates  at 
upper  tropospheric  levels  are  erroneously  low  is  substantiated  by  the  independent  results  of 
Table  4,  obtained  from  the  collocation  study  of  Norquist. 12,28  Table  4  shows  a  substantially 
dry  bias  in  humidity  estimates  Inferred  from  applying  the  upper-level  CPS  curves  of  Figure  8 
to  the  upper-level  February  1979  3DNEPH  cloud  amounts,  relative  to  measured  RH  from 
collocated  RAOBs  for  the  same  February  1979  winter  period  as  Figures  2  and  3. 


Table  4.  Statistics  of  Differences  Between  Relative  Humidity 
Inferred  from  Cloud  Amount  and  that  Measured  with  Radio¬ 
sondes  (In  %)  (from  Norquist1228). 


LEVEL 

AFGWC 

ECMWF 

NMC 

(mb) 

BIAS 

RMS 

BIAS 

RMS 

BIAS 

RMS 

850 

1 

21 

10 

22 

-2 

20 

700 

2 

26 

15 

29 

-4 

25 

500 

-7 

26 

22 

32 

-5 

24 

400 

-12 

30 

26 

34 

-5 

23 

300 

-26 

37 

31 

37 

-5 
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The  results  are  based  on  observations 
for  the  period  2/3/79  -  2/22/79  00Z 

Finally,  Figure  10a  shows  the  zonal  average  of  the  experimental  humidity  analysis  of  MW, 
wherein  the  operational  NMC  RH  analysis  of  00Z  14  January  1982  was  completely  replaced  by 
an  RH  analysis  derived  from  the  00Z  14  January  1982  3DNEPH  analysis  following  the 
procedure  of  Figure  7.  For  further  details  on  the  derivation  of  this  experimental  analysis,  see 
MW6  where  It  was  first  presented.  Figure  10b  gives  the  plot  of  diflerences  between  Figures  10a 
and  2a.  It  is  clear  that  Figure  5a  for  the  AFGWC  RH  analysis  resembles  Figure  10a  far  more 
closely  than  it  resembles  Figures  2-4  —  this  despite  the  fact  that  RAOBs  of  RH  were  used  in  the 
AFGWC  analysis  of  Figure  5  to  correct  a  largely  RTNEPH-based  5 LAYER  first-guess  humidity 
forecast  (within  the  octagon). 
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Figure  10.  (a)  Zonal  Average  of  RH  from  the  Experimental 
Analysis  of  3DNEPH-inferred  RH  Derived  by  Mitchell  and 
Warburton6  for  00Z.  14  January  82.  (b)  Difference  Between 
Zonal  Average  of  RH  in  Figures  10a  and  2a. 
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The  analysis  of  this  subsection  has  led  to  several  conclusions.  First,  on  average,  the 
AFGWC  Phase  I  global  humidity  analyses  are  substantially  too  dry  in  the  upper-half  of  the 
extratroplcal  troposphere.  Second,  the  cause  of  this  upper-level  dry  bias  is  the  overly  dry 
upper-level  3DNEPH  inferred  humidity  estimates  used  to  initialize  the  5LAYER  three-hour 
first-guess  humidity  forecast.  Lastly,  the  upper-level  3DNEPH  inferred  humidities  are  too  dry 
because  1)  the  zero-cloud  critical  humidity  values  in  the  upper-level  CPS  curves  are  too  low 
and  2)  the  mean  upper-level  3DNEPH  cloud  amounts  are  too  low.  We  contend  that  factor  1) 
dominates  factor  2).  To  support  this  contention,  Table  4  shows  that  Norquist1228  found  only  a 
small  upper-level  dry  bias  in  humidity  values  estimated  from  the  RTNEPH  using  an 
alternative  NMC  cloud-to-humidity  inference  scheme  —  a  scheme  originally  developed  for  use 
with  surface  cloud  reports.  This  latter  scheme  is  considered  further  in  Sections  3.1.2  and 
3.2.2. 

It  is  important  to  emphasize  that  the  height,  wind,  and  temperature  fields  of  the  AFGWC 
analyses  in  Table  1  were  also  examined  and  compared  with  those  from  other  sources.  From 
this  comparison  they  were  deemed  to  be  of  good  quality.  Though  no  details  will  be  presented 
here  for  sake  of  brevity,  the  50  kPa  AFGWC  height  analyses  in  Table  1  compared  surprisingly 
well  with  those  published  by  the  European  Centre  for  Medium  Range  Weather  Forecasts 
(ECMWF),32  especially  considering  the  relatively  poor  resolution  of  the  AWAPS  Phase  I  AFGWC 
analyses  and  GSM  forecast  first-guess  fields. 


2.2  Global  Humidity  Forecasts 

This  section  will  examine  the  characteristics  of  the  4-day  global  humidity  forecasts  of  the 
GL  baseline  GSM.  Recall  from  Section  1  that  the  term  "baseline"  here  refers  to  a  particular 
configuration  of  the  GL  GSM  that  parallels  the  currently  operational  AFGWC  GSM.  Humidity 
fields  from  four  GSM  forecasts  out  to  96  hours  are  examined  here.  One  GSM  forecast  is 
presented  for  each  of  the  four  types  of  GSM  initial  analysis  sources  listed  in  Table  1 .  Thus  all 
four  GSM  forecasts  are  NH  winter  season  cases.  The  basetime  and  source  of  the  four  initial 
analyses  are  shown  in  Figure  1 1 . 


32  Staff  (1986)  Daily  Global  Analyses:  Operational  Data  Assimilation  System,  January  - 
March  1985,  European  Center  for  Medium  Range  Weather  Forecasts,  Reading,  Berkshire,  U.K. 
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Figure  11.  NH  Mean  RH  at  Six  Mandatory  Pressures  as  a  Function  of  Forecast  length  in  Four 
GSM  Forecasts  Initialized  from  Four  Separate  Initial  Analyses  as  Labeled. 


Like  the  AFGWC  GSM,  the  GL  baseline  GSM  explicitly  forecasts  moisture  on  only  the  lowest 
7  of  the  12  model  sigma  layers,  that  Is,  on  those  sigma  layers  that  are  essentially  contained  In 
the  troposphere  as  shown  In  Table  5.  The  GSM  postprocessor  eventually  outputs  the  GSM  RH 
forecast  on  the  lowest  six  mandatory  pressure  surfaces.  Figure  1 1  shows  the  96-hour  temporal 
trend  of  the  NH  mean  RH  at  each  of  these  six  pressure  surfaces.  The  examination  In  this 
section  shall  emphasize  NH  behavior,  as  the  5 LAYER  cloud  forecasts  used  to  compare  with 
GSM  cloud  forecasts  In  Section  3.2  are  exclusively  NH  forecasts. 

Note  In  Figure  1 1  that  the  lower  (upper)  two  plots  are  for  the  two  forecasts  initialized  from 
humidity  analyses  that  directly  or  Indirectly  utilize  (Ignore)  3DNEPH/RTNEPH  cloud  data. 
Figure  12,  generated  from  the  data  In  Figure  11,  shows  the  O-hour  and  96-hour  vertical  profiles 
of  NH  mean  RH  for  the  four  GSM  forecasts  in  Figure  11,  Also  as  a  reference,  Figure  12  shows 
the  vertical  profile  of  the  NH  mean  RH  obtained  from  the  FGGE  n-B  RAOB-based  data  in 
Figure  2a. 


Table  5.  The  o-structure  of  the  Baseline  GL  GSM. 
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Figure  12.  Vertical  Profiles  of  the  NH  Mean  RH  at  the  (a)  O-hr  Initial  Time  and  (b)  96-hr 
Forecast  Time  for  the  Four  GSM  Forecast  Cases  in  Figure  11.  Also  given  as  a  reference 
is  the  vertical  profile  of  NH  mean  RH  as  derived  from  the  FGGE  n-B  zonal  average 
RH  values  in  Figure  2. 


Several  Important  features  are  apparent  in  Figures  11  and  12.  First,  reiterating  several 
O-hour  characteristics  noted  earlier  in  Section  2.1.2,  the  two  3DNEPH/RTNEPH-influenced 
O-hour  humidity  analyses  are  much  drier  in  the  mean  at  30  and  40  kPa  than  the  other  two 
analyses.  Of  these  two  dry  upper-level  analyses,  the  NMC/3DNEPH  analysis  is  the  driest, 
presumably  because  unlike  the  AFGWC  analysis,  it  has  not  been  "corrected"  by  actual  humidity 
RAOBs.  At  50  kPa,  the  O-hour  mean  RH  values  of  all  but  the  NMC/3DNEPH  analysis  agree 
fairly  closely  with  each  other.  At  the  70  and  85  kPa  levels,  the  two  3DNEPH/RTNEPH- 
influenced  analyses  are  somewhat  more  moist  (compare  again  Figures  4b  and  5b).  Lastly  at 
the  100  kPa  level,  the  mean  RH  values  of  all  four  analyses  are  in  near  agreement. 

Of  greater  interest  in  this  section  are  the  distinct  GSM  forecast  trends  apparent  in 
Figures  11  and  12.  Most  notably,  the  GSM  quickly  evolves  or  adjusts  at  all  levels  toward  what 
we  shall  refer  to  here  as  the  model's  preferred  mean  humidity  state.  This  asymptotic,  mean 
GSM  humidity  state  is  somewhat  dry  (moist)  at  levels  below  (above)  50  kPa,  compared  to  the 
RAOB  mean  state.  More  importantly  here,  the  96-hour  vertical  profile  of  GSM  hemispheric 
mean  RH  values  is  largely  independent  of  the  O-hour  initial  humidity  profile.  This  type  of 
behavior  in  forecast  mean  RH  fields  is  a  common  characteristic  in  large-scale  forecast  models 
as  a  class.33  There  is  clearly  strong  forcing  of  the  GSM  forecast  humidity  distribution  due  to 
modeled  internal  vertical  motions  and  dlabatlc  physics.  The  majority  of  the  adjustment  of 
GSM  mean  humidity  occurs  in  the  first  24  (48)  hours  in  the  upper  (lower)  levels.  Not 
surprisingly,  the  magnitude  of  this  rapid  initial  humidity  adjustment  is  largest  in  the  upper 
levels  of  the  forecasts  from  the  AFGWC  and  NMC/3DNEPH  analyses,  as  these  levels  depart 
more  from  both  the  observations  and  the  GSM's  preferred  mean  humidity  state. 

It  is  instructive  to  examine  the  GSM  temporal  humidity  trends  in  greater  detail  across  the 
model's  horizontal  spatial  domain.  Corresponding  to  Figures  lib  and  lid.  Figures  13  and  14 
show  the  zonally-averaged  48-hour  RH  forecast  and  forecast  trends  (forecast  minus  analysis) 
for  the  NMC/OP  and  AFGWC  cases,  respectively.  Note  that  Figure  13b  is  the  difference  of 
Figure  13a  and  Figure  4a;  and  similarly.  Figure  14b  is  the  difference  of  Figure  14a  and  the 
corresponding  single  initial  analysis  field  taken  from  the  eleven  analysis  fields  averaged  in 
Figure  5a. 


33  Lejenas,  H.  (1979)  Initialization  of  moisture  in  primitive  equation  models,  Mon.  Wea.  Rev. 
107:1299-1305. 
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Figure  13.  (a)  Zonal  Average  of  RH  for  the  48-hr  GSM  Forecast 
Initialized  from  the  NMC  Analysis  of  OOZ,  14  January  1982. 

(b)  Difference  Between  the  Above  Zonal  Average  and  that  of  the 
Corresponding  Initial  Analysis  {Figure  4a). 
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In  the  NMC/OP  case  of  Figure  13,  the  modest  increase  in  mean  RH  above  50  kPa  is 
pervasive  across  virtually  all  NH  latitudes.  The  general  NH  drying  at  70-85  kPa  in  Figure  13  is 
interrupted  at  a  few  latitudes  by  moistening.  Though  not  shown,  the  zonal  mean  RH  trends  in 
the  GSM  forecast  from  the  17  February  1979  NMC  FGGE  III-A  analysis  are  very  similar  to 
those  in  Figure  13.  The  RH  trends  in  Figures  1  la-1  lb  and  Figure  13  are  very  likely 
representative  of  the  mean  RH  trends  AFGWC  will  observe  in  its  GSM  forecasts  initialized 
from  the  AWAPS  Phase  II  moisture  analyses  after  March  1988,  because  then  the  AFGWC  RH 
observation  sources  (RAOBs)  and  global  data  assimilation  forecast  model  (GSM)  are  very 
similar  to  those  represented  in  Figure  13. 

As  expected  from  Figure  11,  the  Figure  14  RH  trends  in  the  GSM  forecast  from  the  AFGWC 
analysis  are  much  larger  than  those  of  the  NMC/OP  case  in  Figure  13.  The  RH  trends  evident 
in  Figure  14  are  very  significant  to  the  experimental  GSM  cloud  forecasts  presented  in 
Section  3,  as  all  of  these  GSM  cloud  forecasts  are  executed  from  January  1985  AFGWC  analyses 
in  order  to  allow  comparisons  with  the  January  1985  5LAYER  cloud  forecasts  archived  for  this 
study.  In  Figure  14  the  large  moistening  at  upper-level  extratropical  locations  (that  is,  within 
the  octagon  domain)  increases  with  height  and  dominates  all  trends.  Also  of  note,  the  RH 
trends  at  the  70  -  100  kPa  levels  at  latitudes  within  the  NH  octagon  show  a  modest  drying  on 
average. 

Comparing  Figure  14  with  Figure  5,  it  is  important  to  recognize  that  in  the  zonal  mean 
sense,  the  GSM  RH  forecast  is  actually  a  significant  improvement  over  the  initial  AFGWC  RH 
analysis.  Figures  14b  and  5b  show  that  the  GSM  RH  forecast  is  growing  more  moist  (more  dry) 
primarily  in  those  areas  where  the  AFGWC  RH  analysis  is  erroneously  too  dry  (too  moist)  with 
respect  to  the  FGGE  II-B  observed  zonal  mean  moisture  distribution.  That  is,  the  forecast  RH 
mean  state  in  Figure  14a  agrees  more  closely  with  the  observed  RH  mean  state  in  Figure  2a 
than  does  the  AFGWC  analyzed  mean  state  of  Figure  5a.  This  characteristic  whereby  the  0-48 
hour  GSM  RH  zonal  mean  forecast  trends  are  large  (especially  at  upper  levels  within  the 
octagon)  because  the  forecast  needs  to  correct  errors  in  the  mean  state  of  the  AFGWC  RH 
analysis  is  a  serious  detriment  to  the  ability  of  the  GSM  to  yield  short-range  competitive  cloud 
forecasts  in  Section  3.2. 

Next,  Figures  15-18  expand  the  present  focus  on  the  temporal  trends  of  the  GSM  humidity 
forecast  to  the  model's  entire  horizontal  domain.  For  the  NMC/OP  case.  Figures  15  and  16 
provide  NH  shaded  displays  of  the  0-hour  NMC  RH  analysis  and  the  associated  48-hour  GSM 
RH  forecast,  respectively,  for  all  six  moist  mandatory  levels.  Figures  17  and  18  provide  the 
corresponding  displays  for  the  AFGWC  case.  A  comparison  of  Figures  15e-f  and  17e-f  vividly 
depicts  the  pervasiveness  of  the  very  dry  character  of  the  AFGWC  30-40  kPa  analyses  within 
the  AFGWC  octagon  domain.  The  48-hour  RH  forecasts  at  30  kPa  in  Figures  16e-f  and  18e-f 
show  how  quickly  the  GSM  48-hour  RH  forecast  from  the  AFGWC  RH  analysis  has  recovered 
from  the  Initial  dry  bias  and  produced  an  RH  forecast  field  whose  overall  synoptic  character  is 
already  remarkably  similar  to  the  48-hour  RH  forecast  from  the  NMC/OP  case. 
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Figure  15.  NH  NMC  RH  Analysis  for  00Z,  14  January  1982  at 
Levels  (a)  100,  (b)  85,  (c)  70,  (d)  50.  (e)  40  and  (1)  30  kPa.  (The 
shading  key  In  Figure  6  applies  to  all  shaded  displays  of  RH 

fields). 
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Figure  16.  NH  GSM  48-hr  RH  Forecast  at  levels  (a)  100,  (b)  85, 

(c)  70,  (d)  50,  (e)  40  and  (0  30  kPa  tor  the  Forecast  Initialized  From 
the  NMC  Analysis  of  Figure  15 
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Figure  17.  NH  AFGWC  RH  analysis  '  -»  OOZ,  17  January  1985  at 
Levels  (a)  100,  (b)  85,  (c)  70,  id)  00.  !< ;  40  and  (0  30  kPa.  (The  shading 
key  In  Figure  5  applies  to  all  shaded  displays  of  RH  fields). 


Figure  18.  NH  GSM  48-hr  RH  Forecast  at  Levels  (a)  100,  (b)  85, 

(c)  70.  (d)  50,  (e)  40  and  (0  30  kPa  for  the  Forecast  Initialized  from 
the  GSM  Forecast  Initialized  from  the  AFGWC  Analysis  of  Figure  17. 
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But  the  primary  Interest  here  in  Figures  15-18  is  the  100,  85,  and  70  kPa  levels,  where 
some  new  characteristics  not  evident  In  Figures  11-14  are  revealed.  By  comparing  Figures 
15b-c  to  Figures  16b-c  and  Figures  17b-c  to  Figures  18b-c,  one  notes  that  both  48-hour  forecast 
fields  show  higher  "contrast"  within  the  display  field  than  their  0-hour  analysis  counterparts. 
That  is,  the  RH  forecast  displays  show  more  areas  at  the  dry  and  moist  ends  of  the  humidity 
range  and  sharper  gradients  between  these  moist  and  dry  areas  than  the  RH  analysis  displays. 
Most  significant  is  the  vivid  synoptic  patterns  of  narrow  mo<st  fronted  bands  and  post-frontal 
"dry  tongues"  that  have  accompanied  the  emergence  of  more  extreme  wet  and  dry  areas  in  the 
GSM  forecast. 

Quantitatively,  this  Increase  In  contrast  of  the  RH  display  in  the  GSM  forecast  shows  up  in 
terms  of  a  changing  frequency  distribution  of  RH  during  the  forecast.  As  an  example,  for  the 
NH  octagon  domain,  Figure  19  shows  the  temporal  change  toward  an  increase  In  the  percent 
frequency  of  occurrence  of  GSM  RH  forecast  values  in  both  the  moist  and  dry  Intervals  of  the 
RH  range  for  the  85,  70,  and  50  kPa  levels  In  both  the  NMC/OP  and  AFGWC  forecast  cases.  The 
above  occurs,  of  course,  at  the  expense  of  fewer  RH  values  occurring  In  the  Intermediate 
humidity  categories  (not  shewn).  In  Figure  19,  the  tendency  of  the  GSM  to  drive  the  moisture 
distribution  toward  the  dry  and  moist  ends  simultaneously  is  particularly  highlighted  by  the 
85  and  70  kPa  levels,  where  there  is,  in  fact,  a  modest  decrease  in  the  octagon  mean  RH  value 
(not  shown)  during  the  forecast  period. 
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The  above  spin-up  of  the  GSM  toward  a  more  widespread  occurrence  of  extremes  in  the 
moisture  field  is  essentially  accomplished  by  48  hours  in  Figure  19.  This  result  is  a  classic 
manifestation  of  the  well-known  "spin-up':  period  associated  with  large-scale  models.  The  GL 
and  AFGWC  GSMs,  like  virtually  all  present-day  large-scale  models,  begin  with  a  quasi- 
balanced,  quasi-nondivergent  initial  wind  analysis,  owing  to  the  application  of  a  nonlinear 
normal  mode  initialization  to  eliminate  spurious  large- amplitude  gravity  waves.  Because  of 
the  largely  nondivergent  initial  wind  fields,  the  short  range  (0-6  hours)  GSM  vertical  motion 
fields  and  frontal  convergence  zones  are  relatively  weak,  but  subsequently  intensify  as  the 
dynamics  of  the  forecast  model  act  to  increase  the  magnitude  of  the  horizontal  divergence  (see 
Brenner  et  al.7).  As  the  GSM  forecast  proceeds,  the  attendant  increase  in  the  magnitudes  of  the 
vertical  velocities  and  their  organization  into  synoptic  patterns  (at  least  at  extratropical 
latitudes)  results  in  the  aforementioned  development  of  distinct  synoptic  frontal  structures  in 
the  forecast  RH  field.  Indeed,  it  is  the  vividness  of  the  synoptic  frontal  patterns  in 
Figures  16b-c  and  18b-c  that  leads  us  in  the  first  place  to  attempt  cloud  forecasting  with  a 
large-scale  global  forecast  model. 

Further  examination  of  the  moist  frontal  bands  in  Figures  16  and  18  reveals  an 
unexpectedly  deep  vertical  extent  to  these  bands.  The  occurrence  of  a  saturated  band  at  one 
layer  appears  highly  correlated  with  the  simultaneous  presence  of  a  similar  band  at  virtually 
all  other  tropospheric  layers  above  and  below.  The  moist  band  over  the  west  coast  of  Canada 
in  Figure  16  represents  one  of  many  extratropical  examples  in  Figures  16  and  18.  Outside  the 
tropics  in  winter,  we  can  rule  out  deep  convection  as  a  dominant  process.  Hence  we  conclude 
that  the  deep  character  of  the  extratropical  saturated  fronts  in  the  GSM  is  a  systematic  error. 
Such  deeply  saturated  fronts  are  substantially  less  frequent  in  the  initial  RH  analyses  of 
Figures  15  and  17,  even  in  data  rich  areas.  A  similar  systematic  trend  toward  overdeveloping 
the  vertical  extent  of  moist  frontal  bands  has  been  noted  also  in  the  ECMWF  global  model.3 
The  trend  observed  in  the  present  study  likely  contributes  to  the  lack  of  a  vertical  gradient  in 
mean  RH  between  the  levels  of  70  and  30  kPa  in  the  GSM  96-hour  NH  mean  RH  profiles  in 
Figure  12b. 

To  conclude  this  section,  we  focus  on  a  systematic  RH  forecast  trend  apparent  at  100  kPa 
over  tropical  and  subtropical  oceans,  especially  away  from  coastlines  where  air  of  continental 
origin  may  suppress  the  trend.  Comparing  Figure  15a  to  Figure  16a  and  Figure  17a  to  Figure 
18a,  we  observe  In  these  warm  ocean  areas  an  increase  in  RH  to  virtually  saturated  values  in 
Just  48  hours.  Some  of  these  very  moist  low-latitude  ocean  areas  extend  northward  of  the 
octagon  boundary  and  will  affect  the  cloud  forecasts  verified  over  the  octagon  in  Section  3.2. 

The  cause  of  this  systematic  trend  is  the  absence  of  a  representation  of  vertical  diffusion 
above  the  lowermost  model  level  in  the  baseline  GSM.  In  the  bottom  level  of  the  GSM,  surface 
evaporation  over  oceans  occurs  at  the  potential  rate,  as  governed  by  a  standard  bulk 
aerodynamic  drag  formula  and  the  specified  sea-surface  temperatures.®  Over  warmer  low- 
latitude  waters  this  evaporation  rate  can  be  substantial,  if  the  air  is  not  already  near  satura¬ 
tion.  As  the  evaporation  proceeds,  the  lowest  model  layer  continues  to  moisten  rapidly  (until 
the  approach  to  saturation  shuts  the  evaporation  down),  because  no  vertical  diffusion  is 
present  in  the  GSM  above  the  first  layer.  Such  a  vertical  diffusion,  if  it  were  present,  would 
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act  as  a  continual  sink  of  surface  layer  water  vapor  owing  to  the  vertical  turbulent  flux  to 
overlying  layers. 

Over  land  the  GSM  layer-one  mean  RH  trend  Is  reversed,  as  there  the  GSM  ignores  surface 
evaporation,  yet  precipitation  still  acts  as  a  moisture  sink.  These  reverse  trends  at  100  kPa 
over  land  and  water  undoubtedly  cancel  each  other  to  a  significant  degree  in  the  NH  mean 
trends  shown  at  100  kPa  in  Figures  11  and  12.  However,  in  Figures  13b  and  14b,  between 
latitudes  40S  and  20N,  where  the  percentage  of  land  area  is  small,  the  mean  moistening  over 
oceans  at  100  kPa  during  the  forecast  is  evident,  particularly  in  the  January  1985  case  of 
Figure  14.  This  latter  case  is  one  of  the  two  cloud  forecast  cases  presented  later  in  Section  3.2. 


3.  EXPERIMENTAL  GLOBAL  CLOUD  FORECASTS 


3.1  Diagnostic  Cloud  Forecast  Schemes 


There  are  a  myriad  of  challenging  difficulties  and  issues  in  developing  a  cloud  forecast 
scheme  for  a  large-scale  NWP  model.  The  NWP  model  cloud  forecast  studies  by 
J.M.  Slingo3-34  35  provide  excellent  discussions  and  examples  of  these  problems  and  issues  and 
the  reader  is  encouraged  to  review  these  references.  For  the  reasons  emphasized  at  the 
beginning  of  Section  2  and  following  the  recommendation  of  Slingo,  this  study  will  examine 
only  diagnostic  cloud  forecast  schemes;  that  is,  schemes  that  diagnose  cloud  properties  (in  this 
case  layer  and  total  fractional  coverage)  empirically  or  statistically  from  other  model  output 
fields. 

Diagnostic  cloud  forecast  schemes  themselves  encompass  two  broad  types,  which  are  1) 
schemes  that  utilize  only  the  model's  forecast  humidity  field  and  2)  schemes  that  utilize  both 
the  humidity  field  and  additional  model  output  fields  such  as  vertical  velocity,  convective  and 
non-convective  precipitation  rates,  vertical  temperature  and  moisture  gradients  in  the  PBL, 
and  surface  fluxes  of  moisture  and  heat,  among  others.  A  diagnostic  scheme  of  type  1  was  used 
in  the  ECMWF  global  forecast  model  until  May  1985;  and  the  performance  of  this  cloud 
scheme,  known  as  the  Geleyn  scheme1,  is  well  documented  by  Slingo.3  Since  May  1985,  a 
diagnostic  scheme  of  type  2  has  been  used  in  the  ECMWF  model,  and  the  performance  of  this 
cloud  scheme,  known  as  the  Slingo  scheme,  is  also  well  documented  by  Slingo.35 

The  choice  of  which  type  of  diagnostic  cloud  scheme  to  employ  in  a  given  large-scale  model 
is  dictated  for  the  most  part  by  the  complexity  and  veracity  of  the  model's  parameterized 
physical  processes.  Type-2  diagnostic  schemes  require  fairly  reliable  deep  and  shallow 
convection  schemes,  fairly  complex  PBL  schemes  that  employ  at  least  three  or  four  model 


34  Slingo,  J.  (1980)  A  cloud  parametrization  scheme  derived  from  GATE  data  for  use  with  a 
numerical  model,  Q.J.R.  Met.  Soe  106:747-770. 

35  Slingo,  J.  (1987)  The  development  and  verification  of  a  cloud  prediction  scheme  for  the 
ECMWF  model,  Q.J.R  Met.  Soc.,  113:899-927. 
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layers  In  the  lowest  kilometer,  and  radiative  forcing  schemes  with  a  diurnal  cycle.  As  detailed 
in  the  Introduction,  AFGWC's  currently  operational  GSM  (and  hence  GL  s  baseline  GSM  used 
here)  lacks  such  physical  complexity.  Although  the  AFGWC  GSM  does  include  a  deep 
convection  scheme  patterned  after  the  widely  used  parameterization  of  Kuo,10  It  also  Includes 
several  additional  constraints  on  the  occurrence  of  deep  convection  not  found  in  Kuo's  original 
scheme.  Examination  by  GL18  of  the  resulting  convective  precipitation  rates  and  frequency  of 
convective  points  showed  vastly  smaller,  almost  negligible,  time  and  space  averaged  values  for 
these  quantities  In  the  AFGWC  GSM  (or  baseline  GL  GSM)  compared  to  other  models  and 
estimated  climatologies. 

As  a  consequence  of  these  baseline  GL  GSM  characteristics,  all  experimental  GSM  cloud 
forecast  schemes  tested  in  this  study  and  presented  in  this  section  were  purposely  limited  to 
type-1  diagnostic  schemes.  Because  of  the  limitations  in  GSM  parameterized  physics,  we 
cannot  expect  the  GSM  cloud  forecasts  presented  here  to  perform  well  depicting  most  tropical 
cloud  regimes  (which  are  dominated  by  deep  and  shallow  convection  and  have  a  strong  diurnal 
character  over  land)  or  those  low  cloud  regimes  that  are  driven  by  local  PBL  processes  and 
their  Interaction  with  terrain  and  radiation  (such  as  fair  weather  "popcorn"  cumulus  or 
persistent  marine  stratocumulus).  But  the  AFGWC  5LAYER  model,  which  also  uses  a  type-1 
diagnostic  cloud  scheme,  similarly  lacks  an  ability  to  forecast  such  tropical,  convective,  or 
PBL  cloud  regimes,  as  it  too  lacks  physical  treatments  of  convection,  radiation,  and  the  PBL.4 

In  terms  of  cloud  forecast  skill  then,  the  strong  suit  of  both  the  baseline  GSM  and  5LAYER 
models  is  expected  to  be  prediction  of  the  extratropical  large-scale  cloud  organizations 
associated  with  synoptic  circulations  —  an  association  often  observed  in  satellite  imagery. 
Certainly  In  the  case  of  the  baseline  GSM,  synoptic  patterns  were  distinctly  present  In  the  RH 
forecast  displays  of  Figures  16  and  18.  Because  of  some  similarity  in  the  strengths  and 
weaknesses  of  the  baseline  GSM  and  5 LAYER  models  (and  recalling  that  5 LAYER  uses  as  input 
the  GSM's  horizontal  and  vertical  wind  velocity  forecast),  it  is  reasonable  to  expect  baseline 
GSM  cloud  forecasts  to  be  competitive  with  5LAYER  cloud  forecasts. 

3.1.1  GLOBAL  CLOUD  FORECAST  PROCEDURES 

The  use  of  a  type-1  diagnostic  cloud  scheme,  which  again  utilizes  only  the  model's  moisture 
forecast  (plus  possibly  the  temperature  forecast  to  allow  conversion  to  other  moisture 
variables),  offers  one  important  expedient  advantage  over  more  complex  schemes.  Namely, 
since  a  model's  moisture  and  temperature  forecast  are  standard  model  output  fields,  a  type-1 
diagnostic  cloud  scheme  can  be  run  as  a  "cloud  postprocessor",  executed  external  to  the  GSM 
model  (even  on  a  separate  computer  system).  The  type-1  diagnostic  cloud  schemes  tested  here 
Indeed  were  executed  as  a  separate  cloud  postprocessor  after  executing  the  baseline  GL  GSM. 

Specifically,  after  GSM  execution,  the  first  task  was  to  execute  the  standard  non-cloud  GSM 
postprocessor  wherein  (1)  GSM  forecast  spectral  coefficients  of  temperature  (T)  and  specific 
humidity  (Q)  were  transformed  to  a  2.5°  latitude/longitude  grid  (following  the  AFGWC  grid 
convention)  on  the  GSM's  seven  moist  sigma  surfaces,  (2)  T  and  Q  at  each  grid  point  were 
converted  to  RH  (which  due  to  well  known  moisture  spectral  truncation  errors  was  then 
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checked  and  constrained  to  be  between  1  and  100  percent),  and  (3)  T  and  RH  were  vertically 
Interpolated  from  sigma  surfaces  to  the  lowest  six  mandatory  pressure  surfaces.  After  the 
standard  GSM  postprocessor  executed  these  three  steps,  the  cloud  postprocessor  was  executed, 
which  Involved  four  steps.  First,  the  T  and  RH  fields  on  the  six  pressure  surfaces  from  step  3 
above  were  input  and  biquadrallcally  Interpolated  to  AFGWC's  fully  hemispheric  1 /2-mesh 
129  X  129  grid  (It  is  the  RH  fields  from  this  slep  that  are  displayed  in  Figures  15- 18).  Next,  RH 
was  converted,  If  necessary,  to  another  moisture  variable  (using  T  when  needed).  Then,  the 
chosen  type-1  diagnostic  cloud  st  !  •  me  ".as  applied  to  convert  the  moisture  variable  to  a  cloud 
amount  at  each  moist  pressure  surliue.  Finally,  a  vertical  stacking  algorithm  was  applied  to 
the  layer  cloud  amounts  to  obtain  a  i-nal  cloud  amount.  In  the  secund  step,  we  imfAlcltly  view 
the  GSM  forecast  moisture  value  i  given  pressure  level  as  a  mean  layer  value  for  the 
corresponding  layer  centered  around  ihe  given  pressure  level  and  hence  view  the  resulting 
cloud  amount  as  a  layer  cloud  a  n  m  a  The  same  vertical  stacking  algorithm,  details  of  which 
are  given  in  AppendLx  A,  was  u-  >  3  for  all  the  cloud  schemes  (with  minor  variations  described 
in  Appendix  A).  In  the  stacking  algorithm,  whenever  a  given  pressure  level  Is  found  to  be 
below  the  terrain  height  defined  on  AFGWC's  129  X  129  grid,  the  layer  cloud  amount  is  set  to 
zero. 

The  AFGWC  93  X  101  5LAYER  octagon  1/2-mesh  grid  is  an  exact  pointwise  subset  of  the 
129  X  129  hemispheric  1 /2-mesh  grid  cited  above.  All  the  total  cloud  forecast  verification 
statistics  presented  later  were  computed  on  the  5LAYER  1 /2-mesh  grid,  as  the  verifying 
RTNEPH  cloud  analyses  archived  for  tills  study  were  also  represented  on  this  grid. 

In  any  future  follow  on  study,  one  is  strongly  urged  to  reduce  postprocessing  error  by 
reducing  and  modifying  the  above  steps  as  follows:  carry  out  the  o-surface  horizontal  spectral 
transform  to  a  liner  scale  (that  is,  1.0  or  0.5°)  rather  than  2.5°  latitude/longitude  grid  before 
interpolating  to  the  verifying  1/2-mesh  grid  and  remain  on  the  GSM  0-surfaces,  thereby 
eliminating  the  errors  of  vertical  interpolation  to  pressure  surfaces  (also  the  grouping  of 
0-layers  into  cloud  layers  can  occur  then  on  physical  grounds  rather  than  according  to  the 
arbitrary  location  of  pressure  surfaces).  Presently,  AFGWC  does  not  store  the  original  GSM 
forecast  spectral  coefficients  on  0-surfaces  beyond  12  hours  (the  storage  from  0-12  hours  is  for 
the  global  data  assimilation  cycle).  In  the  present  study,  we  chose  to  utilize  only  GL  GSM 
output  fields  corresponding  to  those  routinely  stored  by  AFGWC  GSM  postprocessing. 

To  provide  the  5LAYER  cloud  forecasts  for  GSM  comparison  purposes  in  this  study,  AFGWC 
originally  archived  and  provided  GL  with  12  NH  5LAYER  48-hour  forecast  cases  (three  each 
during  the  AWAPS  Phase  I  months  of  April,  July,  and  October  1984,  and  January  1985).  As 
companion  data  sets,  AFGWC  also  archived  the  corresponding  12  AWAPS  global  objective 
analyses  b  L^.  used  as  the  initial  conditions  (except  for  humidity)  for  the  subsequent  off-line 
GL  GSM  cloud  forecast  experiments  for  tills  study.  Before  examining  these  GL  GSM  cloud 
forecasts  here,  it  is  important  to  emphasize  why  this  study  did  not  simply  use  the  archived 
data  (in  a  type-1  diagnostic  cloud  postprocessor)  from  the  operational  AFGWC  AWAPS  GSM 
temperature  and  humidity  forecasts  for  the  above  12  cases.  The  reasons  were  fourfold: 

1)  The  quality  of  the  AWAPS  Phase  I  humidity  analyses  used  as  initial  moisture 
conditions  for  the  AWAPS  GSM  forecasts  was  poor  for  the  several  reasons  detailed  in 
Section  2.1.1. 
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2)  The  AWAPS  Phase  i  br'dial  humidity  analyses,  unlike  die  blAf  ER  Initial  humidity 
analyses,  neglected  to  utilize  RTNEPH  cloud  data  as  a  surrogate  humidity  data  source.  A 
priori,  this  was  viewed  as  giving  5LAYER  mi  unfair  advantage  over  die  GSM,  especially  in  the 
short  range,  since  the  RTNEPH  was  to  be  used  as  the  cloud  verification  database. 

3)  The  effective  horizontal  spatial  resolution  of  the  AWAPS  Phase  1  GSM  (rhomboldal  20) 
was  considerably  coarser  than  that  of  SLAYER  (1/2-rnesh  or  100  nm  at  60°  N)  or  that  of 
AFGWC's  present  Phase  II  GSM  (rhomboldal  40). 

4)  The  study  of  Mitchell  and  Yang1-1  showed  Lhai  the  methods  of  vertically  inteipolating 
moisture  fields  between  pressure  and  sigma  surfaces  in  the  preprocessor  and  postprocessor  ol 
tiie  AWAPS  Phase  I  GSM  were  sources  of  significant  GSM  moisture  forecast  error,  especially  in 
the  short  range. 

For  the  above  reasons,  this  study  utilized  separate  GL  GSM  forecasts  executed  off-line  on 
the  aKSCC-K  CRAY- 1/5  and  designed  to  mitigate  these  four  factors.  !  he  GL  GSM  horizontal 
resolution  was  Increased  to  rhombokiai-40  and  the  preprocessor  and  postprocessor  were 
changed  to  utilize  the  vertical  moisture  interpolation  methods  recommended  by  Ml'°ii*»)i  and 
Yang.14  Most  importantly,  as  described  in  Section  2.1.1,  the  AWAPS  Phase  I  humidity 
analyses  were  rep1, iced  here  with  AFGWC's  alternative  global  humidity  analyses  that  utilized 
RflNEPH  clouo  data  via  a  SLAYER  ttrst -guess  and  therefore  more  closely  resembled  the  initial 
humidity  analyses  of  5LAYER. 

Despite  the  good  intentions  of  using  AFGWC's  non-AWAPS  global  moisture  analyses,  the 
results  of  Section  2.1.2  show  that  the  latter  analyses  suffer  from  their  own  unique  problems, 
including  a  severe  extra  tropical  upper-level  dry  bias.  More  pertinent  here  Is  the  additional 
problem  described  in  Section  2.1.1  whereby  the  nine  April  through  October  1984  cases  proved 
to  be  unusable  in  the  parallel  GL  GSM  cloud  forecast  experiments  owing  to  the  zero  humidifies 
stored  at  all  tropical  points.  This  in  fact  left  only  the  three  January  19«5  cases  as  candidates 
for  further  GSM  and  5 LAYER  comparisons. 

In  the  end,  only  the  two  51  AYER  cases  for  00 Z  on  17  January  and  24  January  1985  were 
accessible  for  the  comparisons  here  (physical  errors  on  (he  ar<  hive  I  ape  prevented  use  of  the 
third  January  1985  ,  a  sc).  Recall  that  as  background  lor  die  O'.  GSM  <  loud  telecasts  to  be 
examined  for  these  iwo  cases  below,  the  GL  GSM  humidity  forecasts  bn  ,lu  !  7  January  case 
were  examined  in  detail  in  See! Ion  2.2. 

3.1.2  FREVIOLS  DIAGNOSTIC  CLOHD  SCHEMES 

In  the  sCps  ol  the  cloud  postprocessor  described  in  the  previous  •'-■(■lion,  the  <  entral  step 
invokes  the  chosen  type-  i  diagnostic  cloud  scheme.  This  study  vJ.i  lest  and  examine  four  such 
schemes.  The  first  three  schemes,  which  are  considered  in  ibis  sew  ion.  air  previously 
available  schemes  «  urn-iitlv  or  formerly  user)  by  other  .'AVI'  centers  -tieeitu  ally  one  scheme 
each  iron)  AEGWC.  NMf’  anti  El  MWK  The  lourtii  s.  home,  who  h  h.s>  ! a , ,  variants,  is  a  new 
CL  cloud  .scheme  developed  spi-citli  allv  as  a  product  ot  tins  study  and  is  desc  nbed  separately  In 
Section  3. 1 .3. 


Since  the  experimental  GSM  cloud  forecasts  are  to  be  compared  with  the  cloud  forecasts  of 
the  AFGWC  5LAYER  model,  an  obvious  first  candidate  for  a  GSM  cloud  scheme  Is  the 
diagnostic  scheme  used  by  the  5LAYER  model.  The  latter  scheme  Is  the  previously  Introduced 
AFGWC  CPS-cloud  scheme  given  In  Figure  8.  In  the  5LAYER  model  code,  the  CPS-cloud  curves 
In  Figure  8  are  actually  represented  as  four  look-up  tables  (one  for  each  level  in  Figure  8) 
specified  in  Increments  of  one  percent  cloud  amount  (see  Appendix  B  of  Crum4). 

In  its  treatment  of  moisture,  precipitation  and  clouds,  the  5LAYER  model  Includes  explicit 
prognostic  equations  only  for  the  moisture  variable  CPS  (defined  earlier  in  Section  2.1.2). 

CPS  is  forecast  at  five  layers  whose  msd-  points  are  given  by  the  four  constant  pressure  levels 
of  85.  70.  50,  and  30  kPa  plus  a  so-called  "gradient"  level,  which  is  a  terrain  following  lower 
level  defined  to  be  at  6  kPa  above  the  surface.  At  each  of  these  five  layers,  the  5LAYER  model 
applies  the  respective  curve  in  Figure  8  to  empirically  convert  each  CPS  forecast  value  to  a 
cloud  amount.  For  the  gradient  laye: ,  the  85  kPa  CPS-cloud  curve  is  utilized.  Finally,  total 
cloud  is  derived  in  5LAYER  by  applying  a  vertical  stacking  algorithm  as  given  in  Appendix  A. 

In  the  GSM  application  here  of  the  AFGWC  cloud  scheme,  the  GSM  cloud  postprocessor 
derives  forecast  CPS  values  at  all  six  moist  mandatory  pressure  levels  from  the  GSM  forecast 
values  of  RH  and  T.  Then  the  AFGWC  CPS-cloud  curves  are  applied  to  obtain  the 
corresponding  cloud  amounts.  The  30  and  85  kPa  CPS-cloud  curves  also  are  applied  to  the  40 
and  100  kra  levels,  respectively.  (In  this  regard,  recall  the  very  close  similarity  in  Figure  8 
between  the  50  and  30  kPa  curves  and  the  5 LAYER  model  use  of  the  85  kPa  curve  with  the 
gradient  level.) 

Looking  ahead  to  the  quality  of  the  GSM  cloud  forecasts  obtained  with  the  AFGWC  scheme, 
we  recall  from  Figure  9  that  this  scheme  yields  very  low  values  of  the  critical  relative 
humidity,  RHc,  in  the  range  10  <  RHc  <  30  percent  for  the  30-50  kPa  regime.  As  will  be  shown, 
these  critical  values  are  far  below  those  obtained  from  the  other  cloud  schemes  to  be 
considered  here  and,  as  previously  shown  in  Section  2.1.2,  contributed  substantially  to  the 
significantly  dry  nature  of  the  AFGWC  upper-level  humidity  analysis  within  the  NH  octagon 
domain.  In  contrast,  Figure  1 1  showed  that  the  GSM  NH  mean  RH  forecast  value  in  the  30-50 
kPa  regime  is  between  40  and  50  percent  after  24  hours,  even  when  the  GSM  is  initialized  from 
the  overly  dry  upper-level  AFGWC  analysis.  More  importantly.  Figure  12  established  that  the 
vertical  profile  of  the  NH  mean  RH  of  the  GSM  forecast  agrees  fairly  closely  with  the 
observation-based  vertical  profile,  although  showing  a  modest  moist  bias  at  the  upper  levels. 

Not  surprisingly  then,  when  applied  to  the  GSM  the  AFGWC  CPS-cloud  scheme  yields  very 
large  upper-level  (and  hence  total)  forecast  cloud  amounts  showing  a  severe  positive  bias.  This 
is  illustrated  later  in  Section  3.2  and  in  the  forerunner  study  of  MW.6  The  AFGWC  CPS-cloud 
scheme  was  the  only  cloud  scheme  tested  in  the  forerunner  GSM  cloud  forecast  experiments  of 
MW.  Despite  its  poor  performance  in  that  study,  it  is  tested  again  here  for  the  sake  of  com¬ 
pleteness  and  to  serve  as  one  of  several  examples  demonstrating  a  major  point  to  unfold  in  the 
present  study  --  namely,  a  type  1  diagnostic  cloud  scheme  developed  and  tuned  for  use  in  one  NWP 
model  will  not  perform  well  when  applied  to  an  independent  NWP  model  unless  the  two  models  have 
closely  similar,  level  by  level,  humidity  forecast  characteristics  for  l)  the  domain  mean  value  and 
2)  the  frequency  distribution  of  humidity.  To  further  demonstrate  the  above  point,  we  proceed  to 
consider  and  apply  two  additional  previous  type- 1  cloud  schemes  from  other  NWP  centers. 
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The  next  such  scheme  is  that  scheme  used  operationally  by  the  ECMWF  until  May  1985. 3 
This  scheme  relates  cloud  amount.  Cl^,  at  a  given  layer  k.  to  the  layer  relative  humidity,  RHk, 

by 


CLk  = 


Max 


( RH,  - 


RHCk 


1  -  RHc, 


(1) 


v/here  RHck,  the  layer- dependent  critical  relative  humidity  for  the  onset  of  nonzero  cloud,  is 
specified  by: 


RHck=l-2ak  +  2o“  +  f3ok  (l  -  3ok  +  2<Q 


(2) 


in  which  ok=Pk/P.,  the  ratio  of  the  mid-layer  pressure,  Pk,  to  the  surface  pressure,  P.. 

Equations  (l)-(2)  assume  that  RH  and  CL  are  expressed  as  fractional  values  in  the  range  0  to  1. 
To  apply  Eqs.  (l)-(2)  in  the  cloud  postprocessor,  we  set  Pk  to  the  mandatory  level  pressure  and 
P.  to  100  kPa.  Figure  20  shows  the  RH-cloud  curves  thus  obtained  from  Eqs.  (l)-(2)  for  four 
mandatory  levels. 
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Figure  21  compares  the  critical  relative  humidities  of  the  AFGWC  and  ECMWF  schemes  as 
a  function  of  pressure.  Figure  21a  compares  the  standard  critical  value,  RHc,  for  the  zero 
cloud  threshold,  while  Figure  21b  compares  the  critical  RH  value  for  the  30  percent  cloud 
threshold.  Figure  22  compares  the  AFGWC  and  ECMWF  humidity-cloud  curves  for  the  85  and 
50  kPa  levels  over  the  full  range  of  cloud  amounts.  To  express  the  AFGWC  curves  in  terms  of 
RH  for  the  purposes  of  Figures  21  and  22,  U.S.  Standard  Atmosphere  temperatures  were  again 
used  to  convert  CPS  to  RH.  Also  represented  in  Figures  21  and  22  are  the  NMC  and  GL  cloud 
schemes  described  later. 
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Figure  21.  Critical  Relative  Humidity.  RHc,  as  a  Function  of  Pressure  for  the  Cloud  Amount 
Thresholds  of  (a)  0  Percent  and  (b)  30  Percent  for  the  ECMWF,  AFGWC,  NMC,  and  GL  Humidity- 
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In  Figure  2ia.  (be  AFGWC  scheme  vlelcls  the  smallest  RHc  values  of  all  the  four  cloud 
schemes  M  nearly  ev.-.v  c-vel  in  contrast,  the  ECMWF  scheme  generally  yields  the  largest  RHc 
values,  most  particularly  at  the  lowest  ami  best  levels.  The  type-1  cloud  scheme  In  Eqs. 
(l)-(2)  was  finalized  atiei  adjustments  to  give  reasonable  agreement  between  the  total  cloud 
forecasts  >f  the  ECMWF  model  and  observed  total  cloud  amounts  (in  terms  of  cloud  clima¬ 
tologies  and  subjecti.  ■  comparisons  with  satellite  imagery).1  One  can  speculate  that  the 
adoption  of  large  RHc  values  at  low  and  high  levels  in  the  ECMWF  scheme  was  a  response  to 
significant  moist  lama  ;  in  'he  model's  humidity  forecasts  at  those  levels.  Recall  from 
Section  2.2  that  Stingo  and  Ritter3  found  the  ECMWF  model  to  be  considerably  too  moist  In 
the  model's  lowest  a-  layer,  especially  over  the  oceans.  Further  examination  of  Figures  20-22 
shows  the  ECMWF  scheme  yields  the  largest  critical  humidities  for  virtually  all  cloud  amount 
thresholds  (30,50.70).  Not  surprisingly  then,  in  Section  3.2  the  layer  and  total  GSM  cloud 
forecasts  based  on  the  ECMWF'  scheme  will  depict  too  little  layer  and  total  cloud. 

The  large  ertth  d  homldffies  in  the  ECMWF  scheme  are  also  manifest  In  the  results  of  the 
study  by  Norquist 1  -  ?H.  Using  Eqs.  { 1 )- (2)  In  the  Inverse  sense  to  infer  humidity  estimates  for 
an  objective  humidity  analysis  from  3DNEPH  analyzed  layer  cloud  amounts,  Norquist  found 
the  ECMWF  interred  humidities  to  be  too  high  (substantial  positive  bias)  compared  with 
collocated  RAOBs  as  shown  in  Table  4  --  this  despite  the  fact  that  at  low  and  middle  levels,  at 
least  d«e  3i.)NRPH  mean  cloud  amounts  during  the  Norquist  study  period  were  shown  to  be  in 
fair  agreement  .villi  other  cloud  climatologies. 3 1 

Tito  «nf  ier  expei  Toce  illustrates  the  following  important  point:  an  empirical  humidlty- 
cioud  relationship  originally  derived  for  use  in  a  forecast  model  to  infer  cloud  amounts  from 
forecast  humidity  will  likely  nol  perform  well  in  the  inverse  sense  to  infer  humidity  estimates 
from  observed  cloud  amounts  for  use  In  an  objective  humidity  analysis.  One  reason  for  this  is 
that  ‘die  mean  values  and  frequency  distributions  of  humidity  in  model  forecasts  can  depart 
substantially  from  those  obtained  from  observations.  Such  departures  occur  because  of  1) 
difference::  in  size  and  time  span  of  the  atmospheric  volume  or  grid  cell  represented  by  the 
forecast  -er-ns  ,.>bsoj rations,  that  is,  sampling"  differences,  and  2)  erroneous  systematic 
errois  n;  (he  ‘o  recast  model. 31  For  these  same  reasons,  an  empirical  humidity-cloud  scheme 
developed  origan. My  ‘or  the  purpose  of  inferring  humidity  estimates  from  observed  cloud 
amounts  lot  use  in  an  objeclb’e  humidity  analysis,  likely  will  not  perform  well  when  inverted 
to  inlet  cloud  amounts  from  humidity  forecasts  in  a  forecast  model.  The  next  and  last 
previous  d 'agnostic  numidify  ■*  ioud  scheme  to  lie  considered  here,  namely  the  NMC  scheme,  is 

Mfll  >.U  '  XU  !■  '[  >i< 

i  ho  '.Mr  on  ii.y-..  in  ad  ■  heme  was  iiist  dt'seribed  In  the  NMC  report  by  Chu  and 
Eon  !;  at:  .  n;;T  i ,  ■  ECMWF  as  described  by  Tibairii.37  This  scheme  consists  of  two 


id  l'u!T!-.u  t ).  Il'C/V)  i  tumidity  analyses  for  operational  prediction  models  at  the 
.-■a  t  u-n.ter.  NMC  Office  Note  140.  National  Meteorological  Center, 


i'ibaidi.  5.  i . 
and  on  the  hive 
Arrona.iiiro.  42: 
Wa(t‘-  Vajnn  Il'IK 


he  ECMWF  lunidliy  analysis  and  11  s  general  impact  on  global  forecasts 
r  i  mi  the  Mediterranean  area  in  particular,  Rivasta  di  Meteorologia 
un  328,  (portion  used  in  this  report  also  available  in  the  book  Atmospheric 

’))  Academic  Press.  ISBN  0-12-208440-3). 


parts  that  respectively  infer  relative  humidity  >- <•  •  ,d  -a-  ..!  fiver  •  fmi  ;•)  o!>s  ned 
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Table  6.  Parameters  of  the  NMC/Tlbaldl  Humidity-Cloud  Scheme. 


LAYER 

LAYER 

INTERFACES 

FOR  P*=101  kPa 

(kPa) 

LEVELS  OF 

APPLICATION  IN 

GL  GSM  POST¬ 
PROCESSING 

(kPa) 

RHck  = 
RHmk  = 

Mk'Ak 

Mk  +  Ak 

Mk 

(%) 

Ak 

(%) 

RHck 

(%) 

RHmk 

(%) 

PBL 

101-96 

100 

80 

20 

60 

100 

Low 

96-74 

85 

75 

15 

60 

90 

Middle 

74-52 

70 

60 

15 

45 

75 

High 

52-30 

50,40,30 

55 

10 

45 
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To  apply  Eq.  (4)  in  the  forecast  context  of  the  GSM  cloud  postprocessor,  one  must  invert 
Eq.  (4)  and  solve  for  CL^  as  a  function  of  RHk.  which  yields 


Cl1t  =  (100/jt)  {ARCCOS[(Mk-RHk)/Ak]}< 


Mk  -  =  RHck  <  RH^  <  RHn^  =  +  Ak 


Since  the  function  ARCCOS(  a)  is  defined  only  for  the  argument  range  -1<  a<l,  the  input  value 
of  RHk  in  the  above  is  constrained  to  the  given  range.  In  practice  then,  we  default  to  CL^O 
when  RHk<RHck  and  CL^lOO  when  RHk>RHmk. 

In  the  GSM  application  here,  where  we  are  constrained  to  operate  with  forecast  RH  on  six 
mandatoiy  pressure  levels,  we  set  P.=101  kPa  in  Eq.  (3).  Table  6  shows  the  specific  layer 

interfaces  that  result  and  the  mandatory  pressure  levels  that  each  layer  spans.  To  implement 
this  scheme  in  the  cloud  postprocessor,  we  chose  to  apply  Eq.  (5)  for  a  given  layer  k  to  each 
mandatory  pressure  level  spanned  by  that  layer. 

It  is  instructive  to  compare  the  NMC  scheme  with  the  previous  AFGWC  and  ECMWF 
schemes.  A  review  of  Figures  8,  20,  and  23  shows  the  NMC  scheme  is  the  only  scheme  that 
yields  100  percent  cloud  cover  for  RHclOO,  that  is,  the  only  scheme  in  which  RHmk,  the 
critical  value  for  100  percent  overcast,  satisfies  RHmk<100.  From  a  conceptual  or  intuitive 
view,  this  is  not  necessarily  unreasonable,  as  cloud  layers  in  the  atmosphere  are  often  quite 
thin  and  thus  vertically  "subgrid-scale"  with  respect  to  the  larger  layer  thicknesses  typical  in 
forecast  or  objective  analysis  models.  Equivalently  stated,  thin  cloudy  layers  can  be  imbedded 
in  thicker  atmospheric  layers  whose  vertically  averaged  RH  is  below  100  percent.  However, 
the  degree  to  which  RHmk  falls  below  100  percent  in  Table  6  is  disturbingly  large,  especially 

for  the  middle  and  high  layers.  In  this  regard  it  is  not  surprising  that  in  the  results  of  the 
study  by  Norquist  in  Table  4,  the  bias  of  the  humidities  inferred  from  the  NMC  scheme  is 
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negative  (;>ll  ■■■■  :  .;  •  *.*.  ad  -  •  - ;  .  ••(!  ,'i,  .  .*  ihe  ! it  k-s  iniened  at  middle 

and  upper  1«  veis  *.i om  da-  NMC  v- ■!:  ■  .<_•  n  ■-.y.ii-vA  7:>  and  (>f>  percent  respectively. 
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3.1.3  THE  NEW  GL  DIAGNOSTIC  CLOUD  SCHEME 


The  disappointing  performance  of  the  previous  diagnostic  cloud  schemes  (demons.;!. tied 
later  in  Section  3.2)  led  to  the  conclusion  that  a  diagnostic  humidilv  cloud  scheme  tor  a  t*iven 
forecast  model  must  be  developed  specifically  lor  that  model.  The  develoumen:  .1  h  n<  ■■■■_.■ 
must  account  for  two  major  properties  ol  die  iorecast  model  s  ho miditv  fields:  .•'t.mv.n- 

mean  humidity  and  b)  the  frequency  distribution  of  humidity  in  parti*,  uko  8.  :iou  :c. 
demonstrated  that  these  two  humidity  properties  t i J  change  significamh  'A  ring 
spin-up  period.  (2)  achieve  by  the  end  oi  the  spin-up  period  iairlv  tble  mode.  t.  ••  o  •  simes 
that  are  largely  independent  of  initial  conditions,  (3.1  vary  with  height  md  d  •  •.>.  . 
geographic  region  (especially  land  versus  ocean  at  low  level*!,  '[]  •  ,  a  -o 

here  will  be  developed  specifically  ior  the  baseline  GL  GEM  and  w4  ;  ,  uc 

factors  (l)-(3)  above.  This  new  scheme  can  also  account  for  die  t..  u  s.  c  a  ' 

of  forecast  cases  available  to  tills  study  was  too  limited  to  a-xrcive  d.i;  jv-c.-nt'... 

An  instructive  approach  to  describing  the  scheme  is  to  t  [  it.  U:  I  «.  i  .  '.V-iJS  ,  r  *  i  F 

led  to  its  development.  Of  the  previous  cloud  schemes,  the  ECMWI-  s,  he  >r  .<■  ■ 

give  the  best  (though  still  disappointing)  cloud  forecast  results  when  tea  :  with  ;u*  baseline  cl 
GSM.  This  result  led  to  a  number  of  attempts  to  adjust  or  turn-  the  ::  U.VAYE  ■-*' a.-me.  i  ■■■,;■  •uis 
purpose,  one  can  express  Eq.  (1)  in  a  more  general  form  as  follows: 


I  RH,  -  R1  ic, 

CLuJ  k 


RHm,  -  RHc. 
k  k 


.  0  <  RHc  <  RH,  r  RHm. 


where  k  Is  the  level  Index.  RHck  and  RHm,.  are  the  /em-eloml  wt 
and  Pk  is  the  exponent  of  nonlinearity  (die  special  case  ot  f  • 

serves  to  illustrate  that  there  arc  at  least  tour  d<  gree>  oi  hvedom 
cloud  curves,  namely,  RIIck.  Rllmk.  Pk.  and  k  (vertical  level;  On 
degrees  of  freedom  wheiein  i’k  varies  within  (lie  ramie  oi  Rtnq  • 

changes  in  nonlinearity  in  Figure  8)  and  all  parameters  van  as  c 
in  response  to  temporal  changes  in  the  statistical  character  d  < 

The  initial  tuning  attempts  with  the  ECMWI-'  s,  heme  wen-  oi 
nature  and  aimed  only  at  decreasing  the  RH* y  values  In  the  EC  M 

substantial  negative  bias  in  the  cloud  forecasts  with  that  sihcmt  ; 
verifications  from  one  somewhat  sun  cssl'u!  such  trial  m-  i aiwm 
compares  the  RHc  values  for  the  original  and  modne-d'  '  CMWF 
and  error  modifications  to  the  ECMWF  scheme  also  tested  .  u,i*!r 
Eq.  (6)  (for  example,  P,  --  1,  2.5.  4).  but  without  non  worthy  mi 
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Table  7.  Comparison  of  Critical  Relative  Humidity  Values, 
RHc,  In  Original  and  Modified  ECMWF  Humldty-Cloud  Scheme 
(assume  P*  =  100  kPa). 


p 

(kPa) 

ECMWF 

ORIGINAL 

RHc 

ECMWF 

MODIFIED 

RHc 

30 

0.73 

0.65 

40 

0.61 

0.55 

50 

0.50 

0.50 

70 

0.44 

0.40 

85 

0.61 

0.55 

100 

1.00 

0.85 

Because  the  number  of  degrees  of  freedom  In  Eq.  (bj  Is  formidable,  we  sought  objective 
approaches  to  specifying  them.  We  were  crucially  aided  here  by  the  fact  that  the  KTNEPH  Is  an 
objective,  gridded,  global  3-D  cloud  analysis  database  (as  opposed  to  merely  a  satellite  image  or 
a  zonal  average  cloud  climatology).  The  key  objective  application  here  was  to  calculate  and 
compare  frequency  distributions  of  1)  RH  from  the  GSM  forecasts  and  2)  layer  cloud  amounts 
from  the  verifying  RTNEPH  cloud  analysis.  The  top  of  Figure  24b  shows  as  an  example 
frequency  distribution  for  the  GL  GSM  RH  forecast  at  70  kPa  over  the  NH  octagon  for  the 
24-hour  forecast  initialized  from  the  AFGWC  analysis  of  00Z,  17  January  1985.  The  bottom  of 
Figure  24b  shows  the  corresponding  cumulative  frequency  distribution  of  RH  (that  is,  for  a 
given  RH  value  r,  the  percent  frequency  of  RH  values  in  the  range  0  <  RH  <  r).  The  top  of 
Figure  24a  shows  the  cloud  frequency  distribution  from  the  verifying  RTNEPH  70  kPa  layer 
cloud  analysis  for  the  NH  octagon.  The  bottom  of  Figure  24a  shows  the  corresponding 
cumulative  cloud  frequency  distribution.  Both  the  RH  and  cloud  distributions  in  Figure  24 
were  computed  in  terms  of  101  categories  at  1 -percent  intervals  from  0  to  100. 
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Figure  24.  Category  (top)  and  Cumulative  (bottom)  Frequency  of  Occurrence  (at  one  percent 
intervals)  of  RTNEPH  Analyzed  Cloud  Amounts  (Part  a)  and  GSM  24-hour  Forecast  RH 
Amounts  (Part  b)  at  70  kPa  Over  the  NH  Octagon.  Valid  at  00Z.  18  January  1985.  At  bottom, 
the  one-to-one  mapping  of  the  RTNEPH  observed  cloud  frequency  onto  the  GSM  forecast 
RH  frequency  yields  one  of  the  GL  humidity-cloud  conversion  curves  (solid)  in  Figure  25a. 


The  first  use  of  these  RH  and  cloud  frequency  distributions  at  each  level  was  to  determine 
an  appropriately  tuned  critical  value  of  RHek,  obtained  objectively  by  "mapping"  the  value  of 

the  zero-cloud  cumulative  frequency  to  the  same  cumulative  frequency  value  in  the  RH 
distribution,  as  illustrated  by  the  lowermost  horizontal  directional  arrow  in  Figure  24.  A 
similar  mapping  of  the  cumulative  frequency  values  for  each  1 -percent  cloud  interval  onto  the 
RH  cumulative  frequency  yields  a  complete  quasi-eontinuous  empirical  humidity-cloud  curve 
relating  the  entire  range  of  fractional  cloud  cover  to  an  appropriate  range  of  RHck  to  RHmk  --  a 

procedure  that  can  be  repeated  level  by  level  over  the  desired  GSM  forecast  intervals.  We  shall 
henceforth  call  this  procedure  the  GL  humidity-cloud  scheme. 

Figure  25a  shows  the  explicit  70  kPa  RH-cloud  curve  (solid  curve)  that  results  from  the 
mapping  illustrated  in  the  bottom  of  Figure  24.  The  dashed  curve  in  Figure  25a  shows  the 
corresponding  70  kPa  curve  that  results  from  the  24-hour  GSM  forecast  (and  verifying 
RTNEPH  layer  cloud  analysis)  initialized  from  the  00Z,  24  January  1985  AFGWC  analysis. 

The  dotted  curve  in  Figure  25a  shows  the  "aggregate”  24-hour  70  kPa  RH-cloud  curve  that 
results  from  first  combining  the  individual  samples  of  the  70  kPa  17  and  24  January  GSM 
24-hour  RH  torecast  and  then  verifying  RTNEPH  layer  cloud  analysis  into  a  single  aggregate 
sample,  from  which  the  requisite  cumulative  frequency  distributions  of  RH  and  cloud  is  then 
computed.  Figure  25b  shows  the  corresponding  three  curves  derived  from  the  24-hour  GSM  RH 
forecasts  and  layer  RTNEPH  cloud  analyses  at  the  30  kPa  level.  As  expected,  the  aggregate 
curves  fall  between  the  two  separate  curves  for  the  Individual  forecast  cases. 
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Figure  25.  Curves  of  the  GL  Humiditv-clouci  Conversion  Scheme 
Derived  for  the  GSM  24-hour  Forecasts  at  (a)  70  kPa  and  (b)  30  kPa 
by  Applying  the  Methodology  of  Figure  24  to  Two  Separate  Samples 
(and  their  aggregate)  of  Uie  Frequency  Distributions  of  GSM  Forecast 
RH  and  RTNEPH  Verifying  Cloud  Over  the  NH  Octagon. 
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In  practice,  the  GSM  forecast  RH  was  converted  to  Rtf'.  Then  the  calculation  of  frequency 
distributions,  mapping  of  those  distributions,  and  resulting  humidity-cloud  curves  were  all 
defined  in  terms  of  RH*.  However,  all  related  plots  pieseru.v  Vm  in  this  study  are  depicted  to 
runs  of  the  normal  RH,  sim  «*  the  trader  t*  I  <uit •••#■'•.  i-.jon  »m  >vlth  its  range  and  since 
It  :  takes  it  easic  to  '-•>■;  with  ptev'-v-..  '.I.  humidity-cloud  seller ’« 

.vlth  tlie  GSM,  this  study  found  the  1  '  ;.,t>  j  .-•ore.*  '■<  he  slightly  better  usm,. 

'he  curves  derived  ftoni  RH*  than  using  those  ■  t  •  1  '•<  •  FH.  Given  only  two  forecast  cases, 
this  slight  skill  improvement  may  be  foituitous.  Nevertheless,  we  proceeded  to  use  RH* 
routinely.  All  the  GL  RH -cloud  curves  presented  to  this  ropr> '  were  derived  using  RH*- 
We  next  extend  the  examples  of  the  GL  <  loud  scheme  io  utn-c  km-mst  times  and  the 
remaining  levels,  The  real  va'ue  of  (he  Ol  cloud  scheme  is  It  ability  1°  respond  and  adjust  t< 
the  GSM  spin  up  period.  Hv  /'•prating  the  piocedme  ■ : i  ;  .me  24  for  the  GSM  RH  frequency 
distribution  valid  for  each  t-SM  forecast  Inten  ,i  (12 •  hourly  hen  .  one  can  derive  RH-cloud 
curves  that  vary  as  a  luncHon  of  GSM  forecast  i-ugth  P'liiir  hows  the  12-hourly  70  kPa 
arid  30  kPa  RH-cloud  cm  »•»<•  thus  oblaMK-d  .  ,  .-pec ted.  ih-  hi,  >vst  change  occurs  in  tire  first 
■  ?  hours,  but  by  48-licm  ei m *-s  h  •-,-••!  umhea1'-  .•.•cited  a  fairly  constant 
•emporal  confirmation  -king  th-  ♦-•>>  1  o  •••.  ,  <ealtngly.  the  early  chanr* 

in  the  70  KPa  •  u  es  u  I  igm  '  '’(>;•  i  >  p.  <  i  nonlinearity  of  the 

PH  i  cquenev  distribution 
■  •  "s  by  only  a  modest  6 
1  value  at  70  kPa 
■  urves  in  Figure  26b  N 
*•  to  the  rapid 
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Figure  26.  Curves  of  the  GL  Humidity-cloud  Conversion  Scheme 
as  a  Function  of  GSM  Forecast  Length  for  (a)  70  kPa  and  (b)  30  kPa 
Mandatory  Pressures.  (Curves  based  on  the  aggregate  NH  octagon 
GSM  RH/RTNEPH  cloud  frequency  distributions  for  the  00Z, 

17  and  24  January  1985  GSM  forecast  cases). 
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Figure  26  offers  a  vivid  demonstration  of  the  GL  humidity-cloud  scheme  adaptability  to  the 
temporal  changes  in  the  gross  character  of  the  GSM  RH  forecasts.  As  demonstrated  In 
Section  3.2,  a  chief  outcome  of  this  adaptability  is  the  capability  of  the  GL  scheme  to  provide 
short-range  total  cloud  forecasts  whose  octagon  domain-mean  cloud  amounts  are 
substantially  Invariant  in  time  during  the  GSM  spin-up  period  and  whose  octagon  bias  with 
respect  to  the  verifying  RTNEPH  total  cloud  analysis  is  small.  (What  total  cloud  bias  there  is 
enters  from  the  vertical  stacking  algorithm  of  Appendix  A,  as  the  GL  scheme  effectively 
guarantees  that  the  octagon  domain-mean  layer  cloud  forecast  will  have  virtually  zero  bias 
against  the  RTNEPH.  The  occurrence  of  some  total  cloud  bias  prompted  the  vertical  stacking 
experiments  presented  later  in  Section  3.2,  as  well  as  the  variant  of  the  GL  cloud  scheme 
described  next.) 

Figure  27  shows  the  GL  RH-cloud  curves  obtained  for  the  O-hour  and  48-hour  GSM  RH 
forecasts  at  all  six  moist  mandatory  levels.  It  must  be  reiterated  here  that  this  study  utilized 
RTNEPH  fields  archived  In  the  spatially  compacted  format  of  the  5LAYER  grid  (see  Figure  7), 
that  is,  on  the  1 /2-mesh  NH  octagon  grid  at  the  vertical  levels  (layers)  of  30,  50,  70,  and  85  kPa 
and  the  gradient  level.  Thus  to  derive  the  GL  RH-cloud  curves  for  the  GSM  levels  of  40  and  100 
kPa  In  Fig.  27,  we  had  to  accept  some  mismatch  between  these  two  GSM  output  levels  and  the 
RTNEPH  compacted  (5LAYER)  levels.  As  a  compromise,  we  used  the  RTNEPH  cloud  frequency 
distributions  at  the  gradient  and  30  kPa  levels  to  derive  the  100  and  40  kPa  GSM  RH-cloud 
curves,  respectively. 
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i><uie  '?,!  Curves  of  tin  tjl.  Humidity-cloud  Conversion  Scheme 
,  ;  ,1  i  •'!)•  !!  •;  >  '  •  ■i'  ’  v  (‘p-.Mii «•  lor  tlie  GSM  O-hour  and 
-1^  Mc»i!  !  o.  "casts  'based  on  aggregate  sample  as  In  Figure  26). 
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The  above  "pairings"  dilemma  was  an  unavoidable  consequence  of  the  structure  of  the  data 
bases  at  our  disposal.  In  a  follow-on  study,  which  wculd  anticipate  the  development  of  wh 
cloud  curves  following  the  GL  procedure,  one  would  prefer  to  work  with  GSM  RH  toreca  '  c 

directly  on  GSM  o-layers  and 


(1)  ignore  the  lowest  o-layer  (much  too  moist  over  oceans  as  cited  in  Section  2.2), 


(2)  group  the  remaining  six  moist  GSM  o-layers  into  three  pairs  corresponding  to  low, 
middle,  and  high  cloud  layers. 


(3)  choose  the  maximum  RH  in  each  pair  of  o-layers  to  represent  the  GSM  RH  forecast  for 
the  combined  layer,  and  finally 


(4)  map  the  RH  frequency  distributions  from  these  low,  middle,  and  high  GSM  layers  to  the 
frequency  distributions  for  low,  middle,  and  high  cloud  layers  obtained  from  a  suitably 
matched  vertical  compaction  of  the  original  RTNEPH  layers 


The  modest  real  vertical  resolution  of  the  RTNEPH  and  the  high  vertical  correlation  m  he  RH 
forecasts  of  adjacent  GSM  layers  (as  shown  in  Section  2.2)  does  not  encourage  more  than  the 
three  layers  recommended  here. 

The  present  study  could  not  follow  the  above  procedure.  Nevertheless,  an  appropriate 
rendition  was  pursued  by  developing  a  four-layer  variant  of  the  six  curves  in  Figure  27  These 
four  curves  are  given  in  Figure  28.  The  50  and  70  kPa  curves  in  Figure  28  are  unchanged  fr  om 
those  in  Figure  27.  The  RH-cloud  curve  designated  high"  in  Figure  28  was  developed  by 
mapping  the  frequency  distribution  of  the  RTNFPH  30  kPa  cloud  analysis  to  the  Rid  frequency 
distribution  obtained  by  taking  the  maximum  of  the  30  and  40  kPa  GSM  RH  forecast  values  at 
each  grid  point. 
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Figure  28.  Curves  of  the  Four-layer  Version  of  the  GL 
Humidity-cloud  Conversion  Scheme  (see  text)  as  a  Function 
of  Mandatory  Pressure  for  the  GSM  O-hour  and  48-hour 
Forecasts  (based  on  aggregate  sample  as  in  Figure  26). 
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The  RH-cloud  curve  designated  "low"  in  Figure  28  was  developed  by  taking  the  maximum 
RTNEPH  cloud  amount  of  the  gradient  and  85  kPa  levels  at  each  grid  point  and  mapping  the 
resulting  cloud  distribution  to  the  RH  distribution  of  the  85  kPa  GSM  RH  forecast.  The  100 
kPa  GSM  RH  forecast  is  purposely  ignored,  because  Figures  16  and  18  show  it  suffers  from  a 
severe  moist  bias  over  oceans.  The  vertical  cloud  stacking  algorithm  for  the  four-layer  GL 
cloud  scheme  is  described  in  Appendix  A.  The  six-layer  and  four-layer  GL  cloud  schemes  will 
henceforth  be  designated  as  GL6  and  GL4,  sometimes  with  an  additional  subscript  to  denote  a 
specific  GSM  forecast  length  in  hours  (for  example,  GL624  or  GL424).  The  performance  of  both 

the  GL6  and  GL4  schemes  will  be  examined  in  Section  3.2. 

In  concluding  this  section,  it  is  useful  to  consider  the  GL  humidity-cloud  scheme  as 
represented  in  Figure  27  in  relation  to  the  free  parameters  introduced  earlier  in  presenting 
Eq.  (6).  Figure  27  (and  its  counterparts  for  the  other  GSM  12-hourly  forecast  times,  not  shown) 
explicitly  specifies  RHck  and  RHmk  (the  latter  always  turns  out  to  be  100  percent)  as  a  function 
of  level  k  and  forecast  length.  Furthermore,  the  nonlinearity  exponent,  Pk,  and  its  variation 
within  the  range  RHck-  RHmk,  though  not  explicitly  obtained  by  the  GL  procedure,  is 
nonetheless  implicitly  embodied  in  the  specification  of  the  GL  curves;  that  is,  in  the  six  (GL6) 
or  four  (GL4)  look-up  tables  defined  for  each  desired  GSM  forecast  length  in  terms  of  101  one- 
percent  cloud  categories.  Thus  the  GL  scheme  embraces  all  six  degrees  of  freedom  cited  in  the 
discussion  following  Eq.  (6). 

To  conclude  this  section,  we  examine  in  Figures  2 1  and  22  the  comparison  of  the  GL  RH- 
cloud  scheme  to  the  three  previous  diagnostic  cloud  schemes  of  Section  3.1.2.  These  figures 
show  the  GL6  scheme  at  the  24-hour  forecast  time  (GL624),  which  is  after  the  majority  of  the 

GSM  spin-up  is  completed  (at  least  on  the  extra  tropical  octagon  domain).  The  zero-cloud 
critical  humidity  profile  of  the  GL  scheme  shows  moderate  values  that  do  not  vary  drastically 
with  height.  These  values  lie  between  those  of  the  AFGWC  and  ECMWF  schemes  and  fairly 
closely  follow  the  profile  of  NMC  zero-cloud  critical  humidity  values.  Figure  22  shows  that,  in 
agreement  with  the  AFGWC  and  ECMWF  schemes  (but  unlike  the  NMC  scheme),  the  overcast- 
cloud  threshold  of  the  GL624  scheme  is  RHmk=100. 

Rather  interestingly  in  Figure  22,  the  overall  shape  of  the  GL624  curves  is  fairly  similar  to 
that  of  the  ECMWF  curves,  and  thus  reasonably  described  oy  a  quadratic  function.  Thus  after 
the  GSM  spin-up  at  forecast  times  of  24  hours  and  beyond,  the  GL6  curves  could  be  analytically 
approximated  by  Eq.  (6)  with  RHmk=l,  Pk=2,  and  RHck  as  given  for  the  GL624  scheme  in 

Figure  21a  (normalized  by  100).  Furthermore,  it  is  interesting  to  compare  the  objectively 
derived  values  of  RHck  for  the  GL624  scheme  in  Figure  21a  to  the  subjectively  (trial  and  error) 
derived  values  of  RHck  for  the  "modified"  ECMWF  scheme  in  Table  7.  The  fixed  analytical 
representation  of  Eq.  (6)  was  not  further  pursued  here  for  the  GL  schemes,  as  we  wanted  to  tap 
the  full  capability  of  the  GL  schemes  to  adjust  to  the  short-range  GSM  spin-up  behavior. 
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fhe  previous  section  examined  and  compared  the  formulations  of  six  diagnostic  cloud 
.->•  .  :  x  ltd  to  as  the  .u-YWC,  NMC,  ECMWF,  modified  ECMWF,  GL6,  and  GL4  schemes, 

i  ids  section  examines  the  cloud  forecast  performance  of  these  six  schemes  when  applied  to  the 
humidity  forecasts  of  the  baseline  GL  GSM.  The  GSM  results  from  these  six  schemes, 
respectively,  shall  be  denoted  in  tables  and  figures  ac  GSM/A,  GSM/N,  GSM/E,  GSM/EM, 
GSM/GLG,  and  GSM/GL4.  Cloud  forecast  performance  will  be  examined  quantitatively  and 
subjectively  for  two  NH  winter  GE  GSM  forecast  cases  --  namely,  the  GSM  forecasts  initialized 
from  the  opera  Ilona1  \1  GWC  global  analyses  for  00Z,  17  January  1985  (designated  Case  I)  and 
00  Z,  24  January  1985  (designated  Case  ll).  As  a  standard  for  comparison,  this  section  also 
examines  the  cloud  forecast  performance  of  AFGWC's  operational  5LAYER  model  for  Cases  I 
th!  F  V:  backup  -  ind  for  the  GSM  cloud  forecasts.  Sections  2.1.2  and  2.2  »»ar!i»t 

evnutncd  tbc-  rot  responding  relative  humidity  analyses  and  forecasts  for  Case  I. 

The  quantitative  performance  of  the  GSM  and  5LAYER  total  cloud  forecasts  (as  well  as 
RTNEPH  persistence  -  denoted  PERS)  are  measured  here  following  AFGWC  operational 
practice:  that  Is,  by  polutwise  grld-to-grld  verifications  against  the  AFGWC  operational 
WiNEPH  total  cloud  analysis.  This  verification  is  performed  on  the  93  X  101  AFGWC  NH 
i  ,>.n.c5!»1  >c  !a*;  'u  grid  This  grid  Is  1)  a  one-to-one  subset  of  the  129  X  129  1/2-mesh 

Jap!  'ilc  output  gr  !  of  the  GSM  cloud  postprocessor  (Section  3.1.1)  and  2)  the  grid  of  the 

<  ;  ti;  *f  ’-d  i  d'  impacted  RTNEFH  cloud  analysts,  produced  by  the  spatial  compaction 

st;..  •?.'<■■  ■'  ■  (Jie  left  ,1  Figure  7. 

More  r.otably  t  his  grid  is  also  the  operational  grid  of  the  5 LAYER  model.  We  recall  from 
Figure  7  that  'he  compacts  d  Ri'NEPH  cloud  analysis  is  Identical  to  the  5LAYER  0-hour  initial 
’> -  i:d  lb'  eh  iayc  r.nd  (o‘al)  Hence  by  definition,  the  O-liour  5 LAYER  cloud  forecast  ' 

tfi>  .an  s  will  b  perfect.  Even  beyond  the  0-hour,  one  might  argue  that  performing 

die  verification  directly  on  the  SLAYER  grid  will  naturally  favor  5 LAYER  over  GSM  forecasts. 
However,  despite  the  area-averaging  In  the  cited  compaction  step,  the  resolution  of  the 
omp,  cv.-d  RTNEPH  used  tor  the  verification  is  substantially  higher  than  the  2.5° 

j.c;lu,d.  s  v  ’pot  grid  of  the  AFGWC  GSM  Hence  in  considering  the  replacement  of 
m  Y!"  cf'fiu  t'o-cc  .s.s  with  GSM  cloud  forecasts,  one  must  demonstrate  the  GSM 
'  '  1  '  •  i  ■  i  :  ••■'r-rc...olu.lio: i  framework  used  to  routinely  assess  5LAYER  skill, 

s  <i.  ■  ip;  •  drg  steps  of  Section  3. 1 . 1  included  an  interpolation  of  the 

<  ■■  T  .  1  !ude  gild  to  1/2-mrsb. 

’in  j  ci.'  i  «<«'  ?lJcs  obtained  here  will  include  1)  the  bias  in  percent  cloud 

i.  •  -*  - ■  ’•  'c;  •  oiv  fin.wii  as  th  ri.!/ 20  score.  This  score  is  the  traditional  AFGWC 

r_r:  id  i  >o  .  ;  :ri,u  de-tiue  -i  as  the  non  enf  r<~  grid  oofiif-  where  the  forecast  of 

r  differ  by  less  than  7n 

percent.  A  pence t  20/20  scon  is  loO  percent,  the  bias  ana  20/20  score  will  be  derived  only 
for  the  total  cloud  lorecaots,  as  tin’  point  wise  reliability  of  the  verifying  RTNEPH  eleod 


analysis  is  substantially  better  for  total  cloud  than  layer  cloud.  Both  total  and  layer  cloud 
forecasts  will  be  examined  here  subjectively  from  example  displays 

All  of  the  diagnostic  cloud  schemes  detailed  in  Section  S.t  provide  only  layer  cloud 
forecasts  at  the  moist  mandatory  pressure  levels  of  the  GSM  output.  To  obtain  total  cloud 
forecasts,  the  layer  cloud  forecasts  are  passed  through  a  ver'ical  stacking  algorithm  to  obtain 
a  total  cloud  amount.  The  vertical  cloud  stacking  is  described  in  detail  in  Appendix  A.  The 
reader  is  urged  to  review  Appendix  A,  as  several  stacking  issues  wail  be  addressed  in  this 
section. 

We  shall  center  the  examination  of  cloud  forecast  periotma.  .  •.»  on  Tables  P  '1  rough  '  i 
Table  9  shows  the  mean  (area- averaged)  total  cloud  amount  as  a  function  of  forecast  length  tn 
Case  I  for  the  NH  octagon  domain  for  ten  GSM  cloud  schemes,  plus  the  51  AYER  model  ant!  the 
verifying  RTNEPH  analysis.  For  Case  I  and  Case  II,  respectively,  Tables  10  and  1  1  show  the 
total  cloud  forecast  bias  as  a  function  of  forecast  length  for  the  NH  octagon  ana  the  regions 
denoted  "U.S."  and  "Europe"  In  Figure  1  for  the  same  ten  GSM  cloud  schemes,  plus  the  5LAVER 
forecast  and  RTNEPH  persistence. 


HI 


Table  9.  Mean  Total  Cloud  Amount  (percent)  as  a  Function  of 
Forecast  Period  Over  the  N.H.  Octagon  Domain  for  00Z, 

17  January  1985. 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24-HR 

36- HR 

. -  - 

48-HR 

Verification 

(RTNEPH) 

48.0 

50.2 

45.1 

47.0 

42.7 

pnn 

48.0 

48.0 

48.0 

48.0 

48.0 

5 LAYER 

48.0 

47.4 

46.0 

43.1 

41.6 

GSM/A 

58.8 

74.4 

79.9 

msm 

78.2 

GSM/N 

59.6 

71.4 

74.3 

72.8 

73.0 

GSM/E 

19.5 

32.8 

32.7 

31.2 

31.6 

GSM/EM 

21.6 

35.3 

36.0 

35.0 

35.4 

GSM/GL6 

50.9 

54.0 

54.4 

52.1 

52.9 

GSM/GL6A 

55.2 

57.9 

58.1 

56.0 

56.6 

GSM/GL6B 

38.8 

43.1 

42.2 

40.9 

39.7 

GSM/GL6C 

43.8 

47.6 

46.5 

45.3 

43.8 

GSM/GL4 

44.5 

46.9 

46.8 

44.6 

43.1 

GSM/GL4A 

49.1 

51.3 

50.8 

48.9 

47.2 
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Table  10.  Total  Cloud  Forecast  BIAS  (In  Percent  Cloud  Cover) 
for  10  Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  SLAYER 
Model,  and  KTNEPH  Persistence  for  the  00Z,  17  January  1985 
Case  Over  the  Three  Areas  of  (A)  NH,  (B)  U.S.,  and  (C)  Europe. 


(A)  NH _ BIAS _ 00Z,  17  January  1985 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24-HR 

36-PIR 

48-HR 

Persistence 

0.0 

-2.3 

2.2 

0.8 

4.1 

5LAYER 

[0.0] 

■ 

1 

[1.1] 

-5.5 

[-4.4] 

GSM/A 

10.8 

24.7 

33.1 

30.9 

34.9 

GSM/N 

11.7 

21.4 

28.1 

25.8 

30.3 

GSM/E 

-28.5 

-17.0 

-13.4 

-15.0 

-10.1 

GSM/EM 

-26.4 

-14.4 

-10.0 

-11.2 

-6.6 

GSM/GL6 

[2.9] 

[4.5] 

8.3 

5.7 

9.9 

GSM/GL6A 

7.2 

8.4 

12.1 

9.7 

14.2 

GSM/GL6B 

-9.2 

-6.6 

1-4.0] 

-5.5 

1-1.2] 

GSM/GL6C 

[-4.2] 

[-2.0] 

[0.3] 

[-0.9] 

[3.6] 

GSM/GL4 

[-3.5] 

[-2.5] 

[0.6] 

[-1.4] 

12.6] 

GSM/GL4A 

[1.1] 

■ 

I 

[4.7] 

[2.9] 

7.3 

Percent  of 
Timely 
Area 

N/A 

84 

77 

86 

69 

Brackets  denote  values  less  than  5  percent  in  magnitude  (Part  A  only). 
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Table  10.  Total  Cloud  Forecast  BIAS  (In  Percent  Cloud  Cover) 
for  10  Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER 
Model,  and  RTNEPH  Persistence  for  the  00Z,  17  January  1985 
Case  Over  the  Three  Areas  of  (A)  NH,  (B)  U.S.,  and  (C)  Europe. 
(Continued). 


(B)  U.S. _ BIAS _ 00Z,  17  January  1985 


PERIOD 

CLOUD 

SCHEME 

00  HR 

12-HR 

24-HR 

36- HR 

48- HR 

Persistence 

0.0 

-2.0 

2.6 

7.2 

8.8 

5  LAVE  t- 

0.0 

mm 

-5.5 

-1.0 

-6.7 

GSM/A 

-4.0 

14.7 

22.0 

31.8 

31.3 

GSM/N 

-2.0 

12.3 

18.0 

26.9 

27.5 

| 

GSM/E 

-40.0 

-23.8 

-19.6 

-13.2 

-10.1 

GSM/EM 

-38.2 

-21.3 

-16.8 

-9.2 

-6.8 

GSM/GL6 

-12.9 

-2.7 

-0.2 

8.1 

7.7 

GSM/GL6A 

-8.7 

1.0 

3.3 

12.1 

11.5 

i 

GSM/GL6B 

-24  6 

-14.7 

-10.8 

-3.5 

-1.3 

GSM/GL6C 

-19.7 

-10.2 

-6.8 

1.1 

3.0 

GSM/GL4 

-17.2 

-11.1 

-7.5 

-1.3 

0.7 

GSM/GL4A 

-12.8 

-6.7 

-3.8 

3.1 

4.7 

Percent  of 
Timely 
Area 

N/A 

90 

80 

94 

62 

Table  10.  Total  Cloud  Forecast  BIAS  (In  Percent  Cloud  Cover) 
for  10  Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER 
Model,  and  RTNEPH  Persistence  for  the  00Z,  17  January  1985 
Case  Over  the  Three  Areas  of  (A)  NH,  (B)  U.S.,  and  (C)  Europe. 
(Continued). 


(C)  Europe  BIAS  00Z,  17  January  1985 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24-HR 

36-HR 

48-HR 

Persistence 

0.0 

-1.2 

-3.0 

-1.2 

-0.2 

SLAYER 

0.0 

2.2 

8.6 

-2.8 

2.3 

GSM/A 

5.1 

14.3 

22.4 

20.6 

23.8 

GSM/N 

5.0 

9.6 

18.9 

15.8 

17.7 

GSM/E 

-37.6 

-23.2 

-17.6 

-18.6 

-18.1 

GSM/EM 

-35.5 

-20.8 

■ 

I 

-16.0 

15.3 

GSM/GL6 

-2.6 

-4.6 

-2.5 

3.9 

3.8 

GSM/GL6A 

1.7 

-1.3 

1.7 

0.1 

0.4 

GSM/GL6B 

-19.1 

-14.2 

-7.8 

-9.7 

-9.8 

GSM/GL6C 

-13.6 

-10.2 

-3.1 

-5.2 

-5.1 

GSM/GL4 

-13.9 

-9.9 

-3.7 

-5.5 

-5.9 

GSM/GL4A 

-8.4 

-5.9 

0.7 

-1.5 

-1.3 

Percent  of 
Timely 
Area 

N/A 

79 

71 

85 

69 

85 


Table  1 1 .  Total  Cloud  Forecast  BIAS  (In  Percent  Cloud  Cover) 
for  10  Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER 
Model,  and  RTNEPH  Persistence  for  the  00Z,  24  January  1985 
Case  Over  the  Three  Areas  of  (A)  NH,  (B)  U.S.,  and  (C)  Europe. 


(A)  NH BAS 00Z^24JanuarAM985 


PERIOD 

CLOUD 

SCHEME 

00-HR 

12-HR 

24-HR 

36- HR 

48-HR 

Persistence 

0.0 

-3.7 

CO 

1 

-4.9 

-6.3 

5LAYER 

[0.01 

[-4.21 

-7.4 

-7.9 

-9.8 

GSM/A 

17.8 

27.0 

26.7 

30.4 

28.0 

GSM/N 

18.1 

24.3 

22.6 

26.0 

24.0 

GSM/E 

-22.6 

-16.2 

-17.2 

-15.6 

-14.9 

GSM/EM 

-20.6 

-13.6 

-14.2 

-11.8 

-11.5 

GSM/GL6 

10.4 

5.8 

[2.8] 

5.5 

[3.4] 

GSM/GL6A 

14.5 

10.2 

6.8 

9.8 

7.4 

GSM/GL6B 

1-2.3) 

[-4.9] 

-7.6 

-5.1 

-6.5 

GSM/GL6C 

(2.61 

[0.2] 

1-2.8] 

[-0.3] 

[-2.1] 

GSM/GL4 

13.6] 

[0.0] 

[-2.2] 

1-1.5] 

[-2.9] 

GSM/GL4A 

8.0 

[4.7] 

[2.2] 

[3.1] 

[1.4] 

Percent  of 
Timely 
Area 

N/A 

88 

62 

85 

73 

■  i  . . .  ...  i  ■  .  ■■■  -  - - « 

Brackets  denote  values  less  than  5  percent  in  magnitude  (Part  A  only). 
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Table  1 1 .  Total  Cloud  Forecast  BIAS  (In  Percent  Cloud  Cover) 
for  10  Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER 
Model,  and  RTNEPH  Persistence  for  the  00Z,  24  January  1985 
Case  Over  the  Three  Areas  of  (A)  NH,  (B)  U.S.,  and  (C)  Europe. 
(Continued). 


(B)  U.S.  BIAS  OOZ,  24  January  1985 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24- HR 

36- HR 

48-HR 

Persistence 

0.0 

-0.4 

3.3 

5.2 

9.6 

5LAYER 

0.0 

-7.9 

-6.5 

-10.4 

-7.5 

GSM/A 

-5.1 

10.2 

17.8 

24.5 

27.5 

GSM/N 

-5.9 

6.6 

15.6 

20.1 

21.8 

GSM/E 

-43.9 

-31.4 

-21.2 

-23.0 

-18.4 

GSM/EM 

-42.3 

-29.0 

-18.4 

-19.0 

-15.0 

GSM/GL6 

-12.8 

mm 

-3.3 

-1.3 

2.8 

GSM/GL6A 

-9.4 

-7.3 

1.1 

3.1 

6.7 

GSM/GL6B 

-28.7 

-21.4 

mm 

-12.0 

-10.1 

GSM/GL6C 

-24.3 

-16.8 

-6.7 

-7.1 

-6.0 

GSM/GL4 

-19.5 

-15.5 

-6.3 

-8.4 

-6.3 

GSM/GL4A 

-15.9 

-11.1 

-1.8 

-3.7 

-2.1 

Per'cenT  oT""* 
Timely 
Aica 

N/A 

_ 

98 

76 

99 

87 

87 


Table  11.  Total  Cloud  Forecast  BIAS  (in  Percent  Cloud  Cover) 
for  10  Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER 
Model,  and  RTNEPH  Persistence  for  the  00Z,  24  January  1985 
Case  Over  the  Three  Areas  of  (A)  NH,  (B)  U.S.,  and  (C)  Europe. 
(Continued) . 


jCj^Europe  BIAS 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24- HR 

36- HR 

48-HR 

Persistence 

0.0 

-11.0 

mm 

-15.7 

-12.1 

5LAYER 

0.0 

-3.6 

5.2 

-7.8 

-3.8 

GSM/A 

13.2 

17.2 

26.3 

18.0 

22.0 

GSM/N 

14.4 

14.0 

23.4 

14.2 

18.8 

GSM/E 

-25.1 

-20.3 

-10.9 

-20.5 

-14.8 

GSM/EM 

-23.4 

-18.1 

-8.2 

-17.7 

-12.3 

GSM/GL6 

5.1 

-3.3 

2.6 

-6.8 

-1.7 

GSM/GL6A 

9.0 

0.8 

6.7 

-2.6 

2.0 

GSM/GL6B 

-7.9 

-10.2 

-1.4 

-10.5 

-6.6 

GSM/GL6C 

-2.9 

-5.3 

3.2 

-6.0 

-2.5 

GSM/GL4 

-2.7 

-6.6 

1.4 

-8.0 

-4.2 

GSM/GL4A 

— 

-2.2 

5.5 

-3.6 

-0.2 

Percent  of 
Timely 
Area 

N/A 

78 

63 

78 

75 

88 


Table  12.  Total  Cloud  Forecast  20/20  Accuracy  Score  (Percent 
of  Grid  Points  Correct  to  Within  20  Percent  Cloud  Cover)  for  Ten 
Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER  Model, 
and  RTNEPH  Persistence  for  the  00Z,  17  Januaiy  1985  Case  Over 
the  Three  Areas  of  (A)  U.S.,  (B)  Europe,  and  (C)  NH  (See  text  for 
connotation  of  brackets,  parentheses,  and  underscores). 


(A)  U.S. _  20/20  Score _ OOZ,  17  Januaiy  1985 


PERIOD 

CLOUD 

SCHEME 

00-HR 

12-HR 

24-HR 

36- HR 

48-HR 

Persistence 

100.0 

56.9 

48.8 

44.5 

44.8 

mm 

[100.0] 

56.5 

52.1 

48.8 

46.7 

GSM/A 

44.9 

56.4 

53.9 

42.7 

42.1 

GSM/N 

43.4 

56.3 

55.2 

44.8 

43.5 

GSM/E 

33.1 

45.5 

53.8 

[53.7] 

54.3 

GSM/EM 

34.0 

47.2 

55.3 

52.8 

53.3 

GSM/GL6 

43.6 

[57.2] 

52*3 

46.2 

49.0 

GSM/GL6A 

45.2 

(57.0) 

[60.6] 

46.8 

48.8 

GSM/GL6B 

42.4 

49.9 

57.3 

52.4 

53.8 

GSM/GL6C 

44.8 

51.6 

58.6 

(53.5) 

54.6 

GSM/GL4 

43.8 

50.8 

59.5 

52.1 

[56.0) 

GSM/GL4A 

(45.7) 

52.7 

(60.3) 

52.1 

(55.6) 

Percent  of 
Timely 
Area 

N/A 

_ 

90 

_ 

80 

94 

62 

89 


Table  12.  Total  Cloud  Forecast  20/20  Accuracy  Score  (Percent 
of  Grid  Points  Correct  to  Within  20  Percent  Cloud  Cover)  for  Ten 
Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER  Model, 
and  RTNEPH  Persistence  for  the  00Z,  17  January  1985  Case  Over 
the  Three  Areas  of  (A)  U.S.,  (B)  Europe,  and  (C)  NH  (See  text  for 
connotation  of  brackets,  parentheses,  and  underscores).  (Continued). 


(3)  Europe_ 20/20  Score_ 00Z,  17  January  1985 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24-HR 

36- HR 

48- HR 

Persistence 

100.0 

52.5 

48.1 

39.7 

39.5 

5LAYER 

[100.0] 

[62.4] 

[54.9] 

49.5 

45.5 

GSM/A 

50.2 

(55.9) 

48.9 

50.2 

44.7 

GSM/N 

51.4 

55.7 

49.3 

50.9 

44.6 

GSM/E 

35.6 

44.0 

48.1 

45.7 

46.9 

GSM/EM 

36.1 

45.9 

49.3 

46.3 

47.3 

GSM/GL6 

54.6 

52.0 

52.8 

50.7 

50.1 

GSM/GL6A 

57,2 

53.5 

(53.5) 

[51.0] 

(50.9) 

GSM/GL6B 

49.3 

50.7 

48.8 

47.8 

50.0 

GSM/GL6C 

53.4 

52.5 

49.8 

47.4 

[51.0] 

GSM/GL4 

54.5 

51.6 

50.9 

49.8 

49.7 

GSM/GL4A 

(57.9) 

54.1 

50.6 

49.4 

50.1 

Percent  of 
Timely 
Area 

N/A 

79 

71 

85 

69 

90 


Table  12.  Total  Cloud  Forecast  20/20  Accuracy  Score  (Percent 
of  Grid  Points  Correct  to  Within  20  Percent  Cloud  Cover)  for  Ten 
Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER  Model, 
and  RTNEPH  Persistence  for  the  00Z,  17  January  1985  Case  Over 
the  Three  Areas  of  (A)  U.S..  (B)  Europe,  and  (C)  NH  (See  text  for 
connotation  of  brackets,  parentheses,  and  underscores).  (Continued). 


(C)  NH_ 20/20  Score_ 00Z,  17  January  1985 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24-HR 

36- HR 

48- HR 

Persistence 

100.0 

54.8 

52.3 

47.4 

47.7 

5LAYER 

[100.0] 

[56.6] 

(52.0) 

[50.8] 

(47.5) 

GSM/A 

41.4 

45.2 

39.0 

40.4 

38.0 

GSM/N 

42.5 

44.7 

39.8 

42.4 

37.5 

GSM/E 

38.6 

45.9 

[52.11 

5£Li 

147.91 

GSM/EM 

38.9 

46.6 

51.9 

49.7 

47.0 

GSK/GL6 

44.7 

49.1 

46.8 

46.2 

43.8 

GSM/GL6A 

45.8 

49.3 

46.8 

46.2 

43.7 

GSM/GL6B 

44.7 

47.9 

49.7 

49.9 

47.2 

GSM/GL6C 

46.4 

48.7 

49.8 

50.1 

47.0 

GSM/GL4 

47.8 

49.1 

50.6 

(50.6) 

47.4 

GSM/GL4A 

(49.0) 

(49.7) 

49.9 

50.0 

46.8 

N/A 

84 

77 

86 

69 

91 


Table  13.  Total  Cloud  Forecast  20/20  Accuracy  Score  (Percent 
of  Grid  Points  Correct  to  Within  20  Percent  Cloud  Cov~r)  for  Ten 
Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER  Model, 
and  RTNEPH  Persistence  for  the  00Z,  24  January  1985  Case  Over 
the  Three  Areas  of  (A)  U.S.,  (B)  Europe,  and  (C)  NH  (See  text  for 
connotation  of  brackets,  parentheses,  and  underscores). 


(A)  U.S. 20/20  Score_ 00Z,  24  January  1985 


PERIOD 

CLOUD 

SCHEME 

00- HR 

12-HR 

24-HR 

36-HR 

48-HR 

Persistence 

100.0 

56.4 

49.6 

40.7 

43.6 

5LAYER 

1100.0] 

[51.71 

48.3 

(42.5) 

42.7 

GSM/A 

(49  6) 

(46.4) 

46.2 

41.4 

42.3 

GSM/N 

47.0 

45.7 

46.3 

41.3 

41.6 

GSM/E 

25.7 

33.0 

42.5 

39.5 

47.9 

GSM/EM 

51.5 

34.5 

44.2 

40.7 

49.2 

GSM/GL6 

45.0 

43.3 

46.3 

41.8 

45.7 

GSM/GL6A 

46.2 

46.2 

48.2 

42.3 

44.9 

GSM/GL6B 

39.4 

38.6 

47.9 

42.4 

[50.5] 

GSM/GL6C 

41.6 

42.4 

[50.5] 

[42.6] 

49.9 

GSM/GL4 

45.2 

41.1 

47.5 

41.8 

(50.0) 

GSM/GL4A 

46.6 

42.8 

(50.0) 

42.2 

49.5 

Percent  of 

Timely 

Area 

N/A 

98 

76 

99 

87 

92 


Table  13.  Total  Cloud  Forecast  20/20  Accuracy  Score  (Percent 
of  Grid  Points  Correct  to  Within  20  Percent  Cloud  Cover)  for  Ten 
Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER  Model, 
and  RTNEPH  Persistence  for  the  00Z,  24  January  1985  Case  Over 
the  Three  Areas  of  (A)  U.S.,  (B)  Europe,  and  (C)  NH  (See  text  for 
connotation  of  brackets,  parentheses,  and  underscores).  (Continued). 


(B)  Europe _ 20/20  Score _ 00Z,  24  January  1985 


PERIOD 

CLOUD 

SCHEME 

00-HR 

12- HR 

24-HR 

36-HR 

48- HR 

Persistence 

100.0 

51.3 

53.2 

43.5 

41.3 

5LAYER 

[100.0] 

(60.1) 

HRH 

46.3 

[46.2] 

GSM/A 

51.0 

58.0 

44.3 

(53.0) 

27.5 

GSM/N 

49.3 

58.5 

46.8 

[53.6] 

(45.8) 

GSM/E 

47.1 

54.3 

50.1 

43.4 

40.4 

GSM/EM 

48.3 

48.5 

50.4 

44.8 

40.6 

GSM/GL6 

52.7 

[60.2] 

[52.9] 

50.7 

44.8 

GSM/GL6A 

53.3 

59.9 

51.5 

51.6 

45.2 

GSM/GL6B 

56.0 

57.0 

51.3 

48.4 

42.7 

GSM/GL6C 

56.2 

58.2 

50.3 

49.6 

42.9 

GSM/GL4 

(57.2) 

58.4 

(52.4) 

49.8 

43.8 

GSM/GL4A 

56.6 

58.8 

51.4 

51.5 

44.5 

Percent  of 
Timely 
Area 

N/A 

78 

63 

78 

75 

93 


Table  13.  Total  Cloud  Forecas:  ^0/20  Accuracy  Score  (Percent 
of  Grid  Points  Correct  to  Within  20  Percent  Cloud  Cover)  for  Ten 
Global  Spectral  Model  (GSM)  Cloud  Schemes,  the  5LAYER  Model, 
and  RTNEPH  Persistence  for  the  00Z,  24  January  1985  Case  Over 
the  Three  Areas  of  (A)  U.S.,  (B)  Europe,  and  (C)  NH  (See  text  for 
connotation  of  brackets,  parentheses,  and  underscores).  (Continued). 


(C)  NH_ 20/20  Score_ 00Z,  24  January  1985 


PERIOD 

CLOUD 

SCHEME 

00-HR 

12-HR 

24-HR 

36- HR 

48-HR 

Persistence 

100.0 

56.4 

54.5 

46.3 

48.0 

5LAYER 

(100.01 

[58.1] 

[53.4] 

[46.31 

44.7 

GSM/A 

41.3 

41.1 

42.4 

39.8 

40.9 

GSM/N 

41.1 

41.4 

42.8 

40.3 

40.3 

GSM/E 

42.5 

45.2 

47.7 

44.4 

45.7 

GSM/EM 

42.9 

45.6 

48.1 

44.1 

45.5 

GSM/GL6 

44.4 

47.7 

48,4 

42.3 

43.8 

GSM/GL6A 

45.4 

47.7 

48.0 

42.2 

43.7 

GSM/GL6B 

46.1 

48.2 

49.4 

44.6 

(46.5) 

GSM/GL6C 

47.3 

48.8 

(49.9) 

44,8 

[46.51 

GSM/GL4 

48.7 

48.6 

48.0 

45.1 

46.4 

GSM/GL4A 

(49.1) 

(48.9) 

48.1 

(45.2) 

46.2 

Percent  of 

Timely 

Area 

N/A 

88 

62 

85 

73 
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The  cloud  schemes  designate^  GL6A,  GL6B,  GL6C,  and  GL4A  In  Tables  9-13  represent 
sensitivity  tests  of  variations  in  the  baseline  vertical  stacking  algorithm  (Appendix  A)  used  in 
GL6  and  GL4.  In  the  scheme  j  GL4A.  GL6A,  and  GL6C,  random  vertical  stacking  was  used 
(l.e.  R=1  in  layer  Combinations  I,  II,  and  III  of  Appendix  A),  while  scheme  GL6B  used  the 
standard  values  of  vertical  stacking  (that  is,  0  <  R  <  1  as  given  for  layer  Combinations  I,  II, 
and  III  in  Appendix  A).  In  the  schemes  GL6B  and  GL6C,  the  GSM  layer  cloud  amounts  at  100 
kPa  were  set  to  zero.  The  motivation  for  these  stacking  experiments  will  be  presented  in 
Section  3.2.3. 

Although  the  GSM  forecasts  were  executed  out  to  96  hours,  (for  example,  Figures  1 1  and  12), 
the  verifications  in  Tables  10-13  were  carried  out  to  only  48  hours.  Tne  latter  is  the  forecast 
length  of  the  operational  5LAYER  model  in  the  NH,  and  therefore  AFGWC  archived  KTNEPH 
verifying  cloud  analyses  only  out  to  48  hours  for  this  study. 

To  better  show  overall  systematic  trends,  the  calculation  of  mean  total  cloud  amounts  in 
Table  9  utilized  every  1/2-mesh  grid  point  in  the  NH  octagon.  Tables  10-13  on  the  other  hand, 
focus  on  potntwisc  cloud  forecast  accuracy,  particularly  Tables  12-13.  Hence  the  computations 
yielding  the  statistics  of  Tables  10-13  ignored  grid  points  where  the  verifying  RTNEPH  cloud 
analysis  was  more  than  three  hours  old  (for  example,  in  surface  data  sparse  regions  between 
successive  orbiting  satellite  passes).  For  this  purpose,  the  RTNEPH  database  includes 
pointwise  "time  flags"  indicating  the  timeliness  of  the  analyzed  cloud  amount  at  each  grid 
point. 

Thus  the  sample  of  points  verified  within  a  given  area  differed  somewhat  over  the  various 
forecast  times.  However,  for  a  given  forecast  time,  all  the  cloud  forecast  schemes  in  the  tables 
were  verified  at  exactly  the  same  sample  of  timely  points.  The  percent  of  timely  area  is  given 
in  the  Tables  10-13  to  explain  some  of  the  temporal  idiosyncrasies  in  the  statistics. 

We  consider  first  the  results  in  Table  9.  Over  an  area  as  large  as  the  NH  octagon,  a  good 
cloud  forecast  scheme  should  exhibit  a  mean  cloud  amount  that  (1)  is  nearly  constant  over  the 
entire  forecast  period  (adjusts  to  the  model  s  spin-up  behavior  in  RH)  and  (2)  agrees  closely 
with  the  verifying  mean  cloud  amount  (negligible  bias).  The  first  criterion  is  justified  by  the 
RTNEPH  mean  cloud  amounts  at  the  top  of  Table  9,  which  show  a  quasi-constant  mean  cloud 
cover  of  about  47  percent,  with  12-hourly  deviations  within  about  plus  or  minus  4  percent. 

The  deviations  in  the  bias  of  RTNEPH  persistence  over  the  NH  octagon  in  Tables  10a  and  11a, 
which  represent  two  cases,  also  reflect  small  temporal  changes  in  observed  total  cloud  cover 
over  a  hemispheric  domain.  We  shall  view  temporal  oeviations  in  forecast  mean  cloud  cover 
over  the  octagon  as  acceptable  if  bounded  by  this  observed  variability  of  about  5  percent.  The 
bracketed  bias  values  in  Tables  10A  and  11A  denote  those  instances  in  which  forecast  bias 
over  the  octagon  was  5  percent  or  less.  (The  5  percent  bias  threshold  was  not  highlighted  for 
the  U.S.  and  European  regions;  because  for  a  single  forecast  case  over  a  limited  area,  diurnal 
variations  or  translating  weather  systems  can  affect  the  calculated  bias  irrespective  cf  the 
model's  systematic  trend  —  hence  only  substantially  large  biases  are  meaningful  in  the 
regional  bias  statistics  in  Tables  10-11.) 

Turning  to  Tables  12-13,  we  have  highlighted  certain  values  of  the  20/20  score  in  these 
tables.  Scores  in  brackets  and  parr  a theses  denote  the  first  and  second  highest  forecast  model 
scores,  respectively,  for  each  forecast  time.  Underlined  scores  represent  "honorable  mentions" 
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to  help  identify  cloud  schemes  that  are  consistently  providing  competitive  scores,  albeit  not 
always  the  best  scores.  Bold  persistence  scores  indicate  when  persistence  beat  all  model 
forecasts.  Also  of  note  in  Tables  12  and  13,  we  purposely  first  present  and  henceforth 
emphasize  the  U.S.  and  European  regional  20/20  scores  instead  of  the  NH  octagon  scores. 
Though  the  scores  for  the  GSM  cloud  schemes  are  admittedly  less  competitive  over  the  octagon, 
here  we  Justify  the  emphasis  on  regional  comparisons  for  the  following  reasons: 


(1)  the  greater  availability  of  surface  reports  to  augment  satellite  data,  thereby  increasing 
the  accuracy  of  the  verifying  RTNEPH, 


(2)  the  greater  availability  of  RAOBs  to  correct  the  severe  upper-level  dry  bias  in  the  first- 
guess  forecast  used  by  the  AFGWC  moisture  analysis  over  the  octagon  (Sec.  2.1), 


(3)  the  greater  relevance  of  the  chosen  regions  to  the  preponderance  of  Air  Force 
operations. 


The  results  of  Tables  9-13  and  example  displays  of  cloud  forecast  fields  will  be  assessed  in 
detail  in  the  following  sections  according  to  three  groupings  of  forecasts:  5LAYER  and 
RTNEPH  persistence  in  Section  3.2.1,  GSM  using  previous  diagnostic  cloud  schemes  in 
Section  3.2.2,  and  GSM  using  the  GL  cloud  schemes  in  Section  3.2.3.  As  a  starting  point,  we 
list  below  the  major  conclusions  to  be  drawn  from  the  tables  and  later  cloud  displays: 


(1)  The  margin  of  skill  over  persistence  in  both  the  5LAYER  and  GSM  cloud  schemes  is 
disappointingly  small.  The  accuracy  of  persistence  is  rather  hard  to  surpass  at  24 
hours  or  less.  Even  at  48  hours,  the  margin  of  GSM  and  5LAYER  forecast  skill  over 
persistence  is  lackluster.  The  small  margin  of  5LAYER  and  GSM  skill  over  persistence 
was  also  seen  in  the  nine  winter  forecast  cases  in  the  forerunner  study  by  MW6. 


(2)  The  surprisingly  low  absolute  accuracy  of  5LAYER  and  GSM  cloud  forecasts  at  24-48 
hours,  with  20/20  scores  only  in  the  range  of  40-55  percent,  despite  the  fact  that  the 
definition  of  the  score  allows  a  20  percent  cloud  cover  forecast  error  in  a  "hit"  is 
another  disappointment.  Even  at  the  short  12-hour  range,  accuracy  scores  much  above 
60  percent  are  rare.  These  low  scores  reflect  the  quasi-binary  nature  of  the  verifying 
cloud  fields  (either  mostly  clear  or  mostly  cloudy),  which  increases  the  chances  of  large 
forecast  busts. 
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(3)  The  5LAYER  forecasts: 


(a)  show  an  acceptably  small  but  steady  decrease  in  NH  mean  cloud  amount  over  time 
and  hence  small  but  generally  increasing  negative  biases:  (b)  show  a  consistent 
monotonic  decrease  in  20/20  score  or  accuracy  with  forecast  time,  but  nevertheless 
exhibit  the  highest  12-hour  accuracy  scores  and  often  the  highest  24-hour  accuracy 
scores:  (c)  compared  to  the  GSM/GL4  scheme  in  the  U.S.  and  European  areas,  exhibit 
36-  and  48-hour  accuracy  scores  that  are  always  lower  in  Case  I  and  usually  lower  in 
Case  II. 


(4)  The  GSM  forecasts  using  previous  cloud  schemes: 


(a)  show  a  very  large  spin-up  (increase)  in  cloud  amounts  in  the  first  24  hours;  (b)  show 
unacceptably  large  positive  or  negative  biases;  (c)  do  not  surpass  either  5LAYER  or 
GSM/GL  schemes  in  accuracy  with  any  consistency  whatsoever. 


(5)  The  GSM  forecasts  using  the  GL  cloud  schemes: 


(a)  in  the  case  of  the  GL4  scheme,  show  strikingly  negligible  spin-up  behavior,  nearly 
steady  NH  mean  cloud  amounts,  and  small  NH  octagon  biases  (this  despite  the  huge 
spin-up  noted  in  Section  2.2  in  the  GSM  upper-level  RH  forecasts);  (b)  in  the  case  of 
the  GL6  scheme,  also  show  small  spin-up  behavior,  but  an  unexpected  nontrivial 
bias  (positive)  fThe  latter  result  prompted  the  cloud  stacking  experiments  embodied 
in  the  GL6A,  GL6B,  and  GL6C  schemes);  (c)  consistently  provide  more  accuracy  at  all 
times  than  the  GSM  forecasts  using  previous  cloud  schemes,  and  often  surpass 
5LAYER  accuracy  at  24  hours  and  beyond  over  the  U.S.  and  European  areas; 

(d)  show  such  disappointingly  low  accuracy  at  the  initial  (0-hour)  time  that  accuracy 
actually  increases  during  the  first  12,  24,  or  even  36  hours;  (e)  show  rather  impressive 
accuracy  scores  at  48  hours  relative  to  5LAYER  in  several  instances. 


In  the  following  sections,  we  examine  displays  of  the  5LAYER  and  GSM  layer  and  total 
cloud  fields  over  the  NH  for  Case  I.  These  forecasts  will  be  compared  to  each  other  and  to 
displays  of  the  verifying  RTNEPH  cloud  analyses.  As  our  reference  point,  we  begin  by 
examining  the  latter  analyses  in  Figure  30  (total  cloud)  and  Figures  31  and  32  (layer  cloud), 
which  show  the  RTNEPH  in  Case  I  at  basetime  and  48  hours  later.  fThe  shading  intensity  key 
of  Figure  1  applies  to  these  and  all  later  cloud  displays.)  To  provide  a  further  reference  and  an 
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aid  In  synoptic  Interpretation,  we  provide  In  Figure  29  the  ECMWF  operational  12Z  analyses32 
of  the  50  kPa  height  field  and  sea-level  pressure  field  at  the  12-  and  60-hour  points  of  Case  I. 


b 

Figure  29.  The  ECMWF  Opera tional  NH  12Z  Analyses32  of  the 
50  kPa  Height  Field  (top)  and  S^a-level  Pressure  Field  (bottom) 
for  (a)  17  Januaiy  and  (b)  19  Januaiy  1985, 


Figure  31.  The  RTNEPH  OOZ  Analysis  of  NH  Octagon  Layer  Cloud 
Amount  for  17  January  1985  for  the  (a)  Gradient,  (b)  85,  (c)  70  and 
(d)  30  kPa  Layers.  fThe  shading  key  In  Figure  1  applies  to  all  shaded 
displays  of  cloud  amount). 


Figure  32.  The  RTNEPH  OOZ  Analysis  of  NH  Octagon  Laver  Cloud 
Amount  for  19  January  1985  for  the  (a)  Gradient,  (b)  85,  (c)  70  and 
(d)  30  kPa  Layers.  fThe  shading  key  in  Figure  1  applies  to  all  shaded 
displays  of  cloud  amount). 
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The  distinctive  features  of  the  total  cloud  analyses  in  Figure  30  are  1)  the  distinctive 
frontal  bands  in  certain  areas  (such  as,  the  eastern  Pacific),  2)  the  persistent  cloud-free 
character  of  the  trade  wind  and  desert  areas  of  the  subtropical  Atlantic  and  Pacific,  North 
Africa,  and  the  southwest  U.S.,  and  3)  the  preponderance  of  cloud  over  the  extratropical 
oceans. 

The  foremost  feature  of  the  layer  cloud  analyses  of  Figures  31  and  32  is  the  significantly 
greater  extent  of  cloud  cover  in  the  lower  layers.  (The  RTNEPH  analyses  of  cloud  cover  at  50 
kPa,  not  shown,  appear  similar  to  the  given  30  kPa  analyses).  This  preponderance  of  low- 
level  versus  high-level  cloud  is  observationally  realistic  and  is  a  characteristic  common  to 
other  cloud  databases,  as  described  in  the  study  of  Henderson-Sellers.31  Finally,  Figures 
3 1  -32  show  that  the  larger  total  cloud  amounts  over  the  oceans  noted  in  Figure  30  stem 
largely  from  clouds  in  the  lower  layers.  The  overall  dominance  of  low-level  clouds  in  these 
RTNEPH  analyses  demonstrate  that  for  a  consistently  good  total  cloud  forecast,  it  is  virtually 
essential  to  have  a  good  low-level  cloud  forecast. 

Following  on  this  latter  point,  it  is  crucial  to  note  that  the  vertical  compaction  (Figure  7) 
of  the  originally  retrieved  RTNEPH  cloud  layers  to  the  five  layers  of  the  5LAYER  model 
includes,  as  a  check,  an  application  of  the  cloud  stacking  algorithm  in  Appendix  A.  This 
check  is  followed  if  necessary  by  an  iterative  adjustment  of  the  compacted  layer  cloud 
amounts  until  the  cloud  stacking  gives  a  total  cloud  amount  that  agrees  within  a  few  percent  of 
the  independently  obtained  RTNEPH  1 /2-mesh  total  cloud  amount. 

The  aforementioned  study  of  Henderson-Sellers31  and  the  companion  study  of  Hughes  and 
Henderson-Sellers30  examined  the  AFGWC  cloud  analyses  of  January  and  July  1979  from  the 
3DNEPH  model25  (the  direct  RTNEPH  predecessor  which  used  very  similar  cloud  retrieval 
algorithms).  These  studies  concluded  that  the  3DNEPH  provides  reliable  cloud  distributions 
agreeing  with  known  features  of  the  general  circulation  and  other  cloud  databases,  except  for 
1)  some  misrepresentation  of  cloud  cover  over  the  immediate  polar  regions  (not  included  in  the 
U.S.  and  European  areas  here),  2)  somewhat  underestimation  of  high  clouds  in  general  and 
stratus  clouds  at  low  levels  over  oceans,  and  3)  overestimation  of  clouds  associated  with  the 
NH  summer  monsoon  (not  present  in  the  NH  winter  cases  considered  here).  Despite  the  above 
modest  weaknesses,  daily  operational  experience  with  the  RTNEPH  and  its  regular  comparison 
with  satellite  cloud  imagery  at  AFGWC  show  it  to  be  vastly  more  accurate  than  the  cloud  model 
forecasts  and  hence  a  fundamental  tool  in  the  assessment  of  cloud  model  forecast  accuracy. 

3.2.1  PARALLEL  SLAYER  MODEL  FORECASTS 

The  cloud  forecast  fields  of  the  AFGWC  operational  5LAYER  model  valid  at  48-hours  in 
Case  I  are  given  in  Figure  33  (total  cloud)  and  Figure  34  (layer  cloud).  In  AFGWC  parlance,  and 
strictly  speaking,  the  5LAYER  forecasts  in  Figures  33  and  34  are  45-hour  forecasts,  because 
they  represent  the  results  of  a  quasi-Lagrangian  temporal  integration  over  15  3-hour  time 
steps  starting  from  a  nominal  RTNEPH  basetime  of  03Z  rather  than  00Z  on  17  January  1985. 
However,  as  Justified  below,  these  same  5LAYER  forecasts  can  be  viewed  also  as  48-hour 
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forecasts  starting  from  a  00Z  basetime,  which  Is  the  view  we  chose  to  adopt  for  convenience  In 
labeling  the  5LAYER  statistics  In  Tables  9-12.  An  understanding  of  this  ambiguity  reveals  one 
of  the  two  dominant  factors  leading  to  the  clear  superiority  of  5LAYER  over  the  GSM  cloud 
schemes  In  20/20  scores  over  the  0-12  or  0-24  hour  periods  In  Tables  12-13. 


► 


Figure  33.  The  5LAYER  45-hour  Forecast  of  NH  Octagon  Total 
Cloud  Amount  Valid  at  00Z,  19  January  1985. 
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b 


Figure  34.  The  5LAYER  45-hour  Forecast  of  NH  Octagon  Laver 
Cloud  Amount  valid  at  00Z,  19  January  1985,  for  the  (a)  Gradient, 
fb)  85,  (c)  70  and  (d)  30  kPa  layers. 
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To  sort  out  the  ambiguity,  recall  from  Section  1  that  the  5LAYER  model  calculates  its 
quasi-Lagrangian  air  parcel  trajectories  from  the  3-D  wind  velocity  forecasts  of  a  previous 
execution  of  the  AFGWC  GSM.  The  AFGWC  48-hour  GSM  forecast,  which  is  executed  along 
with  its  predecessor  H1RAS  analysis  on  a  6-hourly  cycle  or  four  times  dally,  finishes  at  about 
5  hours  after  basetime  as  depicted  in  Figures  4-7  of  Stobie,2  that  is,  5  hours  after  0,  6,  12,  and 
18Z.  In  contrast,  the  quasi-Lagrangian  efficiency  and  quick  execution  of  the  5LAYER  model 
allows  it  to  start  by  2  hours  after  a  given  basetime  and  finish  1  /2  hour  later.  Owing  to  the 
latter  execution  efficiency,  the  RTNEPH/5LAYER  pair  execute  on  a  3-hourly  cycle,  that 
is,  finishing  2.5  hours  after  0,  3.  6.  9,  12,  15,  18,  and  21Z.  Therefore,  a  new  03Z  RTNEPH 
analysis  is  available  and  input  to  the  first  5LAYER  execution  able  to  make  use  of  the  GSM 
wind  forecast  from  00Z.  The  official  basetime  given  to  this  5LAYER  execution  is  03Z  --  the 
nominal  time  of  the  input  RTNEPH,  which  is  used  in  the  manner  of  Figure  7  to  directly 
initialize  the  5LAYER  initial  CPS  moisture  field.  However,  the  initial  temperature  field  for 
the  03Z  5LAYER  run  is  a  persistence  of  the  00Z  HIRAS  temperature  analysis  used  by  the  00Z 
GSM.  Finally,  the  octagon-wide  time  average  of  the  RTNEPH  point-wise  time  flags  would 
typically  yield  a  mean  valid  time  of  say  1-2  hours  before  03Z.  In  summaiy  then,  the  45-hour 
forecast  of  the  03Z  5LAYER  has  used  1)  a  48-hour  GSM  wind  forecast  from  00Z,  2)  an  initial 
temperature  analysis  of  00Z,  and  3)  an  RTNEPH  analysis  valid  in  an  area-mean  sense 
somewhat  before  03Z. 

The  key  revelation  of  the  above  discussion  is  that  owing  to  the  quick  execution  time  of  the 
RTNEPH/5LAYER  cycle  (three  hourly)  relative  to  the  HIRAS/GSM  cycle  (6-hourly),  the  5LAYER 
run  most  closely  parallel  to  a  given  GSM  run  is  able  to  begin  with  initial  moisture  fields 
derived  from  an  RTNEPH  6-hours  more  current  than  the  RTNEPH,  influencing  the  initial 
AFGWC  moisture  analyses  used  to  initialize  the  GSM  runs  in  this  study.  To  see  this,  we  recall 
from  Section  2. 1  that  the  first-guess  moisture  forecast  for  the  00Z  AFGWC  moisture  analyses 
used  in  the  00Z  GSM  runs  here  was  the  3-hour  5 LAYER  forecast  from  the  21Z  RTNEPH 
analysis.  Hence  the  one  dominant  factor  in  the  superiority  of  5LAYER  short-range  cloud 
forecasts  versus  the  GSM  is  that  its  superior  execution  efficiency  allows  it  to  utilize  a  more 
current  RTNEPH  analysis.  (Even  if  the  00Z  RTNEPH  were  used  directly  as  a  pseudo  moisture 
observation  source  in  the  00Z  HIRAS  moisture  analysis  of  today,  the  06Z  RTNEPH  for  the  06Z 
5 LAYER  run  would  still  be  6  hours  more  current  than  the  00Z  RTNEPH  input  to  the  00Z 
HIRAS.  It  is  not  until  1 1Z  that  the  next  48-hour  GSM  forecast  would  be  available  from  the  06Z 
GSM  run. 

Furthermore,  as  illustrated  later  in  Section  3.2.2,  an  equally  dominant  factor  in  5LAYER 
short-range  superiority  in  this  GSM  comparison  is  the  extremely  direct  and  fairly  high 
resolution  processing  route  illustrated  in  Figure  7  from  a  1 /2-mesh  RTNEPH  layer  cloud 
analysis  to  the  initial  5LAYER  layer  CPS  humidity  field,  and  the  counterpart  short  and  direct 
inverse  route  from  5 LAYER  CPS  forecasts  to  1 /2-mesh  layer  clouds.  In  the  GSM  case  (recall 
Section  3.1  and  see  later  in  Section  3.2.2),  this  preprocessing/ postprocessing  circuit  involves 
significantly  moie  transformation  steps,  some  occurring  at  resolutions  worse  than  1 /2-mesh. 
In  ..summary  then  .the_paramount  strengths  of  the  5  LAYER  model  for  short-range  cloud 
forecasting  IQ- 12  hours)  are  its  direct  and,  timely  use  of  the  RTNEPH  cloud  analysis  in  deriving 
its  initial  moisture  field. 
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Surprisingly,  this  short-range  advantage  of  5LAYER  Is  realized  despite  the  pervasive  upper- 
troposphere  dry  bias  In  the  initial  5LAYER  CPS  fields  --  a  bias  clearly  demonstrated  In 
Section  2.1.2  and  arising  from  1)  the  somewhat  underestimated  high  cloud  amounts  In  the 
RTNEPH  and  2)  the  very  low  critical  humidity  values  (RHc)  implicit  In  the  CPS  curves 
(Figure  21).  The  answer  to  this  crucial  contradiction  is  that  5 LAYER  sustains  its  Initial  drv 

*  bias  throughout  the  forecast.  Thus,  while  any  user  retrieving  5LAYER's  upper-level  CPS 
hurrUrHfv  forecasts  would  find  them  very  poor  (much  too  dry)  compared  to  RAOBs,  the  5 LAYER 
layer  cloud  forec.  and  hence  total  cloud  forecasts  remain  satisfactory,  because  5LAYER  does 

*  not  significantly  ^jlsten  or  dry  any  of  its  layers.  Hence,  because  the  same  CPS  curves  with 
the  same  critical  humidity  values  are  applied  in  the  Inverse  sense  at  forecast  time  to  get  layer 
clouds  from  forecast  layer  CPS,  the  resulting  cloud  amounts  remain  reasonable  showing  little 
bias  with  respect  to  the  RTNEPH  layer  cloud  amounts. 

These  latter  observations  are  substantiated  by  the  5LAYER  layer  cloud  forecast  displays  of 
Figure  34d,  which  show  that  while  5LAYER  yields  a  more  coherent  frontal  or  banded  structure 
than  the  RTNEPH  analyses  at  the  various  levels.  It  does  not  show  a  significantly  increased  or 
decreased  cloud  amount  over  the  domain  as  a  whole.  Similarly,  the  5LAYER  mean  total  cloud 
amounts  In  Table  9  and  total  cloud  biases  in  Tables  10  and  1 1  show  only  a  small  decrease  in 
total  cloud  cover  with  forecast  time.  Hence  5LAYER  does  not  exhibit  a  large  spin-up  behavior, 
despite  Its  dry  upper-level  Initial  state. 

Although  the  decrease  in  5LAYER  total  cloud  cover  over  48-hours  Is  only  6.4  percent  In 
Table  9.  there  are  well-known  physical  reasons  for  this  small  drying  tendency  in  5LAYER. 
Specifically,  the  5LAYER  model  includes  precipitation  processes  in  saturated  ascending  parcels 
as  a  moisture  sink,  but  It  Ignores  surface  evaporation  as  a  counterpart  moisture  source.4  At 
the  precipltatlon-prone  lower  levels  of  85  and  70  kPa,  a  comparison  of  the  5LAYER  forecasts 
In  Figures  34b  and  34c  with  the  RTNEPH  analyses  at  those  levels  in  Figures  32b  and  32c 
suggests  a  small  decrease  in  cloud  amount  at  these  levels  during  the  5LAYER  forecast. 

Overall,  there  Is  an  Impressive  similarity  in  the  character  of  the  5LAYER  forecast  total 
cloud  field  in  Figure  33  and  the  verifying  RTNEPH  total  cloud  In  Figure  30b.  The  sharp 
gradients  In  cloud  amounts,  the  relative  size  of  large  synoptic  patterns,  the  ragged  appearance 
of  cloud  shield  edges,  and  the  location  of  persistently  clear  areas  In  the  RTNEPH  analysis  are 
reproduced  well  by  SLAYER.  Stated  alternatively,  the  subjective  impression  Imparted  by  the 
5LAYER  cloud  displays  is  that  the  fractal  character  or  geometry  of  .  he  5LAYER  total  cloud 

'  forecast  agrees  well  with  that  of  the  RTNEPH  analysis. 

As  a  final  point,  a  physical  overlay  of  Figure  31b  onto  Figure  31a  (say  via  transparencies) 
highlights  the  persistent  nature  of  severed  large  cloud-free  and  cloudy  regions.  This  illustrates 

*  the  difficult  challenge  of  surpassing  persistence  skill  In  Tables  12-13. 

3.2.2  PREVIOUS  DIAGNOSTIC  CLOUD  SCHEMES 

Tables  9-13  showed  that  cloud  forecasts  from  the  GSM  using  the  previous  diagnostic  cloud 
schemes  suffered  from  a  large  spin-up  behavior  of  increasing  cloud  amounts,  large  biases,  and 
generally  lower  accuracy  scores.  We  will  investigate  example  displays  of  these  GSM  total  and 
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layer  cloud  forecasts  In  Case  I  to  substantiate  this  behavior  and  illustrate  important 
characteristics  t herein. 

1  he  fundamental  basis  lor  the  above  sp.  *p  and  bias  behavior  goes  back  to  the  central 
features  of  the  GSM  relative  humidity  forecasts  documented  in  Section  2.2,  which  the  reader 
is  urged  to  review.  The  central  feature  was  the  rapid  spin-up  of  the  GSM  RH  forecasts  by  48 
hours  to  a  distinct  model  preferred  mean  state  --  a  mean  state  by  and  large  independent  of  the  , 

Initial  mean  state  of  the  input  humidity  analysis.  This  was  clearly  illustrated  in  Figures  11- 
14. 

As  a  close  corollary,  the  magnitude  or  degree  of  the  f  pin-up  will  depend  on  how  closely  the  * 

mean  state  of  the  initial  humidity  analysis  agrees  with  the  model-preferred  state.  The  latter 
also  was  strikingly  demonstrated  in  Figures  11-14.  For  example,  in  Case  I  and  II,  the 
magnitude  of  GSM  spin- up  at  the  30  and  40  kPa  levels  is  notably  large  and  reflects  substantial 
moistening,  because  the  AFGWC  humidity  analyses  at  these  levels  are  much  drier  than  the 
GSM-preferred  humidity  state  (and  drier  than  observations  as  well).  In  contrast,  at  100  kPa, 
where  the  AFGWC  RH  analysis  over  the  octagon  agreed  fairly  well  in  a  zonal  mean  sense  with 
observations  and  two  independent  analyses  (Figures  4,  5,  and  12),  the  GSM  100  kPa  RH 
forecast  develops  a  large  positive  forecast  bias  over  oceans,  due  to  a  systematic  error  in  the 
parameterized  physics  of  the  surface  fluxes  and  vertical  diffusion  of  water  vapor  in  the  GSM 
(Section  2.2,  Figures  17a  and  18a). 

Then  as  a  second  corollary,  it  follows  that  the  model-preferred  forecast  mean  humidity 
state  mav  (1)  not  agree  with  the  obseived  mean  state  and  (2)  not  agree  with  the  preferred  mean 
state  of  a  separate  forecast  model.  Consequently,  the  previous  diagnostic  cloud  schemes  of 
GSM/A  (AFGWC:  CPS  curves  of  Figure  8  from  the  5 LAYER  model)  and  GSM/E  (the  ECMWF 
curves  of  Figure  20  from  the  ECMWF  model  prior  to  May  1985),  which  were  both  developed  for 
other  models,  or  GSM/N  (curves  of  Figure  23  used  In  the  moisture  analyses  of  first  NMC  and 
second  ECMWF)  which  was  developed  for  use  with  surface  cloud  observations  in  an  objective 
analysis  model,  have  strikingly  different  critical  values  of  RH  for  the  onset  of  various  cloud 
amount  thresholds  (Figure  21). 

Because  the  latter  critical  values  were  not  tuned  for  the  particular  model-preferred 
humidity  state  of  the  GL  baseline  GSM  here,  their  use  with  this  GSM  yields  the  large  cloud 
amount  biases  given  for  these  schemes  In  Tables  9-11.  Lastly,  these  previous  cloud  schemes 
are  fixed  schemes,  and  thus  they  cannot  adjust  in  forecast  time  for  the  large  Case  1  and  II  spin- 
up  in  (1)  mean  RH  amount  (Figure  10  and  bottom  right  of  Figure  11)  and  (2)  RH  frequency 
distribution  (Figure  19).  Thus  these  previous  schemes  yield  large  temporal  changes  in  forecast 
values  of  mean  cloud  amount  (Table  9). 

The  cloud  forecasts  in  Case  I  from  the  GSM/ A  scheme  (Figure  8)  are  displayed  for  the  0- 
and  48-hour  times  in  Figure  35  (total  cloud)  and  Figures  36-37  (layer  cloud).  As  shown  by  these 
displays,  one  advantage  of  the  GSM  over  5LAYER  is  the  GSM's  ability  to  provide  fully 
hemispheric  (or  global)  cloud  forecasts,  unlike  *he  octagon  domain  limitation  of  5LAYER  in 
Figures  33  and  34.  However,  we  limit  our  assessment  of  GSM  cloud  forecast  quality  to  the 
octagon,  as  that  is  the  domain  limit  of  the  verifying  compacted  RTNEPH  analyses  in 
Figures  30-32.  In  Figures  36  and  37,  and  all  later  GSM  layer  cloud  forecast  displays,  for 
brevity  we  show  only  tour  of  the  six  output  mandatory  pressure  levels,  as  the  50  and  40  kPa 
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cloud  fields  appear  very  similar  to  the  70  and  30  kPa  fields,  respectively,  especially  at  48 
hours  (for  example,  compare  the  70  and  50  k Pa  levels  and  the  40  and  30  kPa  levels  in  the  GSM 
RH  displays  of  Figure  18). 
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Figure  36.  For  the  AFGWC  Cloud  Scheme,  the  GSM  O-hour 
Forecast  of  NH  Laver  Cloud  Amount  Valid  at  00Z  on  17  January 
1985  for  Layers  Nominally  Centered  at  the  (a)  100,  (b)  85,  (c)  70 
and  (d)  30  kPa  Pressure  Levels. 
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The  GSM/A  scheme  is  characterized  by  a  dramatic  and  erroneous  increase  in  total  cloud 
amount,  vividly  portrayed  in  Figure  35  and  Tables  9,  10A,  and  1 1A.  Figures  36  and  37  show 
that  the  greatest  layer  cloud  increase  occurs  at  the  upper  levels,  since  the  critical  RH  values  for 
these  levels  in  the  GSM/A  (that  is,  CPS)  scheme  are  simply  much  too  low  relative  to  the  GSM- 
preferred  upper- level  RH  mean  state  (Figure  12).  The  85  and  70  kPa  levels,  on  the  other  hand, 
seem  to  show  some  decrease  in  cloud  amounts,  in  conjunction  with  the  modest  GSM  drying  at 
these  levels  apparent  from  Case  I  in  Figure  12.  On  a  positive  note,  the  85  and  70  kPa  levels  in 
Figure  36  also  show  vivid  formation  of  distinct  frontal  bands  in  response  to  the  spin-up  of  the 
divergence  and  vertical  velocity  fields,  yielding  a  more  binary  RH  forecast  distribution,  as 
shown  in  Figure  19.  Finally,  there  is  a  large  increase  in  cloud  amount  at  the  100  kPa  level, 
where  the  low-level  moist  bias  of  the  GSM  forecast  over  oceans  (notably  subtropical  and 
tropical  oceans)  manifests  itself. 

The  study  of  MW6  also  used  the  GSM/A  scheme,  though  arbitrarily  imposing  zero  cloud  on 
the  GSM  forecast  levels  of  100,  40,  and  30  kPa  to  circumvent  the  above  biases.  In  so  doing, 

MW  obtained  total  cloud  forecast  accuracy  scores  rivaling  those  of  the  GSM/GL  schemes  in 
Tables  12-13.  Such  an  approach,  though  enlightening,  is  too  ad  hoc  for  operational  use,  as  it 
fails  in  particular  those  users  who  need  forecasts  of  high  cloud. 

It  is  central  to  this  study  to  further  consider  the  0-hour  layer  cloud  fields  of  the  GSM/A 
scheme  in  Figure  36,  in  conjunction  with  the  0-hour  20/20  scores  for  this  scheme  Tables  12- 
13.  It  is  puzzling  at  first  that  (1)  the  0-hour  20/20  scores  in  Case  I  and  II  are  so  low  for  the 
GSM/A  scheme  and  (2)  the  0-hour  GSM/A  cloud  fields  in  Figure  36  do  not  better  resemble  the 
RTNEPH  layer  cloud  fields  of  Figure  31,  especially  in  RAOB  sparse  areas.  After  all,  the 
hallmark  of  the  AFGWC  initial  moisture  analyses  used  here  to  initialize  the  GSM  moisture  in 
Cases  I  and  n  is  the  use  of  a  3-hour  5LAYER  CPS  forecast,  from  the  21Z  RTNEPH,  as  the  first- 
guess  for  the  00Z  humidity  analysis.  In  fact,  as  the  reader  will  recall  from  Section  2.1.1,  this 
study  effort  went  to  great  lengths  to  use  the  particular  AFGWC  moisture  analysis  that  used  the 
RTNEPH/5LAYER  first-guess,  in  order  to  reduce  the  short-range  forecast  advantage  of  5LAYER 
stemming  from  5LAYER's  RTNEPH-based  moisture  initialization  method  of  Figure  7. 

The  3-hour  5LAYER  CPS  first-guess  forecast  valid  at  00Z  will  be  a  virtual  persistence  of  the 
21Z  initial  5LAYER  CPS  field,  obtained  directly  from  the  21Z  RTNEPH,  also  according  to 
Figure  7.  Most  notably,  here  in  the  GSM/A  scheme  we  are  using  the  same  CPS  curves  to  invert 
the  initial  GSM  humidity  fields  back  to  a  "synthesized"  0-hour  cloud  analysis.  (We  use  the 
term  "synthesized"  here  with  the  GSM  0-hour  "cloud"  fields  to  mean  the  GSM  inferred  or 
derived  cloud  analysis  obtained  by  all  the  static  preprocessing  and  postprocessing  steps 
involved  in  the  humidity/cloud  forecast  calculations  of  the  GSM,  short  of  actually  taking  a 
forward  time  integration  step  in  the  model.  Hence  the  label  "SYN"  shown  with  all  GSM  0-hour 
cloud  "analyses"  in  Figures  35-47).  Therefore,  Insofar  as  we  can  view  the  3-hour  5LAYER  CPS 
forecast  from  21Z  as  a  very  close  surrogate  (similar  to  persistence)  of  the  CPS  fields  obtained 
directly  from  the  00Z  RTNEPH,  the  00Z  synthesized  GSM  cloud  analysis  should  closely 
resemble  the  00Z  RTNEPH  cloud  analysis,  at  least  in  RAOB  sparse  areas. 

But  this  expectation  is  not  borne  out  by  the  distressingly  low  0-hour  20/20  score  for  the 
GSM/A  scheme  or  by  subjective  comparison  of  the  RTNEPH  and  0-hour  GSM/A  cloud  fields. 
This  0-hour  20/20  score  penalty  in  all  the  GSM  cloud  schemes  reflects  the  severe  penalty 
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incurred  from  the  totality  of  the  various  vertical  interpolation  and  horizontal  spectral 
transformation  steps  of  the  GSM  pre-  and  postprocessing  circuit,  which  is  schematically 
illustrated  in  Figure  38,  summarized  in  Section  3.1.1  (for  postprocessing  only)  and  discussed  in 
lull  detail  by  Mitchell  and  Yang. 14 
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Figure  38.  Schematic  Illustration  of  the  Major  GSM  Static 
Preprocessing  and  Postprocessing  Steps,  which  when  Applied 
(Along  with  the  Chosen  Cloud  Scheme  as  the  Last  Step)  to  an 
Input  Humidity  Analysis  on  Mandatory  Pressure  Levels  Yield 
the  So-called  "Synthesized"  GSM  O-hour  Cloud  Forecast. 
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Some  of  the  reduction  of  the  O-hour  GSM/A  20/20  score  Is  also  due  to  some  correction  of 
the  RTNEPH/5LAYER-based  CPS  first  guess  (with  its  upper-level  dry  bias)  by  the  actual  RAOBs 
of  humidity.  The  0-hour  GSM/A  mean  NH  octagon  cloud  amount  of  58.8  percent  versus  the 
lower  verifying  cloud  amount  of  48  percent  surely  reflects  some  of  the  RAOB  correction, 
especially  at  the  upper  levels.  Nevertheless,  we  contend  that  a  major  cause  of  the  low  0-hour 
20/20  scores  for  all  GSM  cloud  schemes  In  Tables  12  and  13  is  the  accumulated  static  error 
from  the  GSM  pre-  and  postprocessing  steps. 

Should  AFGWC  In  the  future  embark  on  an  attempt  to  include  pseudo  or  inferred  moisture 
observations  from  the  RTNEPH  In  Its  present-day  HIRAS  global  moisture  analysis  for  the 
purpose  of  Improving  GSM  short-range  cloud  forecasts,  it  must  drastically  reduce  the  steps  and 
substantially  Increase  the  horizontal  resolution  in  the  pre-  and  postprocessing  of  the  GSM 
moisture  fields.  Details  of  one  such  streamlining  in  postprocessing  were  already  given  in 
Section  3.1.1.  In  streamlining  preprocessing,  an  analog  of  Figure  7  must  be  followed  as  closely 
as  possible,  that  is,  (1)  vertically  and  horizontally  compact  the  RTNEPH  to  the  GSM  Internal 
grid  (GSM  sigma  layers  and  Gaussian  grid  —  where  the  latter  grid  and  corresponding  spectral 
resolution  are  carried  at  sufficiently  high  resolution  to  retain  1 /2-mesh  initial  detail)  and  2) 
directly  convert  this  GSM-compacted  RTNEPH  cloud  analysis  to  humidity  using  the  humidity- 
to-cloud  curves  tuned  specifically  for  the  asymptotic  GSM-preferred  RH  forecast  distribution 
(like  the  GL  scheme  curves  of  Section  3.1.3). 

As  an  alternative  (not  applied  In  this  study),  the  so-called  "trending"  approach 
demonstrated  in  the  study  of  MW  can  be  used  with  substantial  success  to  Increase  the  short- 
range  GSM  cloud  forecast  accuracy.  In  the  trending  approach,  as  applied,  for  example,  to 
obtain  a  12-hour  GSM  cloud  forecast,  a  12-hour  GSM  cloud  change  field  would  be  derived  at 
each  of  the  six  moist  GSM  pressure  levels  by  subtracting  the  synthesized  0-hour  Initial  layer 
cloud  field  from  the  12-hour  forecast  layer  cloud  field.  Then  these  GSM  layer  cloud  trends 
(which  isolate  trends  from  the  forward  time  Integration  of  the  GSM)  would  be  added  to  the 
original  compacted  RTNEPH  layer  cloud  fields,  followed  by  the  stacking  algorithm  to  obtain 
total  cloud. 

Although  usable  with  any  diagnostic  cloud  scheme,  the  trending  approach  is  especially 
suited  to  a  cloud  scheme,  like  the  GL  schemes,  having  little  or  no  short-term  spin-up  in 
forecast  mean  cloud  amounts.  The  trending  approach  In  the  MW  study,  applied  in  conjunction 
with  the  GSM/A  or  CPS  scheme,  yielded  20/20  accuracy  scores  surpassing  those  of  5LAYER  at 
12  hours.  The  trending  approach  in  conjunction  with  the  GL4  cloud  scheme  represents  in  our 
view  the  best  possible  approach  to  short-range  GSM  cloud  forecasting.  The  lack  of  tests  of  the 
trending, approach  here  Is  probably  the  present  study's  most  serious  shortfall. 

We  next  briefly  examine  displays  of  the  Case  I  48-hour  GSM  cloud  forecasts  from  the 
GSM/N  scheme  (Figure  23)  In  Figure  39  (total  cloud)  and  Figure  40  (layer  cloud).  Tables  9,  10A, 
and  1 1A  and  Figure  39  show  that,  like  the  GSM/A  scheme,  this  second,  temporally  fixed 
scheme,  also  suffers  a  dramatic  spin-up  to  total  cloud  forecast  fields  having  large  positive 
biases.  In  the  GSM/N  scheme,  this  positive  bias  occurs  not  from  low  critical  values  of  RHc  for 
the  onset  of  nonzero  cloud,  but  rather  from  the  low  critical  values  of  RHm  for  the  mostly 
cloudy  and  overcast  thresholds,  namely  critical  values  RHm  of  75  and  65  percent  for  the  onset 
of  overcast  at  the  70  kPa  and  30-50  kPa  levels,  respectively.  Thus,  the  GSM/N  cloud  scheme 
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specifies  mostly  cloudy  or  overcast  at  RH  values  well  below  100  percent  at  mid  and  upper 
levels.  Figure  40,  for  example,  shows  widespread  areas  of  virtual  overcast  at  the  70  kPa  and  30 
kPa  levels  (similarly  for  the  50  and  40  kPs  levels  not  shown). 
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Figure  40,  For  the  NMC  Cloud  Scheme,  the  GSM  48-hour  Forecast 
of  NH  Laver  Cloud  Amount  Valid  at  00Z  on  19  January  1985  for 
Layers  Nominally  Centered  at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa 
Pressure  Levels. 
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Finally,  the  small  RH  spread  or  interval  of  only  20  percent  between  the  RHc  and  RHm 
values  (clear  and  overcast  RH  thresholds)  of  the  GSM/N  scheme  (see  Section  3.1.2)  for  the  30, 
40,  and  50  kPa  levels  yields  upper-level  layer  cloud  patterns  in  Figure  40  that  are  overly 
"sharp".  That  is,  the  upper-level  cloud  patterns  are  too  binary  in  character,  dominated  by 
large  pervasive  areas  of  strictly  clear  or  strictly  overcast. 

Turning  now  to  the  final  previous  diagnostic  cloud  scheme,  which  is  also  a  temporally- 
fixed  scheme,  we  consider  displays  of  Case  I  cloud  forecasts  from  the  GSM/E  scheme 
(Figure  20)  in  Figure  41  (total  cloud)  and  Figure  42  (layer  cloud).  A  review  of  Figure  21  reminds 
us  that  this  scheme  has  the  largest  low-  and  upper-level  critical  values  of  RH  for  the  onset  of 
nonzero  cloud.  Figure  42  correspondingly  shows  the  small  amounts  of  cloud  yielded  by  this 
scheme  at  the  100  kPa  and  30  kPa  levels.  The  large  upper-level  critical  values  of  the  GSM/E 
scheme  are  especially  ill-suited  for  the  initial  upper-level  AFGWC  humidity  analyses,  which 
are  erroneously  very  dry  over  most  of  the  octagon.  Correspondingly  then,  the  0-hour  total 
cloud  fields  from  the  GSM/E  scheme  show  very  large  negative  biases.  As  the  low-  and  upper- 
level  moistening  trends  set  in  during  the  GSM  forecast,  this  initial  negative  bias  decreases  to  a 
range  of  -10  to  -20  percent  at  48  hours  in  Tables  10-11,  a  significantly  smaller  bias  magnitude 
than  the  20-35  percent  positive  bias  at  48  hours  for  the  GSM/A  and  GSM/N  schemes.  More 
importantly,  the  GSM/E  scheme  showed  higher  20/20  scores  than  the  GSM/A  and  GSM/N 
schemes  in  more  than  half  the  areas  at  36  and  48  hours  for  all  three  cases. 
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Figure  42.  For  the  ECMWF  Cloud  Scheme,  the  GSM  48-hour  Forecast 
of  NH  Laver  Cloud  Amount  Valid  at  00Z  on  19  January  1985  for  Layers 
Nominally  Centered  at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa  Pressure 
Levels. 
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The  latter  somewhat  encouraging  results  with  the  GSM/E  scheme  led  to  numerous  trial  and 
error  experiments  at  subjectively  tuning  the  critical  RH  values,  RHc,  and  nonlinear  exponent, 
Pk  (see  Eqs.  1  and  6).  The  most  promising  re  .  ’It  was  obtained  with  no  change  In  the  exponent, 

and  modest  changes  In  the  RHc  critica.  values  as  given  in  Table  7.  The  given  changes  In  these 
values  clearly  are  not  radical;  yet  experiments  with  more  radical  changes  yielded  poorer  20/20 
scores,  albeit  not  necessarily  poorer  bias  scores.  The  Impacts  of  the  RHc  changes  In  Table  7  on 
the  GSM/E  bias  and  20/20  scores  are  evident  In  Tables  9-13  under  the  designation  GSM/EM 
(for  modified).  The  Improvement  of  GSM/EM  over  GSM/E  scores  is  not  notable. 

3.2.3  THE  NEW  GL  DIAGNOSTIC  CLOUD  SCHEME 

The  disappointing  performance  of  the  previous  diagnostic  cloud  schemes  led  to  a  resolute 
search  for  a  new  diagnostic  cloud  scheme  specifically  tuned  to  the  present  GSM  forecast  model 
and  verifying  RTNEPH  cloud  analysis.  The  earlier  results  showed  It  was  important  to  find  a 
scheme  with  temporal  degrees  of  freedom  to  compensate  for  the  early  GSM  spin-up  period.  The 
early  attempts  to  modify  the  GSM/E  scheme  for  this  purpose  via  trial  and  error  tests  proved 
almost  endless  owing  to  the  number  of  degrees  of  freedom  in  Eq.  (6).  Hence,  it  became  clear 
that  an  objective  approach  was  necessary.  The  routine  gridded  objective  cloud  analyses  of  the 
RTNEPH  provide  the  crucial  discretized  3-D  cloud  database  needed  for  such  an  effort. 

Section  3. 1 .3  described  the  development  and  procedures  of  the  new  GL  objective  cloud-scheme 
methodology,  yielding  first  the  GSM/GL6  cloud  scheme,  which  utilizes  all  six  moist  pressure 
levels  output  by  the  GSM. 

The  cloud  forecasts  In  Case  I  from  the  GSM/GL6  scheme  (Figure  27)  are  displayed  for  the 
0-  and  48-hour  forecast  times  in  Figure  43  (total  cloud)  and  Figures  44-45  (layer  cloud).  (We 
must  note  that  all  GSM/GL  forecasts  displayed  here  In  Figures  43-47  and  listed  In  Tables  9-13 
used  the  so-called  "aggregate"  curves  -  see  discussion  of  Figure  25  -  derived  from  the  combined 
Case  I  and  II  cumulative  frequency  distributions  of  RTNEPH  clouds  and  GSM  RH  forecasts.) 

The  0-hour  and  48-hour  GSM/GL6  total  cloud  fields  in  Figure  43  appear  at  first  glance  to  be 
rather  cloudy,  especially  In  the  tropics  outside  the  octagon.  However,  we  recall  In  Section  3.1.3 
that  the  GL6  cloud  curves  In  Figure  27  were  derived  by  mapping  RTNEPH  cloud  and  GSM  RH 
layer  cumulative  frequency  distributions  derived  only  from  grid  points  Inside  the  octagon  (the 
RTNEPH  outside  the  octagon  was  not  available).  Within  the  octagon.  Table  9  shows  that  the 
CL-hour  GL6  mean  cloud  amount  (50.9)  agrees  closely  with  that  of  the  RTNEPH  (48.0). 
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Figure  43.  For  the  GL6  Cloud  Scheme,  the  GSM  (a)  O-hour  and 
(b)  48-hour  Forecast  of  NH  Total  Cloud  Amount  Valid  at  00Z  on 
17  and  19  January  1985,  Respectively. 


129 


Figure  44.  For  the  GL6  Cloud  Scheme,  the  GSM  O-hour  Forecast 
of  NH  Laver  Cloud  Amount  Valid  at  00Z  on  17  January  1985  for 
Layers  Nominally  Centered  at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa 
Pressure  Levels. 


Figure  45.  For  the  GL6  Cloud  Scheme,  the  GSM  48-hour  Forecast 
of  NH  Laver  Cloud  Amount  Valid  at  00Z  on  19  Januaiy  1985  for 
Layers  Nominally  Centered  at  the  (a)  100,  (b)  85,  (c)  70  and  (d)  30  kPa 
Pressure  Levels. 


Figure  47.  For  the  GL4  Cloud  Scheme,  the  GSM  48-hour  Forecast 
of  NH  Laver  Cloud  Amount  Valid  at  00Z  on  19  January  1985  for 
Layers  Nominally  Centered  at  the  (a)  Low  and  (b)  High  Cloud. 
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This  points  out  that  in  routine  global  appli'  attmis  of  the  GL6  cloud  scheme,  the  curves 
should  probably  be  stratified  by  region.  Surh  a  stratification  might  include  tropical, 
rnidlatitude,  and  polar  domains,  each  further  divided  according  to  water  and  land  points. 
Access  to  Just  two  independent  cases  here  provided  too  small  a  sample  to  test  areal 
stratification  of  the  curves  in  this  study. 

While  the  octagon-mean  lota!  cloud  amounts  of  tire  O-hour  G1J5  and  RTNEPH  agree  closely, 
the  GL6  cloud  field  lacks  the  "rag del  dec!  appearance  of  the  RTNEPH  analysis  or  5 LAYER 
forecast  fields,  that  is,  it  fails  to  pieser-e  the  perceived  distinct  fractal  geometry  of  the 
relatively  high  resolution,  albeit  comported,  RTNEPH.  Tills  failure  again  is  due  to  the  greater 
pre-  and  postprocessing  penalty  or  synthesis  penalty  encountered  in  the  GSM  versus  SLAYER 
frameworks.  The  O-hour  20/20  cores  for  all  the  GSM/GL  schemes  in  Tables  12-13  suffer  this 
very  substantial  synthesis  pennlb  though  not  as  severely  as  the  previous  schemes.  The 
failure  of  the  GSM  synthesis  to  restive  the  detail  of  definition  in  the  initial  RTNEPH  is  also 
recognized  when  comparing  the  O-hour  GL6  layer  cloud  fields  in  Figure  44  with  those  of  the 
RTNEPH  in  Figure  31. 

Tun.jng  to  a  new  matter  ".e  next  consider  the  O-hour  GL6  and  RTNEPH  30  kPa  cloud  fields 
in  Figures  44d  and  3 Id,  which  illustrate  inherent  limits  to  the  ability  of  the  GL  cloud  scheme 
to  overcome  forecast  or  analysis  systematic  position  error.  From  the  AFGWC  initial  RH 
analysis  at  30  kPa  in  Figure  171,  we  recall  that  the  "boundary-area"  blending  of  the  erroneous, 
very  diy,  5LAYER  first  guess  dependent  moisture  analysis  with  the  independent  and  more 
realistic  tropical  moisture  analysis  (see  Table  3)  yields  an  unnatural  ring  of  maximum 
humidity  (with  respect  to  the  octagon  interior)  along  the  inside  of  the  octagon  boundary.  The 
frequency  distribution  of  RH  derived  from  the  octagon  interior  of  Figure  17f  will  contain  no 
"memory"  of  the  predominant  location  of  this  high  humidity.  Most  importantly,  the  position 
of  this  ring  of  relative  maximum  of  humidity  is  an  analysis  artifact:  and  it  does  not  correlate 
well  with  the  location  of  cloudy  aieas  in  the  RTNEPH  analysis  of  Figure  31d,  because  the  ring 
coin:  ides  with  the  generally  <  lend -tree  subsiding  subtropical  trade  wind  zone. 

Hence,  when  the  underlying  diagnostic  cloud-scheme  assumption  of  high  correlation 
between  locations  of  high  humidity  and  large  cloud  amount  breaks  down  due  to  gross 
mlslocations  of  relative  maximums  of  RH  in  the  analysis  or  forecast,  the  humidity-interred 
location  of  the  clouds  will  break  down  also.  This  limitation  is  true  of  any  type-1  humldity-to- 
cloud  scheme  and  is  thus  a  limitation  not  overcome  by  tile  particular  derivation  method  for 
the  GL  schemes.  However,  the  real  culprit  in  the  wrong  location  of  cloudiness  in  Figure  44d  is 
not  tlie  cloud  scheme,  but  rather  the  total  breakdown  of  the  initial  AFGWC  humidity  analysis 
at  tills  level  in  Figure  1  7f.  We  nevertheless  have  duly  noted  that  whenever  the  distribution  of 
relative  maximums  in  humidity  assumes  an  unnatural  character  in  the  GSM  initial  state  or 
forecast,  then  the  GL  cloud  schemes,  like  any  type-1  scheme,  will  also  reflect  that  unnatural 
character. 

What  may  be  more  significant  here  is  that  even  in  Figure  44,  the  GL6  scheme  does  continue 
to  reproduce  the  correct  octagon -mean  30  kPa  cloud  amount  and  the  correct  octagon-domain 
frequency  distribution  of  the  range  of  cloud  amounts  as  provided  by  the  RTNEPH  analysis  in 
Figure  3 Id.  Furthermore,  by  48- hours,  when  the  GSM  forecast  has  largely  corrected  the 


erroneous  dry  bias  of  the  initial  state,  the  GSM  30  kPa  48-hour  cloud  forecast  and  verifying 
RTNEPH  cloud  analysis  now  appear  much  more  similar  (Figures  45d  ami  32d). 

A  similar  problem  of  this  type  is  observed  again  to  a  lesser,  but  nonetheless  detrimental, 
degree  in  the  GSM  RH  forecast  and  GL6  diagnosed  cloud  amount  at  100  kPa  in  Figure  45a. 

Here,  the  GSM's  large  positive  RH  foreca  t  bias  over  subtropical  oceans  (but  still  inside  the 
octagon)  has  led  to  a  predominance  of  diagnosed  cloud  cover  over  the  subtropical  oceans  at  the 
expense  of  mldlatitude  oceanic  cloud,  almost  opposite  the  oceanic  latitudinal  cloud 
distribution  of  the  RTNEPH  In  Figure  32.  Thus  while  the  GL6  scheme  ensures  the  proper 
octagon-mean  cloud  amount  at  every  layer,  the  GSM  systematic  error  in  the  location  of  high 
RH  has  resulted  in  layer  cloud  at  the  100  kPa  level  in  the  wrong  location. 

Consequently.  In  the  vertical  stacking  steps  of  Appendix  A,  the  derived  total  cloud  cover  in 
the  GL6  scheme  Is  somewhat  positively  biased,  despite  the  enforced  zero  layer  biases,  because 
the  GSM  Initial  and  forecast  cloud  layers  have  an  unobserved  vertical  alignment.  Tables  9, 
10A,  and  1 1A  show  the  tendency  of  the  GL6  total  cloud  analysis  and  the  forecast  to  be 
positively  biased.  The  culprit  early  in  the  forecast  is  the  aforementioned  mislocation  of  high 
RH  In  the  30  kPa  analysis,  and  later  in  the  forecast  the  mislocation  of  high  RH  in  the  100  kPa 
forecast.  Both  the  30  kPa  and  100  kPa  problems  could  be  alleviated  substantially,  we 
speculate,  by  a  stratification  of  the  GL  curves  by  area  as  suggested  earlier. 

Unable  to  pursue  areal  stratification  here  owing  to  the  two-case  limit,  to  overcome  the  100 
kPa  problem  we  tested  various  alternative  vertical  stacking  strategies  in  schemes  GL6A,  GL6B, 
GL6C,  GL4.  and  GL4A.  Scheme  GL6B  simply  imposes  a  forecast  of  zero  cloud  at  100  kPa.  As 
expected,  the  total  cloud  bias  of  the  GL6B  scheme  is  negative.  The  20/20  accuracy  score  of  the 
GL6B  scheme  is  mixed,  occasionally  showing  improvement.  Next,  the  GL6C  scheme  also 
ignores  the  100  kPa  GSM  cloud  forecast,  but  additionally  assumes  random  cloud-layer  overlap 
(R=l  in  stacking  Combinations  I,  II.  and  III  of  Appendix  A).  As  a  random  stacking  reference 
case,  the  GL6A  scheme  also  assume^  random  overlap  in  stacking,  but  retains  the  100  kPa  GSM 
cloud  forecast.  Predictably,  the  GL6A  scheme  shows  a  definite  positive  bias  over  the  octagon. 
Though  mixed,  the  20/20  scores  of  these  GL6  scheme  variations  showed  occasional  superiority 
over  the  original  GL6  scheme. 

Thus  in  a  final  attempt  to  circumvent  the  poor  GSM  100  kPa  cloud  forecasts  over  low- 
latitude  oceanic  areas,  the  four-layer  GL4  scheme  was  developed  (Figure  28)  and  tested 
(Figures  46-47).  As  described  in  detail  in  Section  3.1.3,  this  scheme  also  ignores  the  GSM  100 
kPa  RH  forecast  and  handles  the  30  and  40  kPa  GSM  RH  forecasts  differently  in  deriving 
frequency  distributions  of  forecast  RH  at  upper  levels  (to  better  deal  with  the  lack  of  compacted 
RTNEPH  layer  cloud  at  40  kPa).  Lastly,  the  GL4A  scheme  once  again  assumes  random  layer 
overlap  in  the  derivation  of  total  Joud.  For  overall  consistency  in  terms  of  small  steady  bias 
and  competitive  20/20  scores  in  Tables  9-13,  the  GL4  scheme  appears  to  be  the  best  choice 
among  all  GL  cloud  schemes,  virtually  all  of  which  clearly  surpass  the  previous  GSM 
diagnostic  cloud  schemes. 

From  the  results  of  this  section  on  GL  cloud  schemes,  <t  is  clear  that  the  objective  tuning 
procedure  and  temporal  degree  of  freedom  embraced  in  the  GL  scheme  do  indeed  achieve  the 
desired  goals  of  a  small  forecast  bias  and  negligible  spin-up  in  large-area  mean  cloud 
amounts.  These  strengths  do  not  overcome  the  penalties  in  forecast  accuracy  incurred  by  a 
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poor  initial  humidity  analysis  and  the  complex,  multiple,  and  relatively  low  resolution  of  the 
numerous  static  steps  of  GSM  pre-  and  postprocessing.  These  penalties  are  so  severe  that  the 
20/20  accuracy  scores  of  both  the  GL  schemes  and  previous  schemes  actually  increase  over  12 
hours,  often  over  24  hours  or  longer,  indicatinf  that  the  full  potential  for  GSM  forecast 
accuracy  in  the  12-24  hour  range  has  not  been  achieved.  By  36  hours  and  beyond,  the  present 
GSM/GL  schemes  regularly  surpass  the  cloud  forecast  accuracy  of  5LAYER,  sometimes  by  a 
significant  margin.  Recommended  approaches  for  achieving  GSM  cloud  forecast  superiority 
over  5LAYER  at  shorter  forecast  ranges  are  reiterated  in  the  following  conclusion  section. 


4.0  SUMMARY  AND  CONCLUSIONS 


Automated  objective  cloud  forecasts  (especially  in  the  0-24  hour  range)  have  been 
important  products  of  AFGWC  for  at  least  two  decades,29  in  support  of  various  world-wide  Air 
Force  missions.  During  this  period,  AFGWC's  primary  cloud  forecast  model  has  been  the 
5LAYER  model,  and  its  direct  predecessors.  This  family  of  cloud  models  has  always  required 
input  wind  forecasts  from  a  prior  execution  of  a  separate  dynamic  NWP  model,  which  prior  to 
1985  at  AFGWC  were  dry  models  (non-moisture  bearing).  In  October  1985,  AFGWC 
implemented  the  AWAPS,  including  a  moist  global  spectral  model  from  NMC,8  which  for  the 
first  time  afforded  AFGWC  the  possibility  of  producing  global  cloud  forecasts  directly  from  a 
single,  global,  moist  dynamic  model.  In  response  to  this  opportunity  and  at  the  request  of 
AWS,  GL  initiated  a  program  to  develop  cloud  forecast  schemes  suitable  for  use  in  AFGWC's 
present  and  future  global  (and  regional)  model. 

The  purpose  of  this  study  was  to  develop  and  test  a  skillful  scheme  for  forecasting  layer 
and  total  fractional  cloud  cover  in  AFGWC's  presently  operational  GSM  (also  GL's  baseline 
GSM)  anA  to  compare  these  GSM  cloud  forecasts  with  the  corresponding  forecasts  from 
AFGWC's  traditional  5LAYER  model.  One  specific  objective  was  to  determine  the  shortest 
forecast  length,  T,  if  any,  at  which  the  cloud  forecast  skill  of  the  GSM  equalled  or  surpassed 
that  of  5 LAYER.  Thus  this  study  had  to  be  particularly  cognizant  of  the  preeminent  source  of 
5LAYER's  rather  good  short-range  cloud  forecast  skill  —  namely  5 LAYER'S  timely  and  direct 
use  (via  a  minimum  of  static  preprocessing  steps)  of  AFGWC's  RTNEPH  global  3-D  high- 
resolution  cloud  analysis  to  empirically  derive  initial  CPS  (5LAYER’s  moisture  variable). 

Thus,  to  minimize  the  time  T  cited  above,  the  central  GSM  forecasts  in  this  study  were 
purposely  executed  from  AFGWC  global  moisture  analyses  that  also  heavily  utilized  RTNEPH 
clouds  to  infer  the  first-guess  fields  for  the  moisture  analysis. 

As  a  precursor  to  the  GSM  cloud  forecasts,  we  first  examined  these  AFGWC  relative 
humidity  (RH)  analysis  fields,  along  with  the  follow-on  GSM  RH  forecasts.  We  established  that 
the  above  RTNEPH- Influenced  AFGWC  moisture  analyses  were  much  too  dry  in  the  upper  half 
of  the  troposphere  owing  to  (1)  overly  low  zero-cloud  humidity  thresholds  in  the  clcud-to- 
humidlty  inference  scheme  (that  is.  the  5LAYER  CPS-cloud  scheme)  and  (2)  underestimation  of 
high  cloud  amounts  in  the  RTNEPH.  Unlike  the  5 LAYER  model,  which  maintains  this 
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substantial,  initial  upper-level  dry  bias  during  its  forecast  (and  therefore  still  yields 
reasonable  forecast  mean  cloud  amounts  when  it  applies  the  CPS-to-cloud  scheme  at  forecast 
time),  the  GSM  moistens  these  upper  levels  during  its  forecast  and  largely  corrects  the 
erroneous  initial  dry  bias.  The  large  magnitude  of  this  upper-level  GSM  spin-up,  which  is  an 
artifact  of  the  poor  AFGWC  upper-level  moisture  analysis,  severely  degrades  the  GSM  short- 
range  cloud  forecast  skill. 

Separate  GSM  RH  forecasts  from  more  realistic  NMC  moisture  analyses  demonstrated 
much  less  spin-up.  Additionally,  these  latter  GSM  forecasts,  together  with  the  GSM  forecasts 
from  the  AFGWC  analyses,  showed  that  the  GSM  in  24  hours  evolves  quickly  toward  its  own 
model-preferred  zonal  mean  humidity  state  and  model-preferred  RH  frequency  distribution, 
largely  irrespective  of  the  chosen  initial  moisture  analysis.  As  a  class,  large-scale  moist  NWP 
models  characteristically  show  distinct  drift  toward  a  model-preferred  moisture  state  and  a 
concomitant  model-preferred  hydrological  balance.  However,  these  model-preferred  moisture 
states  differ  from  observations  and  from  model  to  model  to  various  degrees  depending  on  the 
veracity  of  the  model's  physical  parameterizations.  The  GSM  preferred  moisture  state  here 
exhibits  a  small  upper-level  moist  bias  and  a  modest  mid-  and  lower-level  dry  bias,  except  for 
a  substantial  lowest-level  moist  bias  over  tropical  (warm)  oceans. 

To  produce  cloud  forecasts  with  the  GSM,  we  sought  an  appropriate  diagnostic  or  empirical 
scheme  to  infer  cloud  cover  from  standard  model  output.  Because  the  baseline  GSM  here 
substantially  underestimates  convective  precipitation  rate  and  lacks  parameterized  physics 
for  (1)  radiation,  (2)  a  surface  energy  budget,  and  (3)  land  surface  and  PBL  moisture  fluxes,  we 
considered  only  the  simplest  diagnostic  cloud  schemes  (here  called  type-1  schemes),  which 
diagnose  fractional  cloud  cover  using  only  the  model's  RH  forecast.  In  these  simplest 
diagnostic  cloud  schemes,  the  conceptual  physical  picture  is  as  follows:  On  spatial  resolutions 
typical  of  a  GSM,  clouds  are  frequently  subgrid  scale,  both  horizontally  and  vertically.  Hence, 
although  the  forecast  volume-average  RH  for  a  given  grid  box  volume  may  be  somewhere  below 
100  percent,  we  may  expect  some  fractional  cloud  cover  in  the  grid  box  since  the  local  variance 
about  the  average  will  yield  subgrid  regions  of  saturation  (and  hence  clouds). 

In  type-1  diagnostic  cloud  schemes,  the  crucial  empirical  parameter  is  the  critical 
threshold  of  humidity,  or  critical  humidity  (RHc),  at  which  nonzero  fractional  cloud  cover  is 
first  Inferred.  In  more  complex  diagnostic  cloud  schemes  (type-2),  which  we  recommend  be 
tested  in  follow-on  studies  when  the  GSM  has  improved  convection,  radiation,  and  PBL 
physics,  other  model  output  fields  such  as  convective  precipitation  rate,  vertical  velocity,  and 
static  stability  should  be  utilized,  since  the  magnitude  of  the  humidity  variance  within  a  grid 
volume  is  undoubtedly  a  function  of  these  variables  (for  example,  stable  stratification  or 
descending  motion  will  likely  be  associated  with  a  decrease  in  subgrid  moisture  variance  and 
hence  imply  an  Increase  in  the  RHc  parameter). 

In  the  GSM,  we  first  tested  several  pre-existing  type- 1  diagnostic  cloud  schemes  (for 
example,  the  scheme  used  by  ECMWF  prior  to  May  1985).  It  was  found  that  because  type-1 
cloud  schemes  from  another  model  usually  use  RHc  values  tuned  for  that  model's  preferred 
moisture  state,  such  schemes  perform  poorly  when  applied  to  an  independent  model  and  yield 
large  cloud  forecast  biases  and  unacceptable  spin-up  behavior. 
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Thereafter,  we  developed  a  new  methodology,  herein  called  the  GL  cloud  scheme,  to 
objectively  derive  a  type- 1  diagnostic  cloud  scheme  tuned  specifically  for  any  given  NWP 
model.  The  methodology  provides  a  cloud  scheme  that  manifests  (1)  negligible  bias  and  (2) 
negligible  spin-up  in  the  forecast  cloud  amounts.  The  method  derives  explicit  RH-to-cloud 
conversion  functions  or  "curves"  layer  by  layer.  The  method  does  so  by  calculating  the 
cumulative  frequency  distribution  of  the  model's  RH  forecasts  for  a  given  forecast  length  at  a 
given  layer  (over  a  sample  of  forecast  cases  spanning  several  weeks)  and  mapping  this 
distribution  onto  the  cumulative  frequency  distribution  oi  RTNEPH  observed  cloud  amounts  at 
that  layer  for  the  same  time  period.  By  mapping  the  model's  RH  frequency  distribution  to 
observed  cloud  amounts  at  each  regular  interval  of  model  forecast  length  (say  6-hourly),  one 
can  account  for  the  model's  changing  RH  frequency  distribution  during  the  spin-up  period. 

The  resulting  forecast-time  dependent  RH-to-cloud  curves  eliminate  spin-up  biases  in  the 
derived  forecast  cloud  amounts. 

Because  the  cumulative  frequency  distribution  of  model  forecast  RH  or  RTNEPH  cloud 
amounts  can  vary  substantially  over  various  climatological  regions  (owing  to  natural  factors, 
failures  in  model  physics,  or  breakdowns  in  the  RTNEPH  cloud  retrieval  algorithms),  we 
recommend  that  the  above  methodology  be  carried  out  separately  over  large  climatologlcally 
homogeneous  regions  (say  tropical,  midlatitude,  and  polar  areas,  each  stratified  by  land  and 
water).  To  account  for  the  seasonal  cycle  in  these  areas,  we  recommend  a  real-time  update  of 
the  respective  RH  and  cloud  frequency  distributions  based  on  a  continuously  updated  running 
sample  over  a  period  of  several  weeks,  coupled  to  a  corresponding  update  of  the  RH-to-cloud 
curves.  Such  a  continuous  update  would  also  account  for  changing  RH  or  cloud  amount 
distributions  stemming  from  ongoing  developmental  changes  in  model  resolution,  numerics, 
and  physics  or  RTNEPH  cloud  retrieval  algorithms. 

Regarding  the  cloud  forecast  performance  of  the  GL  cloud  scheme,  in  two  separate  winter 
cases,  the  GSM  cloud  forecast  accuracy  as  verified  against  the  RTNEPH  showed  the  GL  cloud 
scheme  consistently  outperformed  three  previous  type-1  cloud  schemes,  including  the  pre-1985 
ECMWF  scheme.  However,  even  with  the  GL  cloud  scheme,  the  GSM  performance  relative  to 
5LAYER  was  less  than  hoped,  but  not  surprising  after  investigating  the  impacts  of  the  poor 
AFGWC  initial  humidity  analyses.  Indeed,  the  accuracy  of  the  GSM  cloud  forecasts  with  the  GL 
and  other  type-1  schemes  actually  increases  in  the  first  24  hours,  as  the  GSM  corrects  and 
recovers  from  the  poor  AFGWC  initial  moisture  state.  Hence,  usually  in  the  24-  to  36-hour 
range,  the  GSM  cloud  forecast  skill  with  the  GL  scheme  surpasses  that  of  5LAYER  over  the 
U.S.  and  European  areas. 

In  addition  to  the  poor  AFGWC  initial  moisture  analysis,  another  cause  of  poor  GSM  short- 
range  cloud  forecast  accuracy  is  the  substantial  O-hour  GSM  error  incurred  in  the  numerous, 
complex,  and  relatively  low-resolution  GSM  static  pre-  and  postprocessing  steps.  Section  3.1.1 
describes  straightforward  ways  of  reducing  the  postprocessing  steps  for  cloud  forecast 
purposes.  The  reduction  of  preprocessing  steps  is  less  straightforward  and  more  challenging. 

The  approaches  one  might  take  to  reduce  preprocessing  error  hinge  on  whether  the 
moisture  forecasts  from  the  GSM  must  support  multiple  users  of  several  model  moisture 
outputs  (RH,  precipitation,  etc.)  or  users  only  interested  in  explicit  short  range  (0-24  hour) 
cloud  forecasts.  In  the  latter  rather  unlikely  situation,  the  GSM’s  initial  moisture  field  could 
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be  initialized  directly  on  the  GSM's  internal  grid  (model  sigma  layers  and  Gaussian  grid)  from 
the  RTNEPH’s  cloud  analysis  by  inverting  the  RH-to-cloud  curves  derived  for  the  model’s  long 
term  (96-hour)  preferred  asymptotic  mean  humidity  state.  In  the  more  likely  event  of  multiple 
users  of  various  GSM  moisture  outputs,  one  must  utilize  a  more  conventional  4-D  assimilation 
of  explicit  moisture  observations  from  RAOBs  (possibly  folded  in  with  moisture  estimates 
derived  from  the  RTNEPH  in  the  manner  explored  by  Norquist12,28),  using  an  OI  (or  other) 
objective  analysis  scheme  and  a  GSM  first-guess  moisture  forecast. 

The  latter  approach  would  again  suffer  the  full  short-range  cloud  forecast  penalty  incurred 
by  the  GSM’s  complex  O-hour  preprocessing  suite.  The  only  viable  method  known  to  us  at 
present  for  overcoming  this  O-hour  penalty  from  the  full  GSM  static  preprocessing  is  the 
trending  technique  described  in  Section  3.2.1  and  demonstrated  in  the  forerunner  study  of 
MW6. 

Since  the  archiving  period  of  this  study,  AFGWC  has  vastly  improved  its  global  moisture 
analysis  (in  HIRAS  -  see  Section  2.1.1),  from  the  point  of  view  of  providing  a  suitable  initial 
moisture  state  for  a  GSM  moisture  forecast  intended  for  multiple  users.  This  HIRAS  RH 
analysis  is  now  based  on  a  4-D  assimilation  using  a  GSM  first-guess,  RAOB  humidity 
measurements,  and  an  OI  objective  analysis.  RTNEPH-based  moisture  estimates  are  presently 
ignored.  Given  the  improved  status  of  the  HIRAS  moisture  analysis,  a  follow-on  GL-funded 
study  is  now  underway  at  AFGWC  to  test  the  GL  cloud  scheme  in  the  AFGWC  GSM  environment 
on  a  continuous  real-time  basis  over  a  period  of  about  9  months.  This  follow-on  study  will 
test  all  the  above  recommendations  regarding  the  type-1  GL  cloud  scheme,  including 
continuous  real-time  updating  of  the  RH-to-cloud  curves,  areal  stratification  of  the  curves, 
reduced  postprocessing  steps,  and  the  trending  approach  to  reduce  initial  O-hour  static 
processing  error. 
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PGGE 
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NH 

NMC 

NWP 

OI 

PBL 

PE 

RAOB 

RH 

RTNEPH 
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List  of  Acronyms 


Air  Force  Global  Weather  Central 

Air  Force  Supercomputer  Center-Kirtland  AFB 

Advanced  Weather  Analysis  and  Prediction  System 

Air  Weather  Service 

condensation  pressure  spread 

European  Centre  for  Medium  Range  Weather  Forecasts 

First  GAP.0  r'irih-’i  Experiment 

Global  Atmospheric  Research  Program 

Global  Data  Assimilation  System 

Geophysics  Laboratory 

global  spectral  model 

High  Resolution  Analysis  System 

lnter-troplcal  convergence  zone 

Multi-level  Analysis  Model 

Northern  Hemisphere 

National  Meteorological  Center 

numerical  weather  prediction 

optimum  Interpolation 

planetary  boundary  layer 

primitive  equation 

radiosonde  observation 

relative  humidity 

Real-time  Nephanalysis  Model 


146 


> 


SH 

TROPUA 

USAF 

3DNEPH 

5LAYER 


Southern  Hemisphere 
Tropical  Upper-air  Analysis  Model 
United  States  Air  Force 
Three-dimensional  Nephanalysls  Model 
Five-layer  Cloud  Forecast  Model 
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APPENDIX  A 

The  Vertical  Cloud  Stacking  Algorithm 

To  obtain  a  forecast  of  total  cloud  from  a  forecast  of  cloud  cover  at  several  layers,  a 
vertical  cloud  stacking  algorithm  must  be  used.  Cloud  stacking  algorithms  usually  fall  Into 
one  of  two  categories,  adhering  either  to  the  assumption  of  (1)  maximum  layer  overlap  or 
(2)  random  layer  overlap.  In  the  more  general  algorithm  applied  here,  taken  from  that  used  In 
the  AFGWC  5LAYER  model,  a  free  parameter  referred  to  as  the  stacking  factor  R  (detailed 
shortly)  is  Introduced  to  permit  maximum  overlap,  random  overlap,  or  any  degree  in  between. 

The  various  cloud  forecast  schemes  and  models  addressed  in  this  study  yield  layer  cloud 
forecasts  at  six  layers  (In  the  schemes  GSM/A,  GSM/N,  GSM/E,  GSM/EM,  GSM/GL6),  five 
layers  (in  the  5LAYER  model),  and  four  layers  (in  the  scheme  GSM/GL4).  To  most  easily 
describe  the  application  of  the  stacking  algorithm  to  these  various  numbers  of  layers.  It  is 
,  convenient  to  present  the  scheme  as  a  four-layer  scheme,  which  In  cases  of  more  than  four 

Input  layers  is  preceded  by  a  preliminary  step  that  combines  one  or  more  pairs  of  neighboring 
cloud  layers  using  maximum  overlap.  Table  A1  shows  the  preliminary  step  in  the  six-,  five-, 

T  and  four-layer  cases  encountered  in  this  study.  In  the  four-layer  case  of  scheme  GSM/GL4,  the 

preUminaxy  step  is  really  a  null  or  "do  nothing"  step,  with  the  exception  of  terrain  masking 
described  next. 
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Table  Al.  Preliminary  Step  of  the  Vertical  Cloud  Stacking 
Algorithm. 


GSM /A,  N.  E.  EM,  GL6 

CLOUD  SCHEMES 

5 LAYER  MODEL 

CLOUD  TEST 

GSM/GL4 

CLOUD  SCHEME 

OUTPUT  LAYER 

CLOUD  AMNT 

100  kPa\  MAX  cld 

85  kPa  J  Amnt 

Gra<Uenl'l  MAX  Cld 
85  kPa  J  Amnt 

LOW  Cld 

Amnt 

Layer  1 

70  kPa  Cld 

Amnt 

70  kPa  Cld 

Amnt 

70  kPa  Cld 

Amnt 

Layer  2 

50  kPa  Cld 

Amnt 

50  kPa  Cld 

Amnt 

50  kPa  Cld 

Amnt 

Layer  3 

40  kPa].  MAX  Cld 
85  kPa  J  Amnt 

30  kPa  Cld 

Amnt 

HIGH  Cld 

Amnt 

Layer  4 

During  the  preliminary  step  for  all  columns  outlined  In  Table  Al,  at  grid  points  where  the 
terrain  height  Is  above  one  or  more  of  the  listed  standard  pressure  surfaces,  the  cloud  amount 
is  set  to  zero  at  each  terrain-filled  surface.  For  this  purpose,  the  5 LAYER  1 /2-mesh  field  of 
terrain  heights  Is  used,  for  both  the  5LAYER  and  GSM  cloud  layers.  This  terrain  field  never 
extends  to  40  kPa,  so  at  least  Layer  4  In  Table  Al  Is  an  active  layer  at  every  grid  point. 

Overall,  veiy  few  points  are  terrain-filled  at  70  kPa  and  above.  Finally,  as  noted  in  the  main 
text,  the  "Gradient"  layer  cited  In  the  5LAYER  column  of  Table  Al  denotes  a  terrain-following 
layer  centered  at  6  kPa  above  the  terrain-surface  pressure.  Thus  In  the  5LAYER  case  of 
Table  Al,  whenever  the  85,  70,  or  50  kPa  levels  are  terrain-filled,  the  MAX  operation  with  the 
Gradient  layer  shifts  to  Layer  2,  3,  or  4  respectively  in  Table  Al. 

Alter  the  preliminary  step,  the  stacking  algorithm  operates  on  the  four  Input  layer  cloud 
amounts  according  to  the  three  successive  pairings  or  combinations  shown  below: 


Layer  1 
Layer  2 
Layer  3 
Layer  4 


} 

} 


Combination  I 

Combination  II 


Combination  III  (Gives  Total  Cloud) 


To  perform  each  Combination  I,  II,  or  III,  one  begins  with  two  layer  cloud  amounts,  identified 
generically  as  amounts  A  and  B,  expressed  in  percent,  with  A  denoting  the  largest.  The 
amounts  A  and  B  are  divided  by  100  to  get  fractional  amounts  CA  and  CB  In  the  range  0  to  1. 
By  way  of  background,  If  cloud  presence  In  one  layer  Is  considered  an  event  Independent  of 
cloud  presence  In  the  other  layer  (that  is,  the  vertical  alignment  is  random),  then  the 
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vertically  superimposed  overlap  of  the  two  layers  Is  the  product  CA*CB.  The  total  or  combined 
sky- cover  fraction,  CC,  of  the  two  layers  Is  then  given  by 


CC  =  CA  +  CB  -  CA*CB, 


(Al) 


or  equivalently 


CC  =  CA  +  (1-CA)*CB 


(A2) 


Eq.  (A2)  shows  that  CC  can  be  viewed  as  the  cloud  fraction  from  layer  A  plus  that  part  of  cloud 
layer  B  randomly  aligned  over  the  clear  area  of  layer  A.  Actual  experience  shows,  however, 
that  cloud  presence  in  one  layer  is  often  positively  correlated  with  cloud  presence  in  another 
layer,  the  correlation  being  greater  for  smaller  layer  separation  distances.  To  account  for  this 
observation,  the  stacking  algorithm  applied  In  this  study  generalizes  Eq.  (A2)  by  introducing  a 
"stacking  factor"  R  according  to 


CC  =  CA  +  (1-CA)*CB*R 
CA  >  CB  ,  0  <  R  <  1 


(A3) 


The  two  limiting  cases  of  R=0  and  R=1  yield  the  combined  cloud  amount,  respectively,  for 
maximum  overlap  (giving  the  smallest  CC)  and  random  overlap  (giving  the  largest  CC). 
Intermediate  values  of  R  give  combined  cloud  amounts  falling  between  these  limiting  values. 

It  can  be  rigorously  proved  that  CC  in  Eq.  (A3)  Is  bounded  in  all  cases  by  0  and  1 . 

It  remains  only  to  specify  the  values  of  R  used  In  this  study  for  the  three  combinations  I,  II, 
and  III  diagrammed  above.  Although  the  details  are  somewhat  tedious,  overall  we  allow  for 
larger  R  values  for  increasing  distances  between  cloud  layers.  In  the  "MAX"  operations  listed 
in  the  preliminary  stacking  step  of  Table  Al,  we  implicitly  applied  Eq.  (A3)  with  R=0,  that  Is 
we  assumed  adjacent  layers  were  maximally  overlapped.  In  Combination  I  of  Layers  1  and  2, 
we  use  R=0.2,  with  the  exception  of  R=0.3  If  the  100  kPa  or  Gradient  layer  is  chosen  In  the 
MAX  operation  in  the  preliminary  step  giving  Layer  1  in  Table  Al .  In  Combination  II  of 
Layers  3  and  4,  we  use  R=0.3,  with  the  exception  of  R=0.4  if  the  30  kPa  layer  cloud  amount  is 
chosen  in  the  MAX  operation  in  the  preliminary  step  giving  Layer  4  In  Table  Al.  In 
Combination  in,  which  combines  the  resulting  cloud  amounts  of  I  and  II,  R  takes  one  of  the 
nine  values  given  by  R  =  RM(I,J)  from  the  3X3  RM  matrix  below: 
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In  assigning  R  =  RM(I,J),  index  I  is  set  to  1,  2,  or  3  depending  on  whether  in  arriving  at 
Combination  I  (including  the  preliminary  step  of  Table  Al),  the  largest  original  layer  cloud 
amount  occurs  in  the  Gradient/ 100,  85,  or  70  kPa  layer,  respectively.  The  index  J  is  set  to  1,  2, 
or  3  depending  on  whether  in  arriving  at  Combination  II,  the  largest  original  layer  cloud 
amount  occurs  in  the  30,  40,  or  50  kPa  layer,  respectively.  In  the  5LAYER  forecast  case,  J=2  is 
excluded  since  an  explicit  40  kPa  layer  cloud  forecast  is  not  provided.  In  the  GSM/GL4  cloud 
forecast  case,  J=2  and  f  =1  are  excluded  similarly. 

Generally,  this  stacking  algorithm,  including  the  preliminary  step,  utilizes  R  values  in  the 
low  to  intermediate  range  (usually  0.6  or  less),  thus  it  tends  somewhat  toward  the  maximum 
overlap  assumption.  Clearly,  the  stacking  factor  R  is  a  tuneable  parameter,  especially  since 
total  cloud  in  the  verifying  RTNEPH  cloud  analysis  is  more  reliable  than  layer  cloud.  In  the 
course  of  the  present  study,  stacking  sensitivity  tests  were  executed,  wherein  the  above  baseline 
R  values  were  replaced  with  (1)  R=0  throughout,  (2)  R=1  throughout,  (3)  the  greater  of  baseline  R 
less  0.2  and  zero,  and  (4)  the  greater  of  baseline  R  less  0.4  and  zero.  Only  the  case  of  R=1 
showed  significantly  improved  results,  which  are  provided  in  Section  3.2.  In  future  studies, 
additional  stacking  experiments  are  warranted. 
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